[R] specifying lme function with a priori hypothesis concerning between-group variation in slopes

2010-10-18 Thread Bill Shipley
how to specify the model in lme but I don't know how to set up the inferential test that b_j=Kj for all j groups versus the alternative hypothesis that b_j is not equal to Kj for at least one group. Any help in explaining how to do this using the mle function in R is appreciated. Thanks. Bill

[R] compiling, linking and including multiple FORTRAN subroutines in R

2010-07-13 Thread Bill Shipley
of the two. Thanks for any help or leads. Bill Shipley bill.ship...@usherbrooke.ca __ 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

[R] avoiding termination of nls given convergence failure

2009-03-12 Thread Bill Shipley
when this happens, so that the loop can continue? Bill Shipley North American Editor, Annals of Botany Département de biologie Université de Sherbrooke Sherbrooke (Québec) J1K 2R1 Canada (819) 821-8000, poste 62079 (819) 821-8049 FAX http://pages.usherbrooke.ca/jshipley/recherche

[R] superposing barplots having different scales

2008-05-28 Thread Bill Shipley
different scales to define the numerical y values for each subgroup. This graph would have one scale on the left-hand y-axis and a second scale on the right-hand y-axis. I cannot simply superimpose two bar plots because I have to make sure that the subgroup bars are beside each other. Bill Shipley

[R] solution to differences in sequential and marginal ANOVA using a mixed model

2008-05-07 Thread Bill Shipley
Yesterday I posted the following question to the help list. Thanks to John Fox (copied below) who pointed out the solution. Original question: I have come across a result that I cannot explain, and am hopingthat someone else can provide an answer. A student fitted a mixed model usingthe

[R] Type I or III SS with mixed model function lme

2008-05-06 Thread Bill Shipley
) give the same output since they correspond to Type I and III sums of squares in the SAS terminology. At least, this is the case with normal (i.e. not mixed) linear models. However, he finds very different results of these two types of ANOVA tables. Why? Bill Shipley North American Editor