similar to: lmer, p-values and all that

Dear community, I'm having this silly problem. I've a linear model. After fixing it, I wanted to know which data had studentized residuals larger than 3, so i tried this: d1 <- cooks.distance(lmmodel) r <- sqrt(abs(rstandard(lmmodel))) rstu <- abs(rstudent(lmmodel)) a <- cbind( mydata, d1, r,rstu) alargerthan3 <- a[rstu >3, ] And suddenly a[rstu >3, ] has

Hi, A following function was kindly provided by GGally?s maintainer, Barret Schloerke. function(data, mapping, ...) { p <- ggplot(data = data, mapping = mapping) + geom_point(color = I("blue")) + geom_smooth(method = "lm", color = I("black"), ...) + theme_blank() + theme(panel.border=element_rect(fill=NA, linetype =

Hello there. I uses the following codes for the purpose of variable selection. > lmModel <- lm(y~.,data.frame(y=y, x=x)) > step <- stepAIC(lmModel, direction="both") > step$anova Stepwise Model Path Analysis of Deviance Table Initial Model: y ~ x.Market.Price + x.Quantity + x.Country + x.Incoterm + x.Channel + x.PaymentTerm Final

Hi, Thanks for your response, it nearly worked! But it only wrote one coloumn of data and not the three columns I need. Using fixef(m1) doesnt give the same results as coef(m1) when you are using more than one random effect. I need the coefficients for each individual so I use coef(m1) to get this which results in an object of class lmer.coef, 3 columns by 700 rows. as.data.frame() wont work on

Hi, all. In setting up my package for post-processing regression models, I am trying to use standard notation as much as possible: thus, I use coef() to access estimated coefficients. I wrote a function called se.coef() to grab standard errors, and se.fixef() and se.ranef() to grab se's from coefficients estimated from lmer(). I also need a function to access sigma-hat (the residual sd

I just got this error: uninitialized constant PriceCurf After much head scratching I remembered the Rails plural to singular mapping rules. My table is price_curves, my object is PriceCurve ... whereas I am guessing Rails thinks the singular ought to br PriceCurf ?! I''m sure there is a simple way around this linguistic blunder that doesn''t involve the use of pigeon English.

Hi, I'm trying to set up a Samba 4.1.16 file server for our (Samba 4) domain. But for some reason I can't seem to be able to get it to work. Even worse, I have run out of ideas on how to debug this. I've been tampering with this for quite some time now, but I've come to realize I can't figure this out on my own. Can someone give me some pointers on where to look in order

I just got a new server with a Dell MD-1000 SAS unit and 6-750 gigabyte drives which are now initializing in RAID 10 which will give me just about 2 terabytes. I vaguely recall reading that fdisk wasn't suitable for partitioning and wonder if I shouldn't be using partd instead. I am also wondering if I should use lvm or just mkfs to create the filesystem. Anyone have suggestions before I

We switched our mail server over from mbox to maildir a few months ago and couldn't be happier - the performance under normal load is incredible. However we now have a problem with backup.. Typically we would run tar on the mail server, sending its output to another server via ssh.. Like this: tar cfp - /home | gzip | ssh other.server (cd /backups && tar xfp -).. I have tried

Dear Community, I want to identify outliers in my data. I don't know how to use identify command in the plots obtained. I've gone through help files and use mahalanobis example for my purpose: NormalMultivarianteComparefunc <- function(x) { Sx <- cov(x) D2 <- mahalanobis(x, colMeans(x), Sx) plot(density(D2, bw=.5), main="Squared Mahalanobis distances, n=nrow(x),

Hi, R people. I'm shy reporting this, as a bug in this area sounds very unlikely. Did I make my tests wrongly? I'm still flaky at all this. Let me dare nevertheless, who knows, just in case... Please don't kill me! :-) Not so long ago, I wrote to this list: > (For now, [the library code] works only for me when I do _not_ use `-l > MY/OWN/LIBDIR' at `R CMD INSTALL'

Can anyone tell me where my blunder is in the following program? llvm-as reports: llvm-as: testit.ll:11: Can't store 'opaque *' into space of type 'opaque *'! Which doesn't seem to make sense to me. What is it that is illegal about storing a pointer to opaque in a space that is of type pointer to opaque? Is it just that you can't store pointers to opaque? %path =

On 24/03/23 2:52 pm, Rowland Penny via samba wrote: > > > On 24/03/2023 09:16, Anantha Raghava via samba wrote: >> Hi, >> >> How can I be so foolish? I just can't imagine. >> >> The issue was in smb.conf file. the the server's "netbios name" was >> reading as AD instead of DC, which was the actual host name. >> Basically I had

hi, I am still batteling with winbind, and frankly, i am getting desperate enough to long for a weekend without any computers in sight (that is bad for an addict) I know I probably make a small stuid blunder but still, I cant find it. The attachment is my smb.conf file. I created the needed paths, I even made sure the ermissions are set! I created asmbusers file with contains the following

Hi folks, sorry if this has been answered before, I searched long and hard before deciding to make a thread. I'm trying to include multiple variables in a non-parametric analysis (hah!). So far what I've managed to figure out is that the NPMC package from CRAN MIGHT be able to do what I need, but I can't get it to. First I created a dataset as NPMC calls for. > Ind=Individual

I have to disagree with both Peter and Martin on this. The underneath issue is that the automatic conversion of characters to factors by the data.frame functions was the single most egregious design blunder in the Statistical Models in S book, and we are still living with it.? The stringsAsFactors option was a compromise to let users opt out of that mistake (one I had to fight hard for).??? In

Hi all I have a hopefully not too stupid question about multi-level / mixed-effects modeling. I was trying to test a strategy from Crawley's 2013 R Book on a data set with the following structure: - dependent variable: CONSTRUCTION (a factor with 2 levels) - independent fixed effect: LENGTH (an integer in the interval [1, 61]) - random effects with the following hierarchical structure: MODE

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

I can't find a way to retrieve z values and p-values from the output from lmer in the lme4 package. How is this done? Rick B.

Hi I know that p-values doesn't appear anymore in the summary of a linear mixed-model with lmer. However, if I do a mixed logistic regression with lmer using family=binomial, the summary includes a p-values for fixed effects. Is it normal, could I use those p-values to interpret the fixed effects or should I use mcmcsamp to obtain 95% confidence interval? Thanks Julien