[R] mixed models lmer function help!!

2017-06-19 Thread Shelby Leonard via R-help 
Hi,I have tumor growth curve data for a bunch of different mice in various 
groups. I want to compare the growth curves of the different groups to see if 
timing of drug delivery changed tumor growth.I am trying to run a mixed models 
with repeated measures over time with each mouse as a random effect with linear 
and quadratic terms for time.This took me a long time to figure out and I just 
wanted to make sure I did it correctly and I am interpreting it correctly. This 
is the code I ran
Rtumor<-lmer(volume~Group+Time+(1|Subject), 
data=Rtumor)summary(Rtumor)Rtumor.null=lmer(volume~Time+(1|Subject), 
data=Rtumor, REML=FALSE)Rtumor.full=lmer(volume~Group+Time+(1|Subject), 
data=Rtumor, REML=FALSE)anova(Rtumor.null,Rtumor.full)
Here is my output Rtumor<-lmer(volume~Group+Time+(1|Subject), data=Rtumor)> 
summary(Rtumor)Linear mixed model fit by REML ['lmerMod']Formula: volume ~ 
Group + Time + (1 | Subject)   Data: Rtumor
REML criterion at convergence: 1541.2
Scaled residuals:     Min      1Q  Median      3Q     Max -1.8006 -0.6348 
-0.0658  0.3903  4.7551 
Random effects: Groups   Name        Variance  Std.Dev.  Subject  (Intercept) 
3.197e-09 5.654e-05 Residual             3.348e+05 5.786e+02Number of obs: 101, 
groups:  Subject, 11
Fixed effects:            Estimate Std. Error t value(Intercept) -495.520    
303.619  -1.632Group         24.350    115.615   0.211Time          79.653      
7.886  10.101
Correlation of Fixed Effects:      (Intr) Group Group -0.933       Time  -0.300 
-0.007
> Rtumor.null=lmer(volume~Time+(1|Subject), data=Rtumor, REML=FALSE)> 
> Rtumor.full=lmer(volume~Group+Time+(1|Subject), data=Rtumor, REML=FALSE)> 
> anova(Rtumor.null,Rtumor.full)Data: RtumorModels:Rtumor.null: volume ~ Time + 
> (1 | Subject)Rtumor.full: volume ~ Group + Time + (1 | Subject)            Df 
>    AIC    BIC  logLik deviance  Chisq Chi Df Pr(>Chisq)Rtumor.null  4 1576.5 
> 1586.9 -784.24   1568.5                         Rtumor.full  5 1578.4 1591.5 
> -784.22   1568.4 0.0457      1     0.8307There were 50 or more warnings (use 
> warnings() to see the first 50)


My questions are1) Did I do this correctly?2) Do I still need to run it again 
with quadratic terms for time?, If so, how do I do this?3) If I am 
understanding these results correctly, they say There is no difference between 
these groups on volume growth curves
[[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.

Re: [R] fftImg() error: fftw_access_func

2016-03-26 Thread Shelby Leonard via R-help 
So do i need to resend the email to someone else? sorry i am just confused 

On Saturday, March 26, 2016 9:21 AM, John Kane  wrote:
 

 It would be helpful if you actually supplied your code and a minimal data set 
to for people to examine.

Please have a look at 
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
 and/or http://adv-r.had.co.nz/Reproducibility.html.

I'm not clear if you are reporting a general R problem or a specific package 
with the ripa package.  If it looks like it is the latter you probably bring it 
to the attention of the package maintainer who may or may not monitor this 
mailing group.


John Kane
Kingston ON Canada


> -Original Message-
> From: ehand...@macalester.edu
> Sent: Fri, 25 Mar 2016 15:09:42 -0500
> To: r-help@r-project.org
> Subject: [R] fftImg() error: fftw_access_func
> 
> Hello-
> 
> My name is Eric Handler and I am an academic technologist supporting
> the Science Division(7 academic departments) at Macalester College in
> Saint Paul, MN. The faculty use R for a variety of teaching and
> research tasks around campus. I administer our RStudio instance and
> have encountered an error I can't resolve. A student working on an
> independent research project has reported that he received the
> following error when attempting to use the ripa function fftImg():
> 
> Error in .C("fftw_access_func", as.complex(img), as.integer(w),
> as.integer(h),  :
>  "fftw_access_func" not available for .C() for package "ripa"
> 
> I've been able to recreate this error in RStudio as well as directly
> in R. I've also recreated the error across different platforms(Ubuntu,
> Mac OS X 10.10 and 10.11). My test platform's sessionInfo() output is
> below:
> 
> R version 3.2.3 (2015-12-10)
> Platform: x86_64-apple-darwin13.4.0 (64-bit)
> Running under: OS X 10.11.3 (El Capitan)
> 
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
> 
> attached base packages:
> [1] parallel  tcltk    stats    graphics  grDevices utils    datasets
> [8] methods  base
> 
> other attached packages:
> [1] fftw_1.0-3 ripa_2.0-2
> 
> loaded via a namespace (and not attached):
> [1] tools_3.2.3 Rcpp_0.12.3
> 
> A google search for "fftw_access_func" doesn't reveal anything modern
> on this topic, only mentions of OS 10.4 and 10.5 and rimage, which
> doesn't seem to exist anymore(perhaps it was a predecessor to RIPA?)
> 
> Can someone help me get the student functional with fftImg() or
> alternately, tell me it is a known issue and an alternative option(if
> available) for the student?
> 
> Thanks,
> Eric
> 
> --
> Eric Handler
> Academic Information Associate - Science Division
> Macalester College - Saint Paul, MN
> Olin-Rice 124
> Office: 651-696-6016
> View my calendar: http://goo.gl/SbxLOu
> 
> __
> 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.


[[elided Yahoo spam]]

__
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.


  
[[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.