search for: agegrp

How do you get the design matrix R used when calculating ANOVA? -- View this message in context: http://www.nabble.com/how-to-get-design-matrix--tp23464638p23464638.html Sent from the R help mailing list archive at Nabble.com.

Hi, I am trying to change a SAS macro to R. here is my code. I get an error at the last line. attach(fram) dset1<-cbind(AGE,BMI,DEATH) > BMIGRP<-cut(BMI,breaks=3,right=TRUE) > AGEGRP<-floor(AGE/10)-2 > dset<-cbind(AGEGRP,BMIGRP,DEATH) > maxage<-max(dset[,1]) > minage<-min(dset[,1]) > #maxcls<-dset[,2] > #mincls<-dset[,2] > nage<-maxage-minage+1 > nclass<-maxcls-mincls+1 > nsub<-nrow(dset) > weight <- matrix(,nage,1) >...

..., sep=",") reading[reading$id %in% c(4, 27, 31, 33, 41, 49, 69, 77, 87), ] xyplot(piat~age | id , data=reading[reading$id %in% c(4, 27, 31, 33, 41, 49, 69, 77, 87), ],panel=function(x,y,*subscripts*){ panel.xyplot(x, y, pch=16) panel.lmline(x,y, lty=4) panel.xyplot(reading$agegrp*[subscripts]*, y, pch=3) panel.lmline(reading$agegrp*[subscripts]*,y) } , ylim=c(0, 80), as.table=T, *subscripts*=T) I just don't know what the subscripts for and the meaning of that. Can someone kindly let me know how it works. Thanks -- View this message in context: http://r.789...

...rug) groups, however my plot shows only one line. What am I doing wrong? My code is reproduced below, my figure is attached to this EMail message. John > #Create simple survival object > GVHDdata<-list(Time=GVHD$Time,Time30=(GVHD$Time)/30, + Age=GVHD$Age,Drug=GVHD$Drug,Died=GVHD$Died, + AgeGrp=cut(GVHD$Age,breaks=c(0,15,25,45))) > > summary(GVHD$Drug) MTX MXT+CSP 32 32 > > > > fit0<-coxph(Surv(Time30,Died)~Drug,data=GVHDdata) > summary(fit0) Call: coxph(formula = Surv(Time30, Died) ~ Drug, data = GVHDdata) n= 64 coef exp(co...

Hello, I am a newbie to the lattice package in R, and I'm trying to make a plot using the xyplot function. I have repeated measures data (2 conditions) for two different groups of subjects (teens and adults). So far, I've made a basic graph using xyplot(y ~x, group=subnum, data=mydata, type="b"). Now I would like to make all the teens' lines one color and the adults'

Hi Folks, I have just started using dplyr and could use some help getting unstuck. It could well be that dplyr is not the package to be using, but let me just pose the question and seek your advice. Here is my basic data frame. head(h) subject ageGrp ear hearingGrp sex freq L2 Ldp Phidp NF SNR 1 HALAF032 A L A F 2 0 -23.54459 55.56005 -43.08282 19.538232 2 HALAF032 A L A F 2 2 -32.64881 86.22040 -23.31558 -9.333224 3 HALAF032 A L A F 2 4 -18.91058 42.121...

Hi, I am trying to runn a age-period-cohort model, but here is what I am having problem with, hope you can help me! This is what I am trying to do: sumzero_a<-((A-min(A))/5+1) - mean((A-min(A))/5+1) where A is my age variable (numeric, the mid-point of a five-year age group), but I got the following error: Error in min(..., na.rm = na.rm) : invalid 'mode' of argument I am pretty

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]