similar to: how to store lme/lmer fit result

Hello, I have a set of data that has a Date column looks like this: 12/9/2007 12/16/2007 1/1/2008 1/3/2008 1/12/2008 etc. I'd like the date to look something like the follow (so that I could sort by date easily). 20071209 20071216 20080101 20080103 20080112 How to do it? Thank you very much Julia -- View this message in context:

Hello, I have a data set that looks like this. ID value 111 5 111 6 111 2 178 7 178 3 138 3 138 8 138 7 138 6 . . . I'd like to calculate the mean and var for each object identified by the ID. I can in theory just loop through the whole thing..., but is there a easier way/command which let me calculate the mean/var by ID? Thanks, Julia -- View this

Hello, I have a data file (.csv) that has a size of about 2.6 GB. I am not able to read in the whole data set because of the memory limit. I actually only need some columns (3 columns) of the data set, is there a way to read in specified columns? I am using windows. Thanks, Julia -- View this message in context:

Hello, I have a vector: x <- c("A", "A", "A", "B", "A", "A", "C") I'd like to compare each of elements of vector x from its previous element (except for the 1st element which does not have previous element). So I'd like to get a vector y (of same length) that looks something like (0, 0, 0, 1, 1, 0, 1) or (F, F, F,

Dear R users: I am running lmer(), and I had the following warning message: "In mer_finalize(ans): iteration limit reached without convergence (9) What is the default limit for lmer()? And how do I increase the limit? Many thanks, Julia -- View this message in context: http://www.nabble.com/How-to-increase-iteration-limit-for-lmer%28%29--tp20252974p20252974.html Sent from the R help

Dear R users: I have a data set that looks something like this: ID time y 1 2/01/2008 4 1 2/09/2008 12 1 9/01/2008 8 2 1/06/2008 3 2 3/01/2008 4 2 3/09/2008 9 2 6/03/2008 4 3 1/02/2008 3 3 1/10/2008 8 3 2/02/2008 7 3 2/10/2008 3 I'd like to sum up the weekly data to the monthly level, so that it looks something

Hello, I have a time series plot drawn using 3 different colored lines, each line corresponds to different category group. I'd like to put legends on the plot. I am using "legend" to do this, however, I can either specify lty or col in legend. Will I be able to do the following in the legend box: a solid black line then g1, a solid red line then g2, a solid blue line then g3?

Dear R users, I need to compare two scatter plots, plot(x1, y1) plot(x2, y2) and would like to plot them in the same figure. How do I do it? Thank you. -- View this message in context: http://www.nabble.com/two-scatter-plots-in-one-tp21963375p21963375.html Sent from the R help mailing list archive at Nabble.com.

Dear R-devel members, I am looking for a fast Cholesky update/downdate. The matrix A being symmetric positive definite (n, n) and factorized as A = L %*% t(L), the goal is to factor the new matrix A +- C %*% t(C) where C is (n, r). For instance, C is 1-column when adding/removing an observation in a linear regression. Of special interest is the case where A is sparse. Looking at the

Hello, and sorry for asking a question without the data - hope it can still be answered: I've run two things on the same data: # Using lme: mix.lme <- lme(DV ~a+b+c+d+e+f+h+i, random = random = ~ e+f+h+i| group, data = mydata) # Using lmer mix.lmer <- lmer(DV ~a+b+c+d+(1|group)+(e|group)+(f|group)+(h|group)+(i|group), data = mydata) lme provided an output (fixed effects and random

Dear R-users, I did do a thorough search and read many articles and forum threads on the lme and lmer methods and their pitfalls and problems. I, being not a good statistician but a mere "user", came to the conclusion that the most correct form of reporting statistics for a mixed linear model would be to report the parameter estimates and SEs, and, if the sample size is considerably

Dear R users, I am running lmer() and having memory problem: "reached total allocation of 1535Mb". The lmer() function was successful (no errors or warnings). However when I do "summary(fit.lmer), I got this "Error: cannot allocate vector of size 19.5 Mb). I used memory.size() and got 1880. I used memory.limit(2500) to increase the memory size to 2500Mb, however I still get

Dear R users, Suppose I have a vector that consists of characters like ABC, A02, RCA, etc., and there are about 700 of possible characters. For example, x <- c("ABC", "ABC", "ABC", "A02", "ABC", "RCA", "ABC", "ABC") I'd like to get a frequency matrix that looks something like this: ABC 6 A02 1 RCA 1 I

Dear R users, I would like to test the randomness in a series of N values (N>=2). I know that runs.test works for dichotomous factor only: x <- rep(c(1,2), 50) runs.test(factor(x)) However it doesn't work for series that can take any N values (N>2): x <- rep(c(1,2,5,4),50) runs.test(factor(x)) Error in runs.test(factor(x)) : x does not contain dichotomous data Are there any R

Dear R users, I have this vector that consists numeric numbers. Is there a command that detects the repeated numbers in a vector and returns the index of the repeated numbers (or the actual numbers)? For example, v <- c(3,4,5,7,4). The command would return me index 2 and 5 (or the repeated number, 4). Thank you very much, Julia -- View this message in context:

Dear R users: I have a vector of dates as follows: t <- c("2007-01-05", "2007-05-14", "2007-12-28", "2008-01-09", "2008-04-24", "2009-02-14") I'd like to calculate number of days between those dates (time interval). How to do it? Thank you, Julia -- View this message in context:

Dear All, After playing with lmer for couple of days, I have to say that I am amazed! I've been using quite some multilevel/mixed modeling packages, lme4 is a strong candidate for the overall winner, especially for multilevel generzlized linear models. Now go back to my two-level poisson model with cross-classified model. I've been testing various different model specificatios for the

[mailed jointly to maintainers@octave.org, r-devel@stat.math.ethz.ch to reach users with relevant experience. watch where your replies go. what would make using LAPACK and ScaLAPACK easier for Octave and R developers? -- ejr, not subscribed] We plan to update the LAPACK and ScaLAPACK libraries and would like to have feedback from users on what functionalities they think are missing and would

Hey, I am having trouble with lmer. I am looking at the presence/absence of water shrews against habitat and other factors e.g so I used this: m1<-lmer(Presencebsence~Habitatype*Width+(1|Sitename))summary(m1) But i keep getting this error up Error in mer_finalize(ans) : Downdated X'X is not positive definite, 16.> summary(m1)Error in asMethod(object) : matrix is not symmetric [1,2] What

Hi - I am doing some monte carlo simulations comparing bayesian (using Plummer's jags) and maximum likelihood (using lmer from package lme4 by Bates et al). I would like to know if there is a way I can flag nonconvergence and exceptions. Currently the simulations just stop and the output reads things like: Error in optim(.Call("lmer_coef", x, 2, PACKAGE = "Matrix"), fn,