[R] R and Data Storage

2005-09-30 Thread rab45+
Where I work a lot of people end up using Excel spreadsheets for storing data. This has limitations and maybe some less than obvious problems. I'd like to recommend a uniform way for storing and archiving data collected in the department. Most of the data could be stored in simple csv type files

Re: [R] repeated - R package - Compilation Error

2005-08-21 Thread rab45+
rab45 == rab45 rab45 on Thu, 11 Aug 2005 09:25:14 -0400 (EDT) writes: . rab45 I'm not sure what your point is. I'm getting a compilation error for a rab45 package that should compile without errors. The error message doesn't say rab45 anything about

[R] glmmPQL and Convergence

2005-08-19 Thread rab45+
I fit the following model using glmmPQL from MASS: fit.glmmPQL - glmmPQL(ifelse(class==Disease,1,0)~age+x1+x2,random=~1|subject,family=binomial) summary(fit.glmmPQL) The response is paired (pairing denoted by subject), although some subjects only have one response. Also, there is a perfect

[R] Mixed Effects Model Power Calculations

2005-08-16 Thread rab45+
Is there an R package available that would facilitate doing a power/sample size analysis for linear mixed effects models? I have seen the Java applets made available by Russell Length which would seem to be able to handle most any lme, but there is little documentation and it's not clear how the

Re: [R] g77 and R

2005-08-12 Thread rab45+
[EMAIL PROTECTED] wrote: Jerry, Here are the two R packages. Thanks. Rick All the fortran codes compile OK so the problem is invoking the compiler you want. One approach which is a work around is to rename the gfortran executable to somthing else like gfortran-save and create a link

Re: [R] g77 and R

2005-08-12 Thread rab45+
[EMAIL PROTECTED] wrote: Jerry, Here are the two R packages. Thanks. Rick All the fortran codes compile OK so the problem is invoking the compiler you want. One approach which is a work around is to rename the gfortran executable to somthing else like gfortran-save and create a link

Re: [R] repeated - R package - Compilation Error

2005-08-11 Thread rab45+
On Wed, 10 Aug 2005 [EMAIL PROTECTED] wrote: Thompson's Manual to Accompany Agresti's book refers to a package named repeated. It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? See the FAQ, Q5.1.5. As to the `best' place, it

Re: [R] repeated - R package - Compilation Error

2005-08-11 Thread rab45+
On Wed, 10 Aug 2005 [EMAIL PROTECTED] wrote: Thompson's Manual to Accompany Agresti's book refers to a package named repeated. It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? See the FAQ, Q5.1.5. As to the `best' place, it

Re: [R] repeated - R package - Compilation Error

2005-08-11 Thread rab45+
When I go into binnest.f, here is what I see: C Calculate First Derivative Do kk = 1, t1 Gradient(kk) = Gradient(kk) + (D1_beta(kk) / Li) End do Do kk = 1, t2 Gradient(t1+kk) = Gradient(t1+kk) + (D1_Sig1(kk) / Li) End do

Re: [R] repeated - R package - Now rmutil Compilation Errors

2005-08-11 Thread rab45+
In my quest to install the repeated package, I have also to install the rmutil package. (BTW, I'm running Fedora Core 4 and R 2.1.1.) But I now get several compilation errors for rmutil and I'm afraid to try to fix them: [EMAIL PROTECTED] R]# R CMD INSTALL rmutil * Installing *source* package

Re: [R] repeated - R package - Compilation Error

2005-08-11 Thread rab45+
You know what spams are? It is highly recommended to try such basic things and to get the depencies right _before_ posting to mailing lists. Regards Michael [EMAIL PROTECTED] wrote: When I go into binnest.f, here is what I see: C Calculate First Derivative Do kk = 1, t1

Re: [R] repeated - R package - Now rmutil Compilation Errors

2005-08-11 Thread rab45+
Please note that you are using gfortran rather than g77 - they are not the same, and very possibly differ in how they react to code. The 1e301 constants do overflow the maximum real value (of the order of 3e+38), but 1d301 is less than the maximum double value of the order of 1.7e+308. So

[R] repeated - R package

2005-08-10 Thread rab45+
Thompson's Manual to Accompany Agresti's book refers to a package named repeated. It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? Rick B. __ R-help@stat.math.ethz.ch mailing list

[R] system is exactly singular

2005-08-10 Thread rab45+
When trying to fit a generalized linear mixed model using glmmPQL: fit0 - glmmPQL(ifelse(response==A,1,0)~gender,data=set1, random=~1|subject,family=binomial) iteration 1 Error in solve.default(pdMatrix(a, fact = TRUE)) : Lapack routine dgesv: system is exactly singular Could this

Re: [R] system is exactly singular

2005-08-10 Thread rab45+
When trying to fit a generalized linear mixed model using glmmPQL: fit0 - glmmPQL(ifelse(response==A,1,0)~gender,data=set1, random=~1|subject,family=binomial) iteration 1 Error in solve.default(pdMatrix(a, fact = TRUE)) : Lapack routine dgesv: system is exactly singular Could

[R] Adding Functionality to stat.table in Epi

2005-08-09 Thread rab45+
The stat.table function in the Epi package won't do standard deviations. It didn't seem that it would be difficult to add an sd function to the stat.table function. Following the example for the mean, I set up a similar function for the sd (and included it as an options) but it just won't work. (I

Re: [R] Adding Functionality to stat.table in Epi

2005-08-09 Thread rab45+
After you copy stat.table to stat.table2 and modify stat.table2 try: environment(stat.table2) - environment(stat.table) (you should only need to do that 1 time after creating/editing stat.table2). hope this helps, Greg Snow, Ph.D. Statistical Data Center, LDS Hospital Intermountain