similar to: Lmer BLUPS: was(lmer multilevel)

My question relates to problems that I'm having matching lme and lmer examples in P&B. using Matix 0.995 In the Oxide example in p167-170 I can't get the level 2 coefficient estimates to match the fm1Oxide model in lme is data(Oxide,package="nlme") lme(Thickness~1,Oxide) which I translate in Lmer syntax to fm3Oxide<-lmer(Thickness~

Hi helpeRs, I've been working through some examples in Pinhiero & Bates( 2000) trying to understand how to translate to the new Lme4 syntax but without much luck. Below is what I think I should do, but either the answers don't come out the same or I get errors. In the Oxide problems I'm particularly interested in obtaining the levels coeficients but this options no longer seems

Dear list: I am trying to merge two files together from output I get based on the coef() command. Here is what I am running into. I have two simple linear mixed models > mod1.lme<-lme(math~year, data=sample, random=~year|childid/schoolid) > mod2.lme<-lme(math~year, data=sample, random=~year|childid) I then call the coefficients and store them in the following objects using >

On Feb 3, 2016, at 16:13, Warren Young wrote: > A dying hard disk can do it. HDDs try to silently paper over I/O errors, but what they can?t hide is the time it takes to do this. If your HDD is constantly correcting errors at the oxide layer, it will be reeeeeallly sllllow. > > You can try running SMART tests on it, though that?s not guaranteed to show the problem. Well, it?s not ?a?

This question comes up periodically, probably enough to give it a proper thread and maybe point to this thread for reference (similar to the 'conservative anova' thread not too long ago). Moving from lme syntax, which is the function found in the nlme package, to lmer syntax (found in lme4) is not too difficult. It is probably useful to first explain what the differences are between the

Dear List: I was able to run the following in nlme successfully, but the same model and code (same dataset) failed to run in lme4 and gave me the error message below. Any thoughts? lme(math~year, data=egsingle, random=~year|schoolid/childid) Error in lme(formula = math ~ year, data = egsingle, random = structure(list( : Unable to invert singular factor of downdated X'X

I?m running CentOS 6.7 on my build servers, and on one of the servers the builds are taking almost an order of magnitude longer than usual. There are no runaway processes and there is plenty of free memory. So I suspected that file I/O might be slow, and sure enough, that appears to be the case. I ran a simple dd test and compared the results to a ?normal? build server (412 MB/s vs. 31.7 MB/s).

Hello, I am looking for a way to obtain standard errors for emprirical Bayes estimates of a model fitted with lmer (like the ones plotted on page 14 of the document available at http://www.eric.ed.gov/ERICDocs/data/ericdocs2/content_storage_01/0000000b/80/2b/b3/94.pdf). Harold Doran mentioned (http://tolstoy.newcastle.edu.au/~rking/R/help/05/08/10638.html) that the posterior modes' variances

I've got the latest asterisk 1.8 running on a Netra X1 with Solaris 10 u10. The system itself is happy and phone calls (between two parties) seem fine. Unfortunately, when a caller listens to a Playback recording, there seems to be moments of stutter - perhaps 1 second of stutter for every 10 seconds of Playback. The stutter is not consistent at the same point of the playback file. To

Hi, I''ve been googling quite a bit - but up until now to no use :( My problem is how to configure Apache conf (and perhaps Passenger) for this setup to work: <VirtualHost *:80> ServerName oxen.company_A.com DocumentRoot /var/www/html/oxen/sites/company_A/public PassengerAppRoot /var/www/html/oxen SetEnv OXID 23 </VirtualHost> <VirtualHost *:80>

On Feb 3, 2016, at 1:30 PM, Alfred von Campe <alfred at von-campe.com> wrote: > > I suspected that file I/O might be slow, and sure enough, that appears to be the case?.What could cause this A dying hard disk can do it. HDDs try to silently paper over I/O errors, but what they can?t hide is the time it takes to do this. If your HDD is constantly correcting errors at the oxide

On Feb 3, 2016, at 2:26 PM, Alfred von Campe <alfred at von-campe.com> wrote: > > On Feb 3, 2016, at 16:13, Warren Young wrote: > >> A dying hard disk can do it. HDDs try to silently paper over I/O errors, but what they can?t hide is the time it takes to do this. If your HDD is constantly correcting errors at the oxide layer, it will be reeeeeallly sllllow. >>

Hi, I've recently had a spate of filesystem errors on my laptop which is running a vanilla RedHat 7.2 (kernel 2.4.7) using the ext3 filesystem on all partitions. I have been using the laptop with no apparent problems for some time but last weekend, when I was trying to compile Evolution overnight, I hit a compile problem; it complained that a Makefile didn't exist in a script directory.

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

I am seeking for a method to calculate, given a fitted lme model and some data for a subject, the random effects predictors for this subject. I can only find predictors for the subjects used in creating the fit. Of course I could just add the subject and redo the fit. But I want to avoid just this refitting. Thanks for help wbk

Hi everyone, I need to get a between-component variance (e.g. random effects Anova), but using lmer I don't get the same results (variance component) than using random effects Anova. I am using a database of students, clustered on schools (there is not the same number of students by school). According to the ICC1 command, the interclass correlation is .44 > ICC1(anova1) [1] 0.4414491

Hi What are the disk IOPS required for Asterisk call recording? I am trying to find out number of disks required in RAID array to record 500 calls. Is there any formula to calculate IOPS required by Asterisk call recording? This will help me to find IOPS for different scale. If I assume that Asterisk will write data on disk every second for each call, I will need disk array to support minimum

Hi I am running Arch x64. I have installed Il2-1946which works perfectly. There are some mods I would like to add but cannot install. When I run 'wine setup.exe' for the mod installer I get the error messages below and a dialogue box telling me that I have the error I reference in the title field. All I've seen so far is that I should delete the installshield professional folder

Hello R-Helpers, I was trying to use the lmList function to get the lmList graphic similar to Pinheiro and Bates (pg 33). I did not have a problem creating the graphic when I used the Orthodont data frame or 2 other data sets when there are no missing values. My data has missing values. Do I need to remove the missing values before the lmList function will work? for a small example: > a

Dear R-list, I?m trying to use the lmer of the lme4 package to fit a linear mixed model of the form Y = Xb + Zu + e and I can?t figure out how to control the covariance structure of u. I want u ~ N(0,sigma^2*I). More precisely I?m trying to smooth a curve through data using the "Penalized Splines as BLUPs" method as described in Ruppert, Wand & Carroll (2003). So I have Z = [Z1