Re: [R] AIC values with lmer and anova function

4th Question: Why have you not posted this on the R-sig-mixed-models list, where it clearly belongs, rather than here? -- Bert On Sat, Dec 29, 2012 at 12:47 PM, Diego Pujoni wrote: > Dear colleagues, > > I have a data from a repeated measures design that I'm analysing through a > mixed model. N

[R] Error in plot.envfit(ef, p.max = 0.1) : (subscript) logical subscript too long

Hello there, I'm trying to plot vectors with p<0.1 in a NMDS ordination plot using p.max. Below the scripts I'm using. I guess I'm missing something! could you please give me a hand? species<-metaMDS(species_matrix)ef<-envfit(species,environmentaldata_file,permu=999,na.rm=TRUE)efplot(species, d

Re: [R] Working with date

Ron, I hope this will be helpful. require(lubridate, quietly = TRUE) require(stringr, quietly = TRUE) test_date_str1 <- "6/3/2012" test_date_1 <- mdy(test_date_str1, quiet = TRUE) test_date_1 test_date_str2 <- "2-3-2011" my_date_strings <- c(test_date_str1, test_date_str2) my_dates <- mdy(my_date

Re: [R] efficiently multiply different matrices in 3-d array with different vectors?

HI, Tx for the fast method. dim(Rt) #b #2  dim(X) #[1] 2 5 Is it possible to obtain the same result as X without converting X to R?   According to OP ". which will give me a (two-dimensional) array or matrix of dimension n x p with the i'th row storing the above. " A.K. --

[R] standard errors in johansen test

Hi, I have the very same problem of this here http://grokbase.com/t/r/r-help/111d386yn9/r-standard-errors-in-johansen-test, I replicated all the steps in the Papar: http://cran.r-project.org/web/packages/vars/vignettes/vars.pdf untill here: R> vecm <- ca.jo(Canada[, c("rw", "prod", "e", "U")],

[R] AIC values with lmer and anova function

Dear colleagues, I have a data from a repeated measures design that I'm analysing through a mixed model. Nine independent sampling units (flasks with culture medium with algae) were randomly divided into 3 groups ("c", "t1", "t2"). There is no need for inclusion of the random effect of the interce

Re: [R] efficiently multiply different matrices in 3-d array with different vectors?

What I had in mind was a tensor multiplication. I think, using tensor arithmetic for multidimensional arrays looks more compacts and efficient (~ 2-3 times). I am guessing that performance will be much more pronounced for n-D arrays (n>3). # Component Wise set.seed(14) Z<-array(sample(1:100,8

[R] I need intercept in plm model

Hi,R 2.15.2 plm() function on Windows 7 when i perform a plm regression, i can't manage to obtain the intercept, but I need it.it gives me just the beta coefficient. my formula: fixed <- plm(deltaS ~ L1.deltaS + L2.deltaS, data=Mody_R, index=c("country_id", "date"), model="within") my output: Co

Re: [R] parallel error message extraction (in mclapply)?

Hi, Ivo, At Fri, 28 Dec 2012 16:34:03 -0800, ivo welch wrote: > so, it isn't obvious to me how I get to the try-error object > list, or how to find out what triggered the abort. A try object is an object that captures the execution context at the time of an error. Consider the following code:

[R] shell function in R

Hi Everyone: I am running NONMEM from R using following command. shell("nmfe72 TOT_EST_DOSE2.CTL DOSE2.OUT",wait=T) This command runs NONMEM successfully but does not show the progress. Usually, intermediate outputs are displayed on CMD console when this same command is run fr

Re: [R] thanks -- Re: syntax for identifying more than one

HI Greg, Sorry, I misunderstand your question. I am not sure whether it works with numSummary() from library(Rcmdr). You could use other ways, such as: test<-read.table(text="  id year incidents  100    1    0  101    1    1  102    1    21  100    1    27  101    1    3  102