search for: mer

...value is possible but rare): Val range -7 to 7 and other Car range 1 to 7. This a code to fake database of research Subj<-rep(1:20,each=60) Sti<-rep(c("G","S","P"),40) SP<-rep(c("S","P"),60) AG<-rep(c("A","P"),60) Mer<-rep(c("M","NM"),60) Car<-round(runif(120,1,7),0) Val<-round(runif(120,-7,7),0) base<-data.frame(Subj,Sti,SP,AG,Mer,Car,Val) In my hypothesis: logit(Sti=="G") ~ SP+AG+Mer+Car+Val+SP*Car+Mer*Val+AG*Val + errorG logit(Sti=="S") ~ SP+AG+Mer+Car+Va...

Hi all I have problems in accessing a mer object called model.01 from a workspace that was created with R 2.8.1 and saved with save into an .RData file (on Windows XP or Ubuntu 8.10, don't remember anymore). Now I want to open it in R 2.9.0 on Ubuntu 8.10. I use # load workspace load("name.RData") which seems to work: ls()...

Hi all, I've now got a problem with some modified lmer code (function lmer1 pasted at end) - I've made only three changes to the lmer code (marked), and I'm not really looking for comments on this function, but would like to know why execution of the following commands that use it almost invariably (but not quite predictably) leads to the R ses...

I use lmer to fit a mixed effect model.It give some warnings.what does this warnings mean? and what should I do? > (fm2.mlm <- lmer(qd ~ edu + jiankang + peixun +hunyin + cadcj + age + age2 + sex + dangyuan + Comp.1 + Comp.2+trust.cz1 +(trust.cz1|commid), data = individual,na.action = "na.exclude&...

Hello, Could somebody please tell me what am I doing wrong in following? I try extract coefficients (using arm-package) from the lmer frunction, but I get the following warning: a<-data.frame(coef(res)) Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class "coef.mer" into a data.fram I think I have done it before the same way and it has worked, but...

I am having a problem extracting from "mer" objects.    I have constructed my problem using existing datasets.   Using the following commands:   require(lme4) fm1 <- lmer(Yield ~ 1 + (1 | Batch), Dyestuff) fixef(fm1) I get the following error message: "Error in UseMethod("fixef") : no applicable method for "fix...

Hi everyone, I'm sure this has been asked before, but I could not find the right search terms to find it. If I want to see what summarizing a model fit with lm() does, I can write "summary.lm". But what if I want to see what summarizing a model fit with lmer() (lme4 package) does? showMethods(summary) Function: summary (package base) object="ANY" object="character" (inherited from: object="ANY") object="mer" object="sparseMatrix" object="summary.mer" tells me there is a "mer" m...

...y=pred), random = ~ (1|vpnr),data=sdata,family=binomial) However this would obviously test each of the smooths aginst a horizontal line which is not what I want. Any help is greatly appreciated, thanks! Andrea Output of 1st model gamm.11: ---------------------------------------- logLik(gamm.11$mer);deviance(gamm.11$mer);attributes(summary(gamm.11$mer))$AICtab[1];gamm.11$mer@deviance["disc"] 'log Lik.' -49054.65 (df=5) ML 98109.3 AIC 98119.3 disc 97600 > summary(gamm.11$mer);summary(gamm.11$gam) Generalized linear mixed model fit by the Laplace approximat...

Folks In trying to use lmer for a hierarchical model, I encountered the following message: Error in UseMethod("VarCorr") : no applicable method for 'VarCorr' applied to an object of class "mer" foo.mer <- lmer(y ~ TP + (TP|M),data=joe.q) > head(joe.q[,1:5]) TP M AB Trt y...

...eftab.default <- function(object,...) { print(class(summary(object))) coef(summary(object)) } The NAMESPACE tries to pull in the necessary bits and pieces from lme4 to extract summaries and coefficients: export("coeftab","coeftab.default") importClassesFrom(lme4,"mer","summary.mer") importMethodsFrom(lme4,"coef","summary","show","print") exportMethods("coef","summary","show","print") exportClasses("mer","summary.mer") S3method(coeftab,default) T...

...911 19 2 0.6 199728100212 20 0 0.4 199820100811 21 1 0.8 199826200611 22 2 0.6 199827203811 23 2 1.0 200038109911 24 0 0.6 200126202511 25 0 0.4 200226100311 26 1 0.6 200226100411 27 1 0.4 200226100611 28 1 0.4 200226126011 29 1 0.4 200226203712 30 2 0.6 200227220313 With the following model; lmer(a~b + (1|c), family=poisson, data=tab), What I want to do is to see if number of recruits (a) is dependent on the brood sex ratio (b) including brood identity (c) as random factor. However, I get the following error message; lmer(a~b + (1|c), family=poisson, data=tab) Error in devAGQ(PQLpars, 1...

...the R-forge maintainers and was redirected to the general R mailing lists] As a side note: The reason for importing lme4's summary is that otherwise my function summary.fixef, which extracts fixed effects from nlme and lme4 models for printing e.g. in reports, doesn't use/find summary.mer on "mer" objects -- even if lme4 is loaded earlier. On the other hand, lme models (from package nlme) are handled correctly despite the fact that I do not import the corresponding summary function. Thank you very much Benjamin

Hi, Has anyone seen any documentation on the Hash Object defined in the prototype library? To be precise the following methods - 1. keys 2. values 3. toQueryString 4. inspect 5. Merge I am trying to play around with them and have come up with the following examples - var h = $H({name: "john doe", email: "john-Ch9RrZxMC0c@public.gmane.org", msg: "say hello to me"}); var h_in = h.inspect(); alert(h_in); //returns #<Hash:{''name'...

...a column of one data.frame to another, my problem is, that i only want to match the data, if 2 columns are equal. i have a example: data.frame1 cu.nr. name value A 1 Evo 100 B 1 Mer 80 C 2 Ford 50 And now i want to match the value to data.frame 2 in dependence of cu.nr. and name. data.frame2 cu.nr. name A 33 Evo B 1...

Dear list, I have to read some clinical data a file coming from Filemaker on Macintosh (Ugh ! But it could be worse and come from Excel...). Exporting via Excel is out of question since the file has 467 columns and 121 lines (+ headers), which is out of reach of Excel. So I received an "mer"" files, which is what Filemaker exports as a text file. It seems to be a semicolon delimited file, with all data item enclosed in double quotes, and header items not delimited. So, first attempt is : In bash : # Convert from apparent Latin-1 to UTF-8 charpent at yod:/media/EMTEC/MARS-...

What is teh reason some functions in the AICcmodavg package do not work with 'mer' class models? One such example would be the 'importance' function. Thanks Ronny -- View this message in context: http://r.789695.n4.nabble.com/AICcmodavg-functions-and-mer-class-models-tp3714534p3714534.html Sent from the R help mailing list archive at Nabble.com.

What is the easiest way to extract a value from a 'mer' object from glmer? The first I need is the trtpair Std.Dev. which in this case is 0.17542? I've managed to get the fixed effects numbers from summary(fednmaout)@coeffs but no luck with the Random effects stuff. TIA Jim > fednmaout Linear mixed model fit by REML Formula: lor...

Dear all, I try to run the example of lmer and get the following error message. > library(lme4) > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) [[1]] Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not found This error message is similar to what I get with othe...

I want to call summary on a mer object (from lme4) within my package but I can't seem to get the namespace to import the necessary method. I've simplified my package to this one function: --------------------------- ss <- function(m) { summary(m) } --------------------------- And my namespace file looks like this,...

Dear all I want to compute Monte Carlo p-values in lmer-models based on sampled data sets. To speed up calculations, I've tried to use internal functions from the Matrix package (as suggested ealier on the list by Doug Bates). So I did: fm2 <- lmer(resistance ~ ET + position + (1|Grp), Semiconductor,method='ML') simdata<-simulate(...