search for: agroecolog

...e "x" by "w" in the formula object "y". How can this be done? Somehow, the methods that can be used in character objects do not work 1:1 in formula objects... Many thanks and best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

..."pairwise.complete.obs") How do I extract only those correlations that are >0.6? w[w>0.6] #obviously doesn?t work, and I can?t find a way around it. I would very much appreciate any help! Best wishes Christoph (using R 2.5.1 on Windows XP) -- Dr. Christoph Scherber DNPW, Agroecology University of Goettingen Waldweg 26 D-37073 Goettingen Germany phone +49(0)551 39 8807 fax +49(0)551 39 8806 homepage www.gwdg.de/~cscherb1

...m1 m2 Do you know a way around this? I guess I would need to use expression() or substitute() but I don?t know how to construct m1, m2... without the quotation marks. Thanks very much for your help! Best wishes, Christoph (using R 2.5.1 on Windows XP SP2) -- Dr. Christoph Scherber DNPW, Agroecology University of Goettingen Waldweg 26 D-37073 Goettingen Germany phone +49(0)551 39 8807 fax +49(0)551 39 8806 homepage www.gwdg.de/~cscherb1

...ysis? What would be your recommendation to the students, given the fact that most of the standard glm models obviously don?t seem to produce good fits here? Many thanks and best wishes Christoph (using R 2.8.0 on Windows XP) -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

...I have read that (2) is the only approach that works. But how could I then still include the nesting information from (1)? Many thanks for your help! Best wishes Christoph (using R 2.9.0 and the nlme library on Windows XP) -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

...for any help or comments. Kind regards Raphael ------------------------------------------------------------------------------------ Raphael Felber, PhD Scientific Officer, Climate & Air Pollution Federal Department of Economic Affairs, Education and Research EAER Agroscope Research Division, Agroecology and Environment Reckenholzstrasse 191, CH-8046 Z?rich Phone +41 58 468 75 11 Fax +41 58 468 72 01 raphael.felber at agroscope.admin.ch<mailto:raphael.felber at agroscope.admin.ch> www.agroscope.ch<http://www.agroscope.ch/> [[alternative HTML version deleted]]

...what I am stuck with; I thought using diag() on the SSPH and SSPE matrices should give me the sums of squares for hypothesis and error; so division should yield the F values? I would be most grateful for any help! Best wishes Christoph -- Dr. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen, Germany www.gwdg.de/~cscherb1

Dear all, I got this error message > library(dprep) > mardia(Savg) Error in cov(data) : 'x' is empty But with the same data, I got > library(mvnormtest) > mshapiro.test(Savg) Shapiro-Wilk normality test data: Z W = 0.9411, p-value = 0.6739 What does the error message "Error in cov(data) : 'x' is empty" mean? Thanks a lot! Jiao

...ael > ------------------------------------------------------------ > ------------------------ > Raphael Felber, PhD > Scientific Officer, Climate & Air Pollution > > Federal Department of Economic Affairs, > Education and Research EAER > Agroscope > Research Division, Agroecology and Environment > > Reckenholzstrasse 191, CH-8046 Z?rich > Phone +41 58 468 75 11 > Fax +41 58 468 72 01 > raphael.felber at agroscope.admin.ch<mailto:raphael.felber at agroscope.admin.ch > > > www.agroscope.ch<http://www.agroscope.ch/> > > >...

...t;, "sa30", "sa31", "sa32", "sa33", "sa34", "sa35", "sa36", "sa37", "sa38", "sa39", "sa40")) Eisenring Michael, Dr. Federal Department of Economic Affairs, Education and Research EAER Agroecology and Environment Biosafety Reckenholzstrasse 191, CH-8046 Z?rich Tel. +41 58 468 7181 Fax +41 58 468 7201 michael.eisenring at agroscope.admin.ch<mailto:michael.eisenring at agroscope.admin.ch> www.agroscope.ch<http://www.agroscope.ch/> [[alternative HTML version deleted]]

...------------------------------------- >> ------------------------ >> Raphael Felber, PhD >> Scientific Officer, Climate & Air Pollution >> >> Federal Department of Economic Affairs, >> Education and Research EAER >> Agroscope >> Research Division, Agroecology and Environment >> >> Reckenholzstrasse 191, CH-8046 Z?rich >> Phone +41 58 468 75 11 <+41%2058%20468%2075%2011> >> Fax +41 58 468 72 01 <+41%2058%20468%2072%2001> >> raphael.felber at agroscope.admin.ch<mailto:raphael.felber@ >> agroscope.adm...

...Would it be sensible, then, to choose the best-fit model as the one that has *both* minimal AIC *and* minimal number of d.f.? Are there any references for such an approach? Many thanks for your help! Best wishes, Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

...ared to the other models. Would it be allowed to select model3 instead, simply because it uses up less df and the delta-AIC between model2 and model3 is just 1? Many thanks for any suggestions/comments! Best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

Dear all, Is there a way to retrieve standard errors from nls models? The help page tells me that arguments such as se.fit are ignored... Many thanks and best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

Hello, We are attempting to use nlme to fit a linear mixed model to explain bird abundance as a function of habitat: lme(abundance~habitat-1,data=data,method="ML",random=~1|sampleunit) The data consist of repeated counts of birds in sample units across multiple years, and we have two questions: 1) Is it necessary (and, if so, how) to specify the repeated measure (years)? As written,

...can I extract the total and residual d.f. from a gnls object? I have tried str(summary(gnls.model)) and str(gnls.model) as well as gnls(), but couldn?t find the entry in the resulting lists. Many thanks! Best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

Dear all, Is it possible to create a pdf file with layers using the pdf() device in R? Many thanks for your help! Christoph (using R 2.15.1 on Windows 7 64-Bit) -- PD Dr Christoph Scherber Georg-August University Goettingen Department of Crop Science Agroecology Grisebachstrasse 6 D-37077 Goettingen Germany phone 0049 (0)551 39 8807 fax 0049 (0)551 39 8806 http://www.gwdg.de/~cscherb1

I have data that looks like YC Age Num 82 11 2 83 10 0 84 9 8 85 8 21 86 7 49 87 6 18 88 5 79 89 4 28 90 3 273 91 2 175 with a program mod1=lm(log(Num+1)~YC, data=box44) plot(log(Num+1)~YC, data=box44, pch=19, xlab="Year Class", ylab="Loge Number at age", ylim=c(0,6), xlim=c(91,82)) abline(lm(log(Num+1)~YC), col="blue", lwd=2) summary(mod1) I need to

...= TRUE, legend =TRUE) plotArrow(dta.plsda2, ind.names = dta.factor, legend =TRUE) plotVar(dta.plsda2, cex = 2) plot(dta.plsda,typeVc = "x-score",parAsColFcVn = dta.factor,parEllipsesL = TRUE) Eisenring Michael, Dr. Federal Department of Economic Affairs, Education and Research EAER Agroecology and Environment Biosafety Reckenholzstrasse 191, CH-8046 Z?rich Tel. +41 58 468 7181 Fax +41 58 468 7201 michael.eisenring at agroscope.admin.ch<mailto:michael.eisenring at agroscope.admin.ch> www.agroscope.ch<http://www.agroscope.ch/> [[alternative HTML version deleted]]

...> >Raphael >------------------------------------------------------------------------------------ >Raphael Felber, PhD >Scientific Officer, Climate & Air Pollution > >Federal Department of Economic Affairs, >Education and Research EAER >Agroscope >Research Division, Agroecology and Environment > >Reckenholzstrasse 191, CH-8046 Z?rich >Phone +41 58 468 75 11 >Fax +41 58 468 72 01 >raphael.felber at agroscope.admin.ch<mailto:raphael.felber at agroscope.admin.ch> >www.agroscope.ch<http://www.agroscope.ch/> > > > [[alternative HTML...