similar to: Question about lm of R

Dear all, I am using the nlm(base) function for a simple least squares problem. The size (order of magnitude) of the minimized function near the final parameter set is normally about 1e-10. Using the fscale=1e-10 argument in nlm is inefficient (no minimization is performed, iterations=0 or 1). Nlm only works when the minimized function is multiplied by 1e10 (without changing the default

Dear members of the R user community, I am pleased to inform you that the program of the useR! 2009 conference (Rennes, France, July 8 to 10, 2009) is available online: http://www.agrocampus-ouest.fr/math/useR-2009/program.html You can register from the link below: http://www.agrocampus-ouest.fr/math/useR-2009/formulaires/registration.php The conference web page:

Dear members of the R user community, I am pleased to inform you that the program of the useR! 2009 conference (Rennes, France, July 8 to 10, 2009) is available online: http://www.agrocampus-ouest.fr/math/useR-2009/program.html You can register from the link below: http://www.agrocampus-ouest.fr/math/useR-2009/formulaires/registration.php The conference web page:

Dear all, I use a nested design with lm and glm, with factor2 nested within factor1. In order to test for the significance of both factors, I use anova tables on the obtained models such as follows: /> mod1<-lm(A~factor1/factor2) > amod1<-anova(mod1, test="F") Analysis of Variance Table Response: A Df Sum Sq Mean Sq F value Pr(>F) factor1

Hi all, I try to build a model using lmer, with 2 crossed fixed factors (poptype and matingtype) and 2 random factors (pop and family) which I want to nest within poptype. Which of these formulae should I use ? 1. lmer(y~poptype*matingtype + (1|poptype/pop/fam)) 2. lmer(y~poptype*matingtype + (poptype|/pop/fam)) thanks in advance Agnes -- Marie-Agn?s Coutellec UMR INRA-Agrocampus Ouest 985

When I add the following instruction : make.link(log) I have an error that I don''t understand. model<-glm(spa~an, family=gaussian, make.link(log)) error in as.vector(x, "double"): cannot coerce to vector) If someone can help, thanks Martial -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

hi everybody I have to do a lot of Anova with R and I would like to have another type of coefficients coding.. I explain. by default if I have 2 temperatures for an experience. 100°C or 130°C and I want to see the temperature effect on the presure I want to estimate the coefficient of each temperature. I will obtain ,with the anova, juste one coefficients for example +3,56 (for 100°C), and the

It is in the doBy package. Thanks From: Bert Gunter [mailto:bgunter.4567 at gmail.com] Sent: Tuesday, February 13, 2018 4:32 PM To: Pius Mwansa <pmwansa at shaw.ca> Cc: R-help <r-help at r-project.org> Subject: Re: [R] LSmeans and lsmeans Always cc the list unless there is good reason to keep your reply private. There is no LSmeans() function in the lsmeans package.

A cursory reading indicates that they are identical; but others more knowledgeable than I need to confirm or deny this. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Feb 13, 2018 at 3:38 PM, Pius Mwansa <pmwansa at

Always cc the list unless there is good reason to keep your reply private. There is no LSmeans() function in the lsmeans package. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Feb 13, 2018 at 3:20 PM, Pius Mwansa

Dear all, I search the mail list about this topic and learn that no simple way is available to get "lsmeans" in R as in SAS. Dr.John Fox and Dr.Frank E Harrell have given very useful information about "lsmeans" topic. Dr. Frank E Harrell suggests not to think about lsmeans, just to think about what predicted values wanted and to use the predict

Dear all, In connection with a stats course where we want to use R, I am trying to find out if there are analogues to the ESTIMATE and LSMEANS ("population" means) statements in SAS (in proc glm, proc genmod and proc mixed). The students are agronomists, vets and that sort of people, and most of them have some familiarity with SAS. The design package seems to go some of the way, but

Hi, I wonder if anybody could help me in converting this easy SAS program into R. (I'm still trying to do that!) PROC IMPORT OUT= WORK.CHLA_italian DATAFILE= "C:\Documents and Settings\carleal\My Documents\REBECCA\stat\sas\All&nutrients.xls" DBMS=EXCEL2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian;

Hi folks, I have a dataset that consists of counts over a ~30 year period at multiple (>200) sites. Only one count is conducted at each site in each year; however, not all sites are surveyed in all years. I need to impute the missing values because I need an estimate of the total population size (i.e., sum of counts across all sites) in each year as input to another model. >

I need help with calculating lsmeans (adjusted means) of different terms in a linear model including the main effect and the interaction effect terms. I use lm to run the linear models...I previously noted from literature that that "effects" package can be used to generate lsmeans. But I tried to use it but could not figure out which option to use to get means. If anyone can give an

Hi, Is there a mean to automatically choose one version of the 'rbind' function ? By default, R chose the 'rbind' from {base}, i would lke to use the one from {IRanges}... Is it possible to set that using 'par' ? > ?rbind Help on topic 'rbind' was found in the following packages: Package Library IRanges

Is there a difference between LSmeans and lsmeans functions in R? Thanks, Pius

You might get better answers if you 1 - break this down into separate issues 2 - tell us what you want to achieve in words rather than SAS, we all read English but few of us speak SAS 3 - post in plain text not HTML as HTML mangles your post On 29/09/2017 13:47, Andrew Harmon wrote: > Hello all, > > My statistical analysis training up until this point has been entirely done > in

Hi everyone, I am interested in using the SSYEV Lapack subroutine in a R package. Unfortunately I am not familiar with the all Lapack library and I was wondering why the routine SSYEV is not part of the R-LAPACK interface module. Is there any plan to extend the R-Lapack module in the next R release? Should I try to directly call the SSYEV routine of the external Lapack library from my R

R users, What is the equivalent of SAS " '/diff' and '/slice'" in R ? I am successful in fitting the same model "PROC MIXED" fits using lme() : ###R code: model <- lme( y ~ g + h + yr + yr:h, random = ~ 1 | b, data = yield) *** SAS code: proc mixed data=yield covtest cl method=reml; class b g h yr; model y = g h yr yr*h; random b; lsmeans g / diff;