[R] Fwd: R package or function for ordinal predictors and continous response

2009-09-03 Thread shirley zhang 
Dear list,

I have a microarray data in which gene expression is the response
(dependent) variable. Can anybody tell me which package/function in R
should I use to model the gene expression with an ordinal predictor
(Independent Variable)

Thanks,
Shirley

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] where to get chr_rpts file for dbSNP human 36.3 assembly

2011-11-02 Thread shirley zhang 
Dear list,

In terms of dbSNP database in NCBI, I can get the chr_rpts files for
the most recent 37.3 assembly from the following FTP site,

ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/chr_rpts/

My question is how/where I can get these chr_rpts files based on the
36.3 assembly

Thanks,
Shirley

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] anova function to test the difference between two coefficients in nlme package

2008-01-22 Thread shirley zhang 
Dear Dr. Bates, and R-help,

I've tried the anova function  to test the difference between two
coefficients, as shown on page 225 of your book "Mixed Effects Models
in S and S-Plus (Statistics and Computing)".

When I type:  anova( fm2BW.lme, L = c(TimeDiet2 = 1, TimeDiet3 = -1) )

I got the following error message:

Error: unexpected '=' in "anova( fm2BW.lme, L = c(Time:Diet2 ="

Could you tell me whether I miss anything?

> sessionInfo()
R version 2.6.0 (2007-10-03)
i386-pc-mingw32

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] stats graphics  grDevices utils datasets  methods   base

other attached packages:
[1] nlme_3.1-85

loaded via a namespace (and not attached):
[1] grid_2.6.0 lattice_0.16-5


Thanks,
Shirley

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] anova function to test the difference between two coefficients in nlme package

2008-01-22 Thread shirley zhang 
Dear Dr. Bates, and R-help,

I just found the reason I got the error message is because I copied
the code directly from the book. When I tried the code in
"~\library\nlme\scripts" under the R installation directory like the
following, I got the same answere as in the book.

anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = -1))

Now my question is instead of test the difference between two
coefficients, can I use anova to test the significance of the sum of
two coefficients, like

anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = 1))

Thanks,
Shirley


On Jan 22, 2008 10:38 AM, shirley zhang <[EMAIL PROTECTED]> wrote:
> Dear Dr. Bates, and R-help,
>
> I've tried the anova function  to test the difference between two
> coefficients, as shown on page 225 of your book "Mixed Effects Models
> in S and S-Plus (Statistics and Computing)".
>
> When I type:  anova( fm2BW.lme, L = c(TimeDiet2 = 1, TimeDiet3 = -1) )
>
> I got the following error message:
>
> Error: unexpected '=' in "anova( fm2BW.lme, L = c(Time:Diet2 ="
>
> Could you tell me whether I miss anything?
>
> > sessionInfo()
> R version 2.6.0 (2007-10-03)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
> States.1252;LC_MONETARY=English_United
> States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
>
> attached base packages:
> [1] stats graphics  grDevices utils datasets  methods   base
>
> other attached packages:
> [1] nlme_3.1-85
>
> loaded via a namespace (and not attached):
> [1] grid_2.6.0 lattice_0.16-5
>
>
> Thanks,
> Shirley
>

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] how to use anova() to test the sum of coefficients in nlme package

2008-01-26 Thread shirley zhang 
Dear R-help,

In nlme package, anova () can be used to test the difference between
two coefficients as shown on page 225
of "Mixed Effects Models in S and S-Plus":

anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = -1))

Now my question is instead of test the difference between two
coefficients, can I use anova to test the significance of the sum of
two coefficients, like

anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = 1))

Thanks,

Shirley

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.