[R] urgent question about Lmer models
Dear r-help, I have a question about the std error of a lmer model with library (LmerTest). I have 12 sites with 6 treatments over each site. I measured some response variable (biomass_of_insects) (with a gaussian distribution). I did : library (LmerTest) model<- lmer((biomass_of_insects) ~ as.factor(treatments) + (1 | sites)) However, the response of the model show always the same std error (see below): Fixed effects: Estimate Std. Error df t value Pr(>|t|) (Intercept) 501.333 80.656 66.000 6.216 3.91e-08 *** as.factor(treat)F 126.667114.065 66.000 1.1100.271 as.factor(treat)I-8.333114.065 66.000 -0.0730.942 as.factor(treat)I+F1/3 -75.000114.065 66.000 -0.6580.513 as.factor(treat)I+F2/3 18.333114.065 66.000 0.1610.873 as.factor(treat)I+F3/3 15.917114.065 66.000 0.1400.889 Do you know why std.error is always the same ? Thank you very much, Aline [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] da.norm function
Hello all I'm running da.norm function in R for climate data rngseed(1234567) theta1=da.norm(mydata, thetahat, steps=1000,showits=T) param1=getparam.norm(mydata,theta1) As I understand the 1000 steps represent the markov chain values. Is there a way to plot them? Something like plot(1:1000, param1$mu[]). I just can't find a way to extract them out of my theta1. Thank you, Andrey. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] da.norm function
Hello all I'm running da.norm function in R for climate data rngseed(1234567) theta1=da.norm(mydata, thetahat, steps=1000,showits=T) param1=getparam.norm(mydata,theta1) As I understand the 1000 steps represent the markov chain values. Is there a way to plot them? Something like plot(1:1000, param1$mu[]). I just can't find a way to extract them out of my theta1. Thank you, Andrey. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] da.norm function
Hello all I'm running da.norm function in R for climate data rngseed(1234567) theta1=da.norm(mydata, thetahat, steps=1000,showits=T) param1=getparam.norm(mydata,theta1) As I understand the 1000 steps represent the markov chain values. Is there a way to plot them? Something like plot(1:1000, param1$mu[]). I just can't find a way to extract them out of my theta1. Thank you, Andrey. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
Re: [R] How to mean for groups of matrix rows?
Dear All Big Thanks! R is wonderful language! -- *Thanks, Andrei* [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] How to mean for groups of matrix rows?
Dear All, For a vector, I use this xu<-1:20 t<-rep((1:4),each=5) tapply(xu,t,mean) 1 2 3 4 3 8 13 18 and for a matrix the only way I may guess is: > xu [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,]14321432 [2,]21432143 [3,]32143214 [4,]14321432 [5,]21432143 [6,]32143214 > t [1] 1 1 1 2 2 2 > y [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,]00000000 [2,]00000000 > for (i in 1:dim(xu)[2]) y[,i]<-tapply(xu[,i],t,mean) > y [,1][,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,]2 2.33 2.6732 2.33 2.673 [2,]2 2.33 2.6732 2.33 2.673 I do not like the need to create a matrix (y) for the result. Is there a better way? Thanks, Andrei. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] Re-installing R
Dear R users How does one completely uninstall R from their machine? Going to control panel>programs does not do it for me. After installing the new version it will still remember my previous workspace and all packages I've installed. Thank you. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] Linear regression interaction terms
Hello I performed a linear regression, my equation is Y = âo+ â1A + â2B + â3AB. Is there a way to separate interaction terms, say â3AB and plot it against a certain variable? Thanks, Andrew [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] Time series plot shift
I have a plot and my x axis is Time in days, I need to shift it by +15 days, how do i do it? Thanks. [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
Re: [R] How to import simple java/mathematica expression to R
2010/7/24 David Winsemius : > > On Jul 23, 2010, at 10:29 PM, jim holtman wrote: > > It has some interesting properties with that identity definition of EE(x). > Local maxima at 1 and 0, "blows up" beyond -1 and 2, and several local > minima nearby: > > Math.Fn <- function(alp) { } > plot( seq(-.3,1.5,by=0.01), Math.Fn(seq(-.3,1.5,by=0.01) ), cex=0.2) > Indeed, EE(x) is not a function with real x, but a symbol with subscript - i.e. EE(1) and EE(2) are two variables e1 and e2. Thanks for Your interest. -- Regards, -Andrey __ 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 commented, minimal, self-contained, reproducible code.
Re: [R] How to import simple java/mathematica expression to R
Hello, 2010/7/24 jim holtman : > Well, I took you equation and put the following at the start: > > Power <- function(x,y) x^y > EE <- function(x) x > alp <- 2 > > x <- 900*Power(-0.2030178326474623 + 0.23024073983368956*(1 - alp) + >... Actually, there are many functions like this (about 18*2), for each situation with input data. So there is an idea to put them in different files (via Export from Mathematica) and then use them somehow (without modifications) in R program. Thank You for the answer. -- Regards, -Andrey __ 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 commented, minimal, self-contained, reproducible code.
Re: [R] How to import simple java/mathematica expression to R
Hello, 2010/7/23 jim holtman : > It would be nice if you could post what the data looks like that you > want to import. R can import any text file and then you have string > manipulation that you can do to parse it. So the basic answer is > probably yes, but we do need to understand the format of the data to > give a more precise answer. I put the target expression to minimize (with some constrains) here: http://analytic-products4you.com/target.txt Is it possible to import it as a function to minimize? > > What is the problem that you are trying to solve? > We solve a problem for parameters estimation with ties. Thank You for the answer. -- Regards, -Andrey __ 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 commented, minimal, self-contained, reproducible code.
[R] How to import simple java/mathematica expression to R
Hello, Is it possible to import some (large enough) polynomial expression from java (or mathematica) to R as a function? I found very interesting package rJava but I did not see this feature there. We used Mathematica as calculation engine for our program with statistical calculations and now we want to sell our package to customers. As we cannot force people to buy the Mathematica only for our program, so we are looking for another calculation engine. R is very good candidate. -- Regards, -Andrey __ 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 commented, minimal, self-contained, reproducible code.
[R] How do I get removed from this mailing list?
How do I get removed from this mailing list? [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.