similar to: 'singularity' between fixed effect and random factor in mixed model

Hi, We released individual birds into a room with 2 trees. We counted the number of visits to each of the 2 tree. One of the trees is always a control tree and the other tree is either treatment 1, treatment 2 or treatment3 or treatment 4. Ind Treat ContrTree ExpTree Total visits 1 1 11 16 27 1 2 6 9 15 1 3 5 13 18 1 4 11 25 36 2 1 2 3 5 4 1 6 7 13 4 3 4 4 8 4 4 2 5 7 6 1 1 1 2 6 4 5 16 21

Hi, I was very much hoping someone could help me with the following. I am trying to convert some SAS NLMIXED code to NLME in R (v.2.1), but I get an error message. Does anyone have any suggestions? I think my error is with the random effect "u" which seems to be parametrized differently in the SAS code. In case it's helpful, what I am essentially trying to do is estimate parameters

Hi all I'm facing an issue while running dovecot inside a singularity (https://sylabs.io/singularity/) container dovecot version is 2.3.4.1 (configuration below) running on debian buster, inside a container made with singularity version 3.4.2 unfortunately, when I try to start dovecot, it gives: Singularity test.sif:~> cat /var/log/mail.log Dec 30 17:23:38 testnode dovecot: master:

Hello. I am analysing data from a mixed model perspective using the lme() function. The fixed effects model is a quadratic (Y~X+X2) where X2 is the square of X and the data have a 3-level structure. I fitted a series of three models with the same fixed effects but differing in the random effects (only intercept, intercept + X, intercept +X +X2). The anova shows that all three parameters vary

Dear R users, I have a problem with something called "NaNs" in a nested mixed model. The background is that I have studied the number of insect nymphs emerging from replicated Willow genotypes in the field. I have 15 replicates each of 4 Willow genotypes belonging two 2 Willow species. Now I want to elucidate the effect of Willow genotype on the number of emerging nymphs. Previously I

Dear all, I am trying to use R to fit mixed models. Take the following example, where ind is a random effect and sample is fixed. I wanted to fit Model 1: values = ind + sample Model 2: values =ind * sample Model 3: values=ind(sample) + sample Tried to use the below for mod1, but it did not work. Can anyone give some help on this. Thanks so much. mod1 <- lme(values ~ sample +

Hi, I am hoping some one can help with this. I am using nlme to fit a random coefficients model. It ran for hours before returning Error: Singularity in backsolve at level 0, block 1 The model is > plavix.nlme<-nlme(PLX_NRX~loglike(PLX_NRX,PD4_42D,GAT_34D,VIS_42D,MSL_42D,SPE_ROL,XM2_DUM,THX_DUM,b0,b1,b2,b3,b4,b5,b6,b7,alpha), + data=data, + fixed=list(b0 +

R friends, In a matrix of 1s and 0s, I'm getting a singularity error. Any helpful ideas? lm(formula = activity ~ metaF + metaCl + metaBr + metaI + metaMe + paraF + paraCl + paraBr + paraI + paraMe) Residuals: Min 1Q Median 3Q Max -4.573e-01 -7.884e-02 3.469e-17 6.616e-02 2.427e-01 Coefficients: (1 not defined because of singularities)

Hi, I am using GAMMs to show a relationship of temperature differential over time with a model that looks like this:- gamm(Diff~s(DaysPT)+AirToC,method="REML") where DaysPT is time in days since injury and Diff is repeat measures of temperature differentials with regards to injury sites compared to non-injured sites in individuals over the course of 0-24 days. I use the following

Hi! First I would like to thank you for the help with my first question. However,I am afraid I am still into trouble. I think I have properly used the getGropus commad but seems to be a fatal error when calculating the fixed effects (see bellow) and I cannot understand why. Any inputs? Thanks in advance, /Jose hongos1=read.table("c:/temp/hongos1.txt", header=T) > attach(hongos1)

Hi R-experts, I have a problem using nlme. I use the following code to group my data: Parameterg <- groupedData( result ~ time | Batch, data = Batchdata, labels = list( x = "Time", y = "analysis") ) and then uses the nlme function to fit a nonlinear mixed model that includes Process as a fixed covariate: nlme.model001epr <- nlme(result ~ A0 * exp(- ( exp(A1)

Hallo! I have a problem with fitting data with nls. The first example with y1 (data frame df1) shows an error, the second works fine. Is there a possibility to get a fit (e.g. JMP can fit also data I can not manage to fit with R). Sometimes I also got an error singularity with starting parameters. # x-values x<-c(-1,5,8,11,13,15,16,17,18,19,21,22) # y1-values (first data set)

This command cdmoutcome<- glm(log(value)~factor(year) > +log(gdppcpppconst)+log(gdppcpppconstAII) > +log(co2eemisspc)+log(co2eemisspcAII) > +log(dist) > +fdiboth > +odapartnertohost > +corrupt > +log(infraindex) > +litrate > +africa >

I am trying to write a function to change the case of all of the text in a data frame to lower case. I do not have foreknowledge of the data frame names or the data types of each column. It seems that if one references the data frame by index, then it returns class "data.frame" but if it is referenced by name, it returns class "factor" or whatever the column actually is:

Dear listers, I am trying to assess variance components for a nested, mixed-effects model. I think I got an answer that make sense from R, but I have a warning message and I wanted to check that what I am looking at is actually what I need: my data are organized as transects within stations, stations within habitats, habitats within lagoons. lagoons: random, habitats: fixed the question is:

Hi, The following function does work optimize.W<-function(W,k,G,Groups,cph.call,z){ n<-length(Groups) grp.wt<-rep(0,n) for(i in 1:(length(G))){ ind<-Groups == G[i] if(G[i]!=k){ grp.wt[ind]<-W[i] } elsegrp.wt[ind]<-1 } z<-data.frame(cbind(z,grp.wt=grp.wt)) #needed to make the case weights #part of the data

I am using two mixed effect models. Firstly, what I am trying to do is to compare green roofs abundance with brownfield, green roof with green space abundance, and finally green space with brownfield abundance. I am unsure if I have done the correct model. I have to use a mixed effect model because my data is nested. This is the code and output >

Hi list, I'm working with gamm models of this sort, using Simon Wood's mgcv library: gm<- gamm(Z~te(x,y),data=DATA,random=list(Group=~1)) gm1<-gamm(Z~te(x,y,by=Factor)+Factor,data=DATA,random=list(Group=~1)) with a dataset of about 70000 rows and 110 levels for Group in order to test whether tensor product smooths vary across factor levels. I was wondering if comparing those two

I am having an issue with lmer that I wonder if someone could explain. I am trying to fit a mixed effects model to a set of longitudinal data over a set of individual subjects: (fm1 <- lmer(x ~ time + (time|ID),aa)) I quite often find that the correlation between the random effects is 1.0: Linear mixed model fit by REML Formula: x ~ time + (time | ID) Data: aa AIC BIC logLik deviance

I am using lme4 to fit a mixed effects model to my data. I have a significant interaction between two variables. My question is what is the correct way to get p-values for single terms involved in that interaction. I have been using stepwise backwards deletion and model comparisons to get p-values,and refitting the model using a REML approach to get estimates.However, presumably to get the p