search for: nestl

...tion = position_dodge(width=.75)) but that did not really help, as all points are now dodged and I just want to have a fixed offset for each subgroup of points such that the boxplot and the points are aligned. Any ideas? Kind Regards, Thorn Thaler Mathematician Applied Mathematics Nestec Ltd, Nestl? Research Center PO Box 44 CH-1000 Lausanne 26 Phone: +41 21 785 8220 Fax: +41 21 785 9486

Hello Fellow R Users,I have spent the last week trying to find a work around to this problem and I can't seem to solve it. I simply want to plot my GEE model result with 95% confidence bands. I am using the geepack package to run a basic GEE model involving nestling weights, to a Gaussian distribution, with "exchangeable" error structure. I am examining how nestling weight varies as a function of distance from a plot boundary. The response variable (WEIGHT) and the explanatory variable (DISTANCE) are both continuous. The clustering factor (MOTH...

...ht=600, pointsize=12, quality=100, bg="white") plot(d) dev.off() } # changing the driver, mist.jpg is empty! # So the problem seems to be some interaction with groupedData objects and graphic drivers. # Any hind is appreciated # Regards, # Dominik Dominik Grathwohl Biostatistician Nestl? Research Center PO Box 44, CH-1000 Lausanne 26 Phone: + 41 21 785 8034 Fax: + 41 21 785 8556 e-mail: dominik.grathwohl at rdls.nestle.com

...a well working Internet Explorer 5 on my NT and it sounds as the perfect solution for my problem to rename internet2.dll to internet.dll, to avoid the entire proxy configuration. Doing so, R crashes! What is wrong with internet2.dll? Kind regards, Dominik Dominik Grathwohl Biostatistician Nestl? Research Center PO Box 44, CH-1000 Lausanne 26 Phone: + 41 21 785 8034 Fax: + 41 21 785 8556 e-mail: dominik.grathwohl at rdls.nestle.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.ht...

...p factor(PRO)0 0.113 1.119 0.304 0.370 0.710 factor(PRE)0 0.702 2.018 0.299 2.350 0.019 factor(PRO)0:factor(PRE)0 -0.770 0.463 0.431 -1.789 0.074 ... What would the experts recommend? Kind regards, Dominik Dominik Grathwohl Biostatistician Nestl? Research Center PO Box 44, CH-1000 Lausanne 26 Phone: + 41 21 785 8034 Fax: + 41 21 785 8556 e-mail: dominik.grathwohl at rdls.nestle.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.ht...

...y below that a line for "Overall Mean" (that is all the three items in one legend)? For that I tried to assign the same names to the legends but this did not work either. So any help would be highly appreciated. Kind Regards, Thorn Thaler Mathematician Applied Mathematics Nestec Ltd, Nestl? Research Center PO Box 44 CH-1000 Lausanne 26 Phone: +41 21 785 8220 Fax: +41 21 785 9486

I have a few questions relating to overdispersion in a sex ratio data set that I am working with (note that I already have an analysis with GLMMs for fixed effects, this is just to estimate dispersion). The response variable is binomial because nestlings can only be male or female. I have samples of 1-5 nestlings from each nest (individuals within a nest are not independent, so the response variable is the ratio of sons to daughters) and some females have multiple nests in the data set (so I need to include female identity as a random effect)....

...istake lies? (that means in particular I am not asking for a solution of how to get the thing done - there are plenty of alternatives - but instead to understand why this very approach does not work) Thanks for your help! Kind Regards, Thorn Thaler Mathematician Applied Mathematics Nestec Ltd, Nestl? Research Center PO Box 44 CH-1000 Lausanne 26 Phone: +41 21 785 8220 Fax: +41 21 785 9486

...n additional variable? Something like z <- list(a=1, b=a) or even z <- list(a=b=1)? Both commands don't work of course, but I'm just curious whether this is possible in principle (a rather academic question I've to admit) Thanks + BR, Thorn Thaler Applied Mathematics Group Nestl? Research Center PO Box 44 CH-1000 Lausanne 26, Switzerland Phone: + 41 21 785 8220 Fax: + 41 21 785 8556

...Max -3.11669139 -0.60963281 0.03963523 0.63074983 3.33520406 Number of Observations: 399 Number of Groups: 57 > VarCorr( fm1Weight ) Subj = pdLogChol(1) Variance StdDev (Intercept) 9.604662 3.099139 Residual 1.187553 1.089749 Dominik Grathwohl Biostatistician Nestl? Research Center PO Box 44, CH-1000 Lausanne 26 Phone: + 41 21 785 8034 Fax: + 41 21 785 8556 e-mail: dominik.grathwohl at rdls.nestle.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.htm...

Hi All, What is the R function for computing multinomial distribution, e.g. f(2,1,3; 2/9, 1/6, 11/18, 6)? That is, a total of 6 trials, event 1's p1=2/9, x1=2, event 2's p2=1/6, x2=1, and event 3's p3=11/18, x3=3. thanks, Johnny [[alternative HTML version deleted]]

...3. This solution should be easily extendable to an arbitrary set of columns and should work with missing cases for the treatments like d <- d[-c(2, 21)] Thanks for your input, I am looking forward to your suggestions. Kind Regards, Thorn Thaler Mathematician Applied Mathematics Nestec Ltd, Nestl? Research Center PO Box 44 CH-1000 Lausanne 26 Phone: +41 21 785 8220 Fax: +41 21 785 9486

Hi, If I set the p.adjust="none", does it meant that the output p values from the pairwise.t.test will be the same as those from individual t.tests (set var.equal=T, alternative="t")? I actually got different p values from the two tests. See below. Is it supposed to be this way? Thanks Johnny > x [1] 61.6 52.7 61.3 65.2 62.8 63.7 64.8 58.7 44.9 57.0 64.3 55.1 50.0 41.0

...o the estimate statement in SAS? E.g. proc mixed data=test; class x1 x2; model y=x1|x2 /s; estimate 'x1' x1 -1 1 x1*x2 -0.5 -0.5 0.5 0.5; estimate 'x2' x2 -1 1 x1*x2 -0.5 0.5 -0.5 0.5; run; Thanks for the help, Dominik Dominik Grathwohl Biostatistician Nestl? Research Center PO Box 44, CH-1000 Lausanne 26 Phone: + 41 21 785 8034 Fax: + 41 21 785 8556 e-mail: dominik.grathwohl at rdls.nestle.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.ht...

Hi, I want to free memory used by R. The usage of rm and gc give no result. I'm running an algorithm consuming a huge memory and I need to recover the memory used by R between 2 call of my algorithm. Thank you in advance for your help. e-mail: sofiane.lariani at rdls.nestle.com Sofiane Lariani

...elp. Matt CODE: require(lattice) trellis.par.set("axis.line", list(col="transparent")) wireframe(volcano, main="Correlated Beta Distributions",type='n', xlab=list("Egg Survival",rot=-18), ylab=list("Nestling Survival",rot=50), zlab=list("Density",rot=95),#zaxt="n", screen=list(z=-30,x=-60,y=0),light.source=c(90,-15,90), scales=list(arrows=F,distance=c(.75,.75,0.25), cex=1,col="black"), aspect=c(1,.4...

Hi, I have to analyse the number of provisioning trips to nestlings according to a number of biological and environmental factors. I was thinking of building a mixed-effects model with species and nestid as random effects, using a Poisson distribution, but the data are overdispersed (variance/mean = 5). I then thought of using a mixed-effects model with negativ...

Dear all, I am trying to fit a model with crossed random effects using lme. In this experiment, I have been measuring oxygen consumption (mlmin) in bird nestlings, originating from three different treatments (treat), in a respirometer with 7 different channels (ch). I have also measured body mass (mass) for these birds. id nest treat year mlmin mass ch hack 1EP51711 17 H 2008 1.401719138 10.7 4 2008:17 1EP51709 17 H 2008 1.25716...

Hello, ? I believe that in R, it is not possible to analyze mixed effect-models when the distribucion is not gaussian (p.e. binomial or poisson), isn't? ? Somebody can suggest me alternative? ? thanks ? xavi ? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Summary: proc mixed vs. lme The objective of this summary is to help people to get more familiar with the specification of random effects with proc mixed or lme. Very useful are the examples of Ramon Littell's book: "SAS System for Mixed Models (1996)" (http://ftp.sas.com/samples/A55235) The same data set's are kindly made available by Douglas Bates in the