search for: m1a

...create in another session or another graph.... Which is why I'm asking for R-help. library(vcdExtra) HEC <- HairEyeColor[,c(1,3,4,2),] # graph 1 mosaic3d(HEC) # hand rotate, scale, then save orientation; make sure rgl window is not hidden #snapshot3d("HEC3D-1a.png", top=TRUE) M1a <- par3d("userMatrix") # try to make this reproducible, across sessions #> dput(M1a) #structure(c(0.890118718147278, 0.105095967650414, -0.44344499707222, #0, -0.454265862703323, 0.282511830329895, -0.84488445520401, #0, 0.0364845432341099, 0.953489363193512, 0.299210667610168, #0,...

Dear R experts. I might be missing something obvious. I have been trying to fix this problem for some weeks. Please help. #data ped <- c(rep(1, 4), rep(2, 3), rep(3, 3)) y <- rnorm(10, 8, 2) # variable set 1 M1a <- sample (c(1, 2,3), 10, replace= T) M1b <- sample (c(1, 2,3), 10, replace= T) M1aP1 <- sample (c(1, 2,3), 10, replace= T) M1bP2 <- sample (c(1, 2,3), 10, replace= T) # variable set 2 M2a <- sample (c(1, 2,3), 10, replace= T) M2b <- sample (c(1, 2,3), 10, replace= T) M2aP1 <-...

Dear listers, One of my former students is trying to fit a model of the negative binomial family with lmer. In the past (two years ago), the following call was working well: m1a<-lmer(mapos~ninter+saison+milieu*zone+(1|code),family=neg.bin(0.451),REML=TRUE,data=manu) But now (R version 2.9.2 and lme4 version 0.999375-32), that gives (even with the library MASS loaded): m1a<-lmer(mapos~ninter+saison+milieu*zone+(1|code),family=neg.bin(0.451),REML=TRUE,data=manu)...

Sorry I had typo in previous email, this typo corrected version: Dear R experts I have simple question, please execuse me: #example data, the real data consists of 20000 pairs of variables K1 <- c(1,2,1, 1, 1,1); K2 <- c(1, 1,2,2, 1,2); K3 <- c(3, 1, 3, 3, 1, 3) M1a <- rep( K1, 100); M1b <- rep(K2, 100) M2a <- rep(K1, 100); M2b <- rep(K1, 100) M3a <- rep(K1, 100); M3b <- rep(K3, 100) mydf <- data.frame(M1a, M1b, M2a, M2b, M3a, M3b) # matrix operation nmat <- matrix(c(paste('M', 1:3, 'a', sep = ''), paste('M&...

...ueFixed" # Choose between restricted and unrestricted estimation cec1=NULL # unrestricted cec1=3 # restrict the coefficient on the dummy to be equal across equations # Estimate the model with `sysGmm` using different weighting matrices: identity, "optimal" and manually specified m1a <- sysGmm(g=ES_g, h=ES_h, wmatrix="ident" , weightsMatrix=NULL, vcov=vc1, crossEquConst=cec1, data=data1); summary(m1a) m1b <- sysGmm(g=ES_g, h=ES_h, wmatrix="optimal", weightsMatrix=NULL, vcov=vc1, crossEquConst=cec1, data=data1); summary(m1b) m1c <- sysGmm(g=ES_g, h=...

...stricted and unrestricted estimation > cec1=NULL # unrestricted > cec1=3 # restrict the coefficient on the dummy to be equal across > equations > > # Estimate the model with `sysGmm` using different weighting matrices: > identity, "optimal" and manually specified > m1a <- sysGmm(g=ES_g, h=ES_h, wmatrix="ident" , weightsMatrix=NULL, > vcov=vc1, crossEquConst=cec1, data=data1); summary(m1a) > m1b <- sysGmm(g=ES_g, h=ES_h, wmatrix="optimal", weightsMatrix=NULL, > vcov=vc1, crossEquConst=cec1, data=data1); summary(m1b) > m1c &lt...

...d, but the solution is to 1) input a series that has already been differenced and not include an order for differencing in the function, or 2) include an xreg argument. I have some sample code: set.seed(66) x = arima.sim(list(order=c(1,1,0), ar=.9), n=100) + 50 y=c(x[-1:-3],x[98:100]) # Method 1 m1a=arima(diff(x),order=c(1,0,0)) m1b=arima(diff(y),order=c(1,0,0)) a=ccf(resid(m1a),resid(m1b)) # Method 2 m2a=arima(x,order=c(1,1,0),xreg=1:length(x)) m2b=arima(y,order=c(1,1,0),xreg=1:length(y)) b=ccf(resid(m2a),resid(m2b)) My question is why do the two methods generate different results for the...

..." <e.robardet@gmail.com> Subject: Re: [R] What happen for Negative binomial link in Lmer fonction? To: r-help@r-project.org Message-ID: <26113408.post@talk.nabble.com> Content-Type: text/plain; charset=us-ascii Thank you for your answers, I have an exemple of that i was using: m1a<-lmer(atpos~ninter+saison+milieu*zone+(1|code),family=neg.bin(0.429),method="Laplace",data=manu) summary(m1a) Generalized linear mixed model fit using Laplace Formula: atpos ~ ninter + saison + milieu * zone + (1 | code) Data: manu Family: Negative Binomial(log link) AIC BI...

..."stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: lme4 lattice Matrix "0.995-2" "0.12-11" "0.995-4" > options(show.signif.stars=FALSE) > m1a<-lmer(cbind(prevented,control.count)~repellant+hour+(1|volunteer)+(1|date), + family=binomial(link='probit'), method='Laplace') Error in dev.resids(y, mu, weights) : argument wt must be a numeric vector of length 1 or length 219 > # probit doesn't converge > m1b<-lme...

...the following additional archive and multimedia filetypes. Several additions reflect the existence of new or more widely used compressed filetypes since the last update >7 years ago, others provide alternate suffixes for existing supported filetypes. Audio ===== .aac (raw AAC) .flac (FLAC) .m1a, .m2a, .mp1, .mp2, .mpa (MPEG audio) .m4a, .m4b, .m4p, .m4r, .m4v (Apple AAC) .opus (Opus) .spx, (Speex) Video ===== .m1v, .m2v, .mpeg, .mpg, .mpv (MPEG video) .qt (Apple Quicktime, additional to .mov) Multimedia containers ===================== .3g2, .3gp (3GPP) .f4v, .flv, (Flash) .m2ts (Bl...

...6; Treatment, 3 Fixed effects: Estimate Std. Error t value (Intercept) 140.500 5.184 27.104 Treatment2 10.500 7.331 1.432 Treatment3 -5.333 7.331 -0.728 Correlation of Fixed Effects: (Intr) Trtmn2 Treatment2 -0.707 Treatment3 -0.707 0.500 > (m1a<-lmer(Glycogen~Treatment+(1|Treatment)+(1|Treatment:Rat)+(1| Treatment:Rat:Liver))) Linear mixed-effects model fit by REML Formula: Glycogen ~ Treatment + (1 | Treatment) + (1 | Treatment:Rat) + (1 | Treatment:Rat:Liver) AIC BIC logLik MLdeviance REMLdeviance 231.6 241.1 -109.8...

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking

...*pJ0PYCW++Ko^n{D#ErwV|#c+-T7{(9q zP~ntzX>r1!M9C2pqgL=gFax9Sg8|Ht1Yjd7xXQ%ft64?)(j=b2x`EqugB!fr8 at xZB z7XNrk;Ev^6ZRSnvrimcYY`t>0{p!rex7i=x{-|I1QBR_c*0*4G`g|ZvWXy&?G==~< zxeMUieb_3-JC!!|OrGu`yI_o~4+`kJ7$P50lz=`6JObC{5Ei3jPEa5`H{325i=`$j zfs+8d^-2 at CRHzq1lxnoLF#Q@>rbfapKRm|?VC?}7UZdm1a- at ziYt}l8q|wMGlPrLk zc7!B_DY3^t_k29IdGyDt$(s(d`w6UVjg}S0C~EGV+L3HN4dvyf;#&@<;+v(P at a2&J z{K{3(6IdEPR9cb+_hy at 0VLDc&qBsO_9xs(ph%ktQ{aECOX^d$!Tqf at gPhCcPS_*(6 zr2t=an^EJdE{P4S_T`~#8FgXD`++beZG5m0fG;DF>#?c^qBX-zS!Qse*npBnGhS0U zA)nQ2*SvdYhNXUmw@^?$%dgT!m*dZ>m}u0>p7)`?p6...