Re: [R] Syntax differences between aov and lmer for 2-way repeated measures design using a mixed model

2019-04-12 Thread Bert Gunter
You should talk with your professor. This list is about R programming. Essentially statistical issues, which this appears mostly to be, are generally off topic. Questions about mixed effects models -- RM and longitudinal designs are typically analysed as such -- and especially using the nlme and/

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread William Dunlap via R-help
You should ask the maintainer of the package about this: bug.report(package="pglm") or maintainer("pglm") should give you contact information. The help file for pglm seems all wrong - it says pglm's output has class "pglm" with components like "fitted.values", but the the example calls to pglm ret

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread Simon Berrebi
Thank you Bert, I wasn’t aware of ?str. The only mention of fitted-values is: $ maximum: atomic [1:1] -9824 ..- attr(*, "fitted.values")= num [1:3460] 1.39 1.3 1.3 1.32 1.27 ... When I try attr(mymodel$maximum, ”fitted.values”), I get the same results as attr(AIC(mymodel), ”fitted.val

Re: [R] Creating a mean line plot

2019-04-12 Thread rain1290--- via R-help
Hi Eric, Ah, I apologize, and thank you for your response!  I just figured out a way to average my x-values, so at least that is solved. I will still include the data for the two variables (1-dimensional) of interest that I was trying to average, just to show what was done: get2.teratons #(90 va

[R] Syntax differences between aov and lmer for 2-way repeated measures design using a mixed model

2019-04-12 Thread Uri Eduardo Ramírez Pasos
Hi everyone, I'm working with the following data frame using R. It consists of measurements obtained from 7 subjects with two independent variables (IV1 and IV2) with two levels each (OFF/ON, ALT/ISO, respectively): >myData Subject DV IV1 IV2 1 2.567839 OFF AL

Re: [R] Help for coding

2019-04-12 Thread Rui Barradas
Hello, Remove ind from the calculations. And maybe substitute length(Variable) for 1312. And standard deviations is the square root of the variance, you are missing the square root. Hope this helps, Rui Barradas Às 17:50 de 12/04/2019, bienvenidoz...@gmail.com escreveu: Dear all, I am tryi

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread Bert Gunter
?fitted Read and follow. As I said, you appear to need to spend time with a tutorial. I do not provide this service, though others may. -- Bert On Fri, Apr 12, 2019 at 10:07 AM Simon Berrebi wrote: > Thank you Bert, > > I wasn’t aware of ?str. The only mention of fitted-values is: > > $ ma

[R] Help for coding

2019-04-12 Thread bienvenidoz...@gmail.com
Dear all, I am trying to write the script to calculate the mean and the standard deviation without using the functions mean and standard deviations since morning. Someone can help These scripts below are what i did but when i run and check with the function mean and standard deviation it is not

Re: [R] Creating a mean line plot

2019-04-12 Thread Eric Berger
I don't have your data. Are the x-values the same in both plots? Does this example cover the situation? f1 <- function(x) { x^3 - 2 } f2 <- function(x) { 2 - x^2 } xV <- seq(from=0,to=2,length=50) y1 <- f1(xV) y2 <- f2(xV) y3 <- .5*(y1+y2) plot(x=xV,y=y1,col="blue",lwd=2,type='l',xlab="x",ylab="y

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread Bert Gunter
?fitted ?predict ## This is what one usually does, but I have not checked pglm. You also need to get friendly with ?str ... and probably also spend time with an R tutorial or two to become familiar with R modeling conventions. Bert Gunter "The trouble with having an open mind is that people kee

[R] Are fitted.values available in pglm?

2019-04-12 Thread Simon Berrebi
Hello everyone, I am using the pglm function in R to fit a Poisson fixed-effects model. According to the documentation , the pglm object should have fitted.values. However, fitted.values(mymodel) returns "NULL". When I run AIC(mymodel) the

[R] Creating a mean line plot

2019-04-12 Thread rain1290--- via R-help
Hi there, I am trying to create a mean line plot that shows the mean of a series of separate line plots that correspond to two climate models. Let's first try getting the mean of two line plots. To create the separate line plots, here is what I did to set up the x and y axis variables: Gett

[R] [R-pkgs] Releasing new package on CRAN

2019-04-12 Thread Roman Minguez Salido
Dear R-users: First, we apologize for possible crossed messages. It is a pleasure for us to have the opportunity to present a new R package, called 'spsur' available in https://cran.r-project.org/web/packages/spsur/index.html This is a specific R package for the estimation of Spatial Seemingl