Search Archives

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091118/c36bb933/attachment-0001.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090428/b6a883b3/attachment-0001.pl>

Thanks a lot for your help! I have tried Peter's code: it works perfectly. Hopefully, this method will also be helpful to others with similar problems. Best regards, Jianghua -- View this message in context: http://n4.nabble.com/How...

Hi, everyone, I met a problem with mixed effect model with lmer. Basically, I want to treat region and mother ID as a nested random effect structure with mother ID nested within a specific region. However, in my dataset, many...

Dear R users, I am frustrated by the residual sum of squares calculation. From Statistical Medhods by Snedecor & Cochran, I know sum(resid(model)^2) should be one method. However, there is a automatic code in R for this operation...

Dear R users, I encounter a problem regarding number of significant digits on y-axis. Below is my basic code: myplotkid<-xyplot(expected_offspringnumber~afr|decade,groups=SES,data1, auto.key=list(space="right"),layout=c(9,1),xlab="",ylab...

Dear R users, I encounter a problem regarding space control in xyplot. Basically, I want to control spaces between label, tick and axis. I remember there is a function called mgp in general plot. Is there a similar function for...

Sorry, I did not notice you were using GAM package. Most R users are using Simon Wood's MGCV package. I recommend you to use it. I have never used GAM package, so I cannot make further comments. Good luck...

Strangely, summary.gam(m1) should give you significance results of parametric terms such as ost, wst, park10, sch50, comm, build and suite. These results should be located above the summary results for smooth terms. Please using summary.gam(m1) to...

Hi, Fabio, I only have an idea on how to calculate deviance explained by the fixed effects. If you remove fixed effects and introduce one null model such as m2<-lmer(response-var ~ 1+(1|Site/Area/Transect),family="binomial...

Sorry for making a mistake! It should be, (m2 at deviance["wrss"]-m1 at deviance["wrss"])/m2 at deviance["wrss"] willow1980 wrote: > > Hi, Fabio, > I only have an idea on how to calculate deviance explained by the fixed > effects. If...

Hi, Dieter, I did add (test="F") in my script, but it does not matter. The following is my whole script. ################################################################# modelsurs_fer13<-gam(sum_surv15~s(FLBS)+SES+s(byear)+s(FLBS,byear),family=quasipoisson) modelsurs_fer14<-gam(sum...

Dear Professor Murdoch, That is exactly the difficulty for me. I don't know how to make a prediction with lmer using "expand.grid"; at the moment, I can use ?model at X%*%fixef(model)? to get predicted values for...

Dear Peter, Many thanks! Your code works perfectly. I hope this method will also be helpful to other users with similar problems. Best regards, Jianghua Peter Ehlers wrote: > > Use the 'scales=' argument together with formatC: > > x <- 1:10 > y <- sample...

Dear R users, I have a problem with multicollinearity in mixed models and I am using lmer in package lme4. From previous mailing list, I learn of a reply "http://www.mail-archive.com/r-help at stat.math.ethz...

Dear R users, I have a problem in plotting 3 dimensional graph using mixed models. My model is sur_prop ~ afr_c+I(afr_c^2)+I(afr_c^3)+byear_c+I(byear_c^2)+I(byear_c^3...

Hi, Simon, I am using mgcv:gam and the version number is mgcv_1.5-2. I also exchanged the order of two models in anova, but this also did not help. >From the differences in DF(0.77246) and...

Dear Simon, Thank you so much! Actually, it seems that Crawley's R book adopted the information from that earlier version and he discussed "by" in the context of that version. I will take a practice according to your suggestion...

I am a little bit confusing about the following help message on how to fit a GAM model with interaction between factor and smooth terms from http://rss.acs.unt.edu/Rdoc/library/mgcv/html/gam.models.html: ?Sometimes models...

Dear R users, I have posted a similar message in the following link: http://old.nabble.com/standard-error-for-the-estimated-value-(lmer-fitted-model)-td26414507.html However, I did not get responses. I guess my question is not...