similar to: R crash with modified lmer code

I'm trying to specify a model using lmer with a binary response and interaction term, but I get an error I can't parse (see below). Here is some sample data: Subject Concord Age Disc SVC999MX148SU-F yes u int TOU999JU030S1 yes u int TOU999JU030S1 yes u int TOU999JU030S1 yes u int TUT578MX037S2 yes g int COL140MX114S2 yes yf

Dear all, I ran the following model without any problem previously to the update of lme4: fm2<-lmer(data=NGud,family=poisson, seed~hab*seedtray +(1|site)+(1|site:hab)) I have 25 sites, 2 habitats ("hab") per site, 8 seedtrays per habitat (4 as control, and 4 as treatment), and I'm interested in comparing the number of seed in seedtrays as a function of the treatment and its

theliver wrote: > > I'm currently doing simple, one way ANOVA using R. > > I'm analyzing a count based upon an Hour value (0-23) in which it occurs. > The count varies greatly over the entire day (usually with a range of 400 > or so), but certain hour groupings could be useful for later data > processing. > > So I've loaded my data: > data.nas <-

Dear R users, I am trying to fit a GLMM to the following dataset; tab a b c 1 1 0.6 199320100313 2 1 0.8 199427100412 3 1 0.8 199427202112 4 1 0.2 199428100611 5 1 1.0 199428101011 6 1 0.8 199428101111 7 0 0.8 199527103011 8 1 0.6 199527200711 9 0 0.8 199527202411 10 0 0.6 199529100412 11 1 0.2 199626201111 12 2 0.8 199627200612 13 1 0.4 199628100111 14 1 0.8

Dear R users I have built several glmm models using glmmPQL in the following structure: m1<-glmmPQL(dev~env*har*treat+dens, random = ~1|pop/rep, family = Gamma) (full script below, data attached) I have tried all the methods I can find to obtain some sort of model fit score or to compare between models using following the deletion of terms (i.e. AIC, logLik, anova.lme(m1,m2)), but I

Hi, I have three questions concerning GLMMs. First, I ' m looking for a measure for the significance of the random variable in a glmm. I'm fitting a glmm (lmer) to telemetry-locations of 12 wildcat-individuals against random locations (binomial response). The individual is the random variable. Now I want to know, if the individual ("TIER") has a significant effect on the model

Hello to all, I recently downloaded R to my PC and am enjoying getting acquainted with it. Thank you to everyone involved in the R-project! I am interested in doing a log-linear analysis with R on a data set with dichotomous variables. There are 11 variables (columns) and around 1000 subjects (rows). How do I aggregate my data, i.e. how do I make a new dataset that includes the variable giving

Sparc Solaris 8 SunStudio 11 compilers R-2.3.1 spatstat 1.10-1 R-2.3.1 and spatstat were installed using SunStudio11 compilers on a Sparc Solaris 8 machine. No errors were reported in the compilations, however, when I try to load spatstat I get an error: > library(spatstat) Loading required package: mgcv This is mgcv 1.3-20 Error in dyn.load(x, as.logical(local), as.logical(now)) :

Hi R helpers, I'm trying to fit a rather complex model to some simulated data using lme and am not getting the correct results. It seems there might be some identifiability issues that could possibly be dealt with by specifying starting parameters - but I can't see how to do this. I'm comparing results from R to those got when using GenStat... The raw data are available on the

There seems to be a minor bug in barplot.default when used with log scale w= here one or more values is NA: dat <- matrix(1:25, 5) dat[2,3] <- NA barplot(dat, beside =3D T) #Plots and appropriate barplot with gaps for m= issing data barplot(dat, beside =3D T, log =3D "y") #Error in if (min(height + offset) <=3D 0) stop("log scale error: at least = one 'height +

Affected: Various 2.x.x releases of samba on Solaris 2.6/8 and Netware 5.1 SP5. Problem: Unable to copy files from a mapped Netware drive to a mapped samba share on a PC desktop (W2K). Error: Several, but commonly "Cannot copy MYFILE. Access is denied. The source file may be in use" or "An unexpected network error has occurred". Has anyone seen similar? It appears that some

Hi, I'm stuck with a fairly basic re-ordering problem. I want to extract part of a matrix as a vector, and reorder it to match with an erratic sequence of x,y coordinates. Here's an example that shows what I want and how close I've got (but for some reason my mind just can't get the final step at the moment!): #Define a grid: xy100 <- expand.grid(1:10, 1:10) # Define a matrix

Hi, I'm trying to make sense of the image3d plots in misc3d. Unfortunately I can't find any authors on the help pages of the function to e-mail directly. In the following toy example, I'm puzzled as to why (1) the positions of the image plot on the axes include negative values when I think I've specified all to be 1:3 and (2) the logic that needs to be adhered to in order to

Hi everyone, I've been having difficulty writing wrapper functions for some functions where those same functions include other functions with eval() calls where the environment is specified. A very simple example using function lmer from lme4: lmerWrapper <- function(formula, data, family = gaussian, method = c("REML", "ML", "PQL", "Laplace",

Dear R helpers   I have following csv file which is an input   id       par_value    coupon_rate     frequency_coupon   tenure    ytm   1            1000             10                      1                     5          12   # Here frequency_coupon is coded s.t. 0 means Daily compounding, 1 means monthly compouding, 2 means Quarterly, 3 means Half yearly and 4 means only once. Thus in the

Hi problem: plot(xxx) results in the graphics device window adding the axes values for what ever xxx is but failing to actually add lines or produce the plot until the consol window is activated and moved over the graphics window. The graphics window then shows what has been plotted. This started happening 4 days ago, up until that time plot(xxx) would produce the plot on the graphics device

# I have written this code in Notepad++ and copied here. ## ONS - PPA    Duration = function(par_value, coupon_rate, freq_coupon, tenure, ytm) { macaulay_duration  =   NULL modified_duration    =   NULL freq_coupon_new    =   NULL if(freq_coupon <= 0) {     freq_coupon_new = 365 } if(freq_coupon > 0 & freq_coupon <= 1) {     freq_coupon_new = 12 } if(freq_coupon > 1 &

Dear R helpers   (I have already written the required R code which is giving me correct results for a given single set of data. I just wish to wish to use it for multiple data.)   I have defined a function (as given below) which helps me calculate Macaulay Duration and Modified Duration and its working fine with given set of data.   My Code -   ## ONS - PPA   duration = function(par_value,

Hi, I'm needing some help finding a function to read a large text file into an array in R. The data are essentially presence / absence / na data for many species and come as a grid with each species name (after two spaces) at the beginning of the matrix defining the map for that species. An excerpt could therefore be: SPECIES1 999001099 900110109 011101000 901100101 110100019 901110019

Adam, With the blog entry[1] you''ve made about gzip for ZFS, it raises a couple of questions... 1) It would appear that a ZFS filesystem can support files of varying compression algorithm. If a file is compressed using method A but method B is now active, if I truncate the file and rewrite it, is A or B used? 2) The question of whether or not to use bzip2 was raised in the