similar to: lme() doesn't converge on IGF example

I'm trying to get code along the following lines to work: temp.name <- paste(TimePt,'df',sep='.') # invent a relevant name/symbol as a character string. assign(temp.name,IGF.df[IGF.df$TPt==TimePt,]) # this works. The relevant variable is now a data frame lm(b ~ Strain+BWt+PWt+PanPix, data=temp.name)) # this gives an error, namely Error in eval(predvars, data, env) : invalid

Dear list Maybe someone can help with the following problem: I have a source file containing the following code for plotting: pdf(file = 'data/mario/boxplot.pdf') bwplot(sig100 ~ target | group, data = fish, main = 'Boxplot der Signale in 100 Zellen nach Gruppe', xlab = '', ylab = 'Anzahl der Signale') dev.off() pdf(file = 'data/mario/xyplot.pdf')

Hi: I am working with someone remotely to allow them access to our data. The follow command using "download.file" works perfectly on my Mac: > > download.file(url="http://oceanwatch.pfeg.noaa.gov:8081/thredds/ > wcs/satellite/AG/ssta/14day? > request=GetCoverage&version=1.0.0&service=WCS&format=NetCDF3&coverage= >

Dear R gurus, There is an interesting problem with accessing specific items in a column of data frame that has incorrectly been given a duplicate name, even though addressing the item by row and column number. Although the column is correctly listed, an item addressed by row and column number gives the item with the correct row and the original not the duplicated column. Here are the

Dear Contributors: I have two vectors x and z, and I want to display the histograms of both vectors in the same graph, x in red bars, z in blue bars. If you have any clue on how to do that, I will be very glad to hear it!!!!!! Thanks in advance again, Juan Pablo Fededa [[alternative HTML version deleted]]

Hola! I have some problems with gam in mgcv. Firts a detail: it would be nice igf gam would accept an na.action argument, but that not the main point. I want to have a smooth term for time over a year, the same pattern repeating in succesive years. It would be natural then to impose the condition s(0)=s(12). Is this possible within mgcv? I tried to obtain this with trigonometric terms, aca:

The .Renviron and .First functions do not seem to work the same way on a Mac OS 10.4 as on a Windows XP machine. From working in Windows I am used to creating a new directory for each data analysis project. In the new directory I place First, an .Renviron file consisting of the following text: R_HISTFILE="history.txt" R_HISTSIZE=1000000 Second, an .RData file containing a .First

I'm transferring 1.1 Mb files over very poor GSM EDGE connection. My rsync command is: rsync --partial --remove-source-files --timeout=120 --append --progress --rsh=ssh -z LOCAL_FILE root at SERVER:REMOTE_PATH File on remote server "grows" in size in steps of 262144 bytes. That is a lot, because system needs to transfer at least 262144 (before compression) every time connection

Hi, I just upgraded to Ruby 1.8.6 on my Mac development machine and here''s what I''m seeing when I try to run mongrel (I have the latest version of mongrel + fastthread installed). It seems to be an OS/X thing, because this same combo works fine on our Linux server. Here''s the error, have you guys seen anything like this?: ** Starting Mongrel listening at

Hi all, I've had one of my most miserable R weeks in memory. I'm trying to deal with huge datasets (>1GB each) but am running up against those pesky memory limits. The libraries filehash and g.data are not very suitable for what I need. I haven't gotten into the sql thing yet. Most recently I've been trying to install the new package ff (not yet on the CRAN repository). I

Dear list, I am having a segfault while working with large lists in R 2.4.1 patched (64-bit compilation) on an apple xserve xeon with 16GB RAM (I don't know if this is reproducible on other configs since this is the only computer with enough memory I have access to...) You will find below a simple somewhat reproducible (usually it crashes around 5.5GB of vsize, sometimes it gets into

Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of

Version 0.9975-11 of the lme4 package has been uploaded to CRAN. The source package should be available on the mirrors in a day or two and binary packages should follow soon after. There are several changes in this release of the package. The most important is the availability of a development version of lmer called, for the time being, lmer2. At present lmer2 only fits linear mixed models.

Version 0.9975-11 of the lme4 package has been uploaded to CRAN. The source package should be available on the mirrors in a day or two and binary packages should follow soon after. There are several changes in this release of the package. The most important is the availability of a development version of lmer called, for the time being, lmer2. At present lmer2 only fits linear mixed models.

I am having trouble printing a table out to the GUI display when the table is created and printed within a loop. I get a "Error: syntax error message" If I comment out the print statement, the loop runs fine and I can print out the last iteration of the table. ...[multiple loops and calculations ending with.....]... +

> (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Error in as.double(start) : Calloc could not allocate (888475968 of 4) memory ************************* > sessionInfo() R version 2.4.1 (2006-12-18) powerpc-apple-darwin8.8.0 locale: C attached base packages: [1] "grid" "datasets" "stats" "graphics" "grDevices"

Hello, Just a quick message to say that, ------------------------------------------------------------------------ ------- when choosing (with or without the tcltk window) "Bern 1 (swiss)" as CRAN mirror to install a package (function install.packages) for R-2.5.0, the list of available packages for osx seems not complete: AMORE, AdaptFit, AlgDesign, AnalyzeFMRI, SemiPar,

I seem to come to this problem alot, and I can find my way out of it with a loop, but I wish, and wonder if there is a better way. Here's an example (lmer1-5 are a series of lmer objects): bs=data.frame(bic=BIC(lmer1,lmer2,lmer3,lmer4,lmer5)$BIC) rownames(bs)=c('lmer1','lmer2','lmer3','lmer4','lmer5') best=rownames(bs)[bs==min(bs)] > best [1]

Hello EVERYONE, I need an URGENT help from you please! How can I see the "estfun" (empirical estimating function) and "df" (degree of freedom) from the following mixed-model please? (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Many thanks in advance for your kind help. Sattar

Dear members, I would like to switch from nlme to lme4 and try to translate some of my models that worked fine with lme. I have problems with the pdMat classes. Below a toy dataset with a fixed effect F and a random effect R. I gave also 2 similar lme models. The one containing pdLogChol (lme1) is easy to translate (as it is an explicit notation of the default model) The more parsimonious