Re: [R] random effects model

2013-01-18 Thread rex2013
those questions in ([hidden > email] <http://user/SendEmail.jtp?type=node&node=4655802&i=0>). I > suggest you to read lme4 book (http://lme4.r-forge.r-project.org/lMMwR/) > #lrgprt.pdf > > A.K. > > > > > > > > - Original Message - > Fr

Re: [R] random effects model

2013-01-16 Thread rex2013
Hi I tried removing the missing values and installing "plyr". Still error message appears with ggplot2 Btw, did you get the attachments with my earlier mail? Ta. On Wed, Jan 16, 2013 at 3:16 AM, arun kirshna [via R] < ml-node+s789695n4655612...@n4.nabble.com> wrote: > > > Hi, > Check these lin

Re: [R] random effects model

2013-01-14 Thread rex2013
Sorry I have corrected the mistakes: BP.stack3 <- reshape(Copy.of.BP_2,idvar="CODEA",timevar="time",sep="_",varying=list(c("Obese14","Obese21"),c("Overweight14","Overweight21")),v.names=c("Obese","Overweight"),times=factor(c(1,2)),direction="long") BP.stack3 head(BP.stack3) tail(BP.stack3) names

Re: [R] random effects model

2013-01-11 Thread rex2013
n(x) > colSums(is.na(x[-1] > > # V1 V2 > #F 1 0 > #M 2 1 > > Regarding the "mi" package, I never used it before. BTW, you didn't > mention what you are going to plot. > Hope it helps. > A.K. > > > > - Original Message - >

Re: [R] random effects model

2013-01-11 Thread rex2013
e, one has to inform the function > about different sep- > arations if there are missing values and other correlation structures than > the independence or > exchangeable structures are used. The waves arguments takes an integer > vector that indicates > that two observations of t

Re: [R] random effects model

2013-01-08 Thread rex2013
0 >> BP.stack7<-BP.stack6[BP.stack6$MaternalAge!=5,] >> >> BP.stack7$MaternalAge<-factor(as.numeric(as.character(BP.stack7$MaternalAge) >> >> fit5<-geese(hibp~MaternalAge*time,id=CODEA,data=BP.stack7,family=binomial,corstr="unstructured",sca

Re: [R] random effects model

2013-01-07 Thread rex2013
<- reshape(seizure, >varying=list(c("base","y1", "y2", "y3", "y4")), >v.names="y", times=0:4, direction="long") > seiz.l <- seiz.l[order(seiz.l$id, seiz.l$ti

Re: [R] random effects model

2013-01-06 Thread rex2013
Hi A.K Regarding my question on comparing normal/ obese/overweight with blood pressure change, I did finally as per the first suggestion of stacking the data and creating a normal category . This only gives me a obese not obese 14, but when I did with the wide format hoping to get a obese14,norm