similar to: lme, spline (revised question)

Is there a way to specify a Z matrix using the lmer function, where the model is written as y = X*Beta + Z*u + e? I am trying to reproduce smoothing methods illustrated in the paper "Smoothing with Mixed Model Software" my Long Ngo and M.P. Wand. published in the /Journal of Statistical Software/ in 2004 using the lme4 and Matrix packages. The code and data sets used can be found at

Hi there, I am reading the 2004 paper "Smoothing with mixed model software" in Journal of Statistical Software, by Ngo and Wand. I tried to run their first example in Section 2.1 using R but I had some problems. Here is the code: library(nlme) fossil <- read.table("fossil.dat",header=T) x <- fossil$age y <- 100000*fossil$strontium.ratio knots <-

Dear All, I have a problem running this program on R. Z is a matrix of spline which is random > fit<-lme(anc~X,random=pdIdent(~Z)) Error in getGroups.data.frame(dataMix, groups) : Invalid formula for groups What I have done wrong?

On Jun 6, 2018, at 4:06 AM, Alice Wonder <alice at domblogger.net> wrote: > > I'm not anti-microsoft but I'm worried they will make changes that I don't like (e.g. requiring ms account, changing billing, etc.) so I figured better take control now. What had you using GitHub in the first place, then? Is it just the web UI, since you?ve apparently decided that none of

I was wondering if one can have the coefficients of VAR model sorted according to variable names rather than lags. If you notice below, the output is sorted according to lags. >VAR(cbind(fossil,labour),p=2,type="const") VAR Estimation Results: ======================= Estimated coefficients for equation fossil: =========================================== Call: fossil = fossil.l1

Hi all - R2.0.1, OS X Perhaps while there is some discussion of lme going on..... I am trying to execute a glmm using glmmPQL from the MASS libray, using the example data set from McCullagh and Nelder's (1989, p442) table 14.4 (it happens to be the glmm example for GENSTAT as well). The data are binary, representing mating success (1,0) for crosses between males and females from two

Hello, my name is Toby, and i live in Sydney Australia. i would like to buy a 25Kg bag of strontium nitrate and have it delivered to an address in sydney. i understand that long distance postage can sometimes be very costly, so can you please give a a quote for the 25Kg bag, plus the postage cost to my address in sydney. thankyou for your time

Hi all, I'm running samba-3.0.25b w/Linux 2.6.16.27 (SUSE 10.1), newly compiled. I've joined the AD domain (which is a subdomain) whose PDC is Win2K3 server and which is running in native mode but allowing NT4 members. From the HOWTO, I followed the instructions - manually created the computer account, then did the "net rpc join ...." and successfully joined the AD

I'm trying to do some smoothing with lme and am having some difficulty bringing commands over from S-PLUS to R. I have the following setup (modified from Ngo and Wand, 2004): set.seed(1) x <- runif(200) y <- sin(3*pi*x) + rnorm(200)*.4 ## library(splines) z <- ns(x, 4) The following runs without error on S-PLUS f <- lme(y ~ 1, random = pdIdent(~ -1 + z)) But in R I get

Dear all, I am trying to fit a model with crossed random effects using lme. In this experiment, I have been measuring oxygen consumption (mlmin) in bird nestlings, originating from three different treatments (treat), in a respirometer with 7 different channels (ch). I have also measured body mass (mass) for these birds. id nest treat year mlmin mass ch hack 1EP51711 17

Dear R users, Does anyone knows how to run a glmm with one fixed factor and 2 random numeric variables (indices)? Is there any way to force in the model a separate interaction of those random variables with the fixed one? I hope you can help me. #eg. Reserve <- rep(c("In","Out"), 100) fReserve <- factor(Reserve) DivBoulders <- rep

Windows XP R 2.10 I am trying to run lme and get the following error: > fitRandom <- lme(values ~ subject, + data=withindata) Error in getGroups.data.frame(dataMix, groups) : Invalid formula for groups my data follows, below which is a copy of all my code > > print(withindata) subject values 1 1 2.3199639 2 1 -8.5795802 3 1 -4.1901241 4

Hello all, In lme4 if you want to model two non-nested random effects you code it like this: mod1 <- lmer(y~x + (1|randomvar1) + (1|randomvar2)) How would you go about to model something similar in nlme? In my database I have two variables for which I have repeated measures, lets call them "individual" and "year". But none of the "individuals" were measured in

I've solved the problem of the detached Eaton_SDK branch. What must have happened here is the branch creator did a non-Subversion cp -r followed by an add of the target directory. This registered a tree that duplicated trunk, *without* ancestry information connecting it to its source revision. The correct fix is to check each new text section against a list of MD5 hashes of old text

Greetings, I might have missed something in the list here.. But could someone point out some introductory material / docs regarding cmirror? Is there something (even beta would do) in CentOS to read and understand about this? OT- On Fossils: While I cannot claim "Fossil state" anywhere near some of the people who were doing cards -- I knew about them though -- when I was born, I have

Dear Rusers, I am trying to feed my database with data from a file. But since my file (2010101000.txt) there`s no headers I am facing problem because the result data.frame is not separated my columns. How could I set variables names for each columns in *dbWriteTable*? I have tried the command below but I don`t know how to fill the *field.type*option. dbWriteTable(con, "b20101010",

I'm getting an INVITE from the Cisco softswitch looking like this (from sip debug trace): Sip read: INVITE sip:9043940358@206.165.120.52;user=phone SIP/2.0 Via: SIP/2.0/UDP sia-ATLCA146.telefyne.com:5060;branch=z9hG4bK_1146_38l5 From: <sip:7278673170@sia-ATLCA146.telefyne.com;user=phone>;tag=1_1146_f153592_732 j To: <sip:9043940358@206.165.120.52;user=phone> Call-ID:

Hi I'm running a glm analysis and I got the following message: glm.i01<-glm(U~y+v+i+b+v:y+v:i+v:b+y:i+y:b+i:b,family=Gamma(link="log")) Error: cannot allocate vector of size 153594 Kb I thougth this was not suppose to happen in version 1.5 ?! How can I solve this problem ?? Can this a problem of swap space ?? My machine is a pentium III 833 with 512MB of RAM and I'm using

Full_Name: Simon Frost Version: 1.0.0 OS: Windows Submission from: (NULL) (129.215.60.79) Hi everyone, I've just installed R for Windows - I've previously been using it on Linux - and I've run into a problem when using lme. The lme routine throws up an error in the getGroups.data.frame command, complaining about an error in the formula: Error in getGroups.data.frame(dataMix,

Hola! I am encountering the following problem, in a multilevel analysis, using glmmPQL from MASS. This occurs with bothj rw1062 and r-devel, respectively with nlme versions 3.1-38 and 3.1-39 (windows XP). > S817.mod1 <- glmmPQL( S817 ~ MIEMBROScat+S901+S902A+S923+URBRUR+REGION+ + S102+S103+S106A+S108+S110A+S109A+S202+S401+S557A+S557B+ + YHOGFcat,