similar to: standardized/studentized residuals with loess

Dear R-Help, I am using the studres function from the MASS package to compute studentized deleted residuals in a oneway anova. I'm having trouble interpreting the results in situations where a factor level has only one observation. Sometimes studres yields an NaN and sometimes it produces a numeric value for cases where a factor level has only one observation. I would think it should

Hi, I must be doing something silly here, because I can't get the studentised and standardised residuals from r output of a linear model to agree with what I think they should be from equation form. Thanks in advance, Jennifer x = seq(1,10) y = x + rnorm(10) mod = lm(y~x) rstandard(mod) residuals(mod)/(summary(mod)$sigma) rstudent(mod)

Is there any reason that rstandard.glm doesn't have a "pearson" option? And if not, can it be added? Background: I'm currently teaching an undergrad/grad-service course from Agresti's "Introduction to Categorical Data Analysis (2nd edn)" and deviance residuals are not used in the text. For now I'll just provide the students with a simple function to use, but I

Dear friends, In s-plus, lm() generates the the studentized residuals automatically for us, and In R, it seems don't have the results: After i fitted lm(), i use attibutes() to see the objects and didn't find studentized residuals . How to get the the studentized residuals in lm(),have i missed something? thanks very much! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of Epidemiology

I contacted John Fox about this first, because parts of the file are attributed to him. He says that he didn't write rstandard.glm(), and suggests asking r-devel. As it stands, rstandard.glm() has summary(model)$dispersion outside the sqrt(), while in rstandard.lm(), the sd is already sqrt()ed. This seems to follow stdres() in VR/MASS/R/stdres.R. Of course for the c("poisson",

I am working on a function that will remove outliers for regression analysis. I am stating that a data point is an outlier if its studentized residual is above or below 3 and -3, respectively. The code below is what i have thus far for the function x = c(1:20) y = c(1,3,4,2,5,6,18,8,10,8,11,13,14,14,15,85,17,19,19,20) data1 = data.frame(x,y) rm.outliers =

Hi, I wonder if anybody could help me in converting this easy SAS program into R. (I'm still trying to do that!) PROC IMPORT OUT= WORK.CHLA_italian DATAFILE= "C:\Documents and Settings\carleal\My Documents\REBECCA\stat\sas\All&nutrients.xls" DBMS=EXCEL2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian;

Hi there, I hope I can get advice regarding the calculation of leverage values or studentized residual values of a non-linear regression model. It seems like rstudent() does not work on a nls object. Many thanks in advance! Best regards, Xingli

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

Full_Name: Heather Turner Version: 2.2.0 OS: Windows XP Submission from: (NULL) (137.205.240.44) Standardized residuals as calculated by rstandard.lm, rstandard.glm and plot.lm are Inf/NaN rather than zero when the un-standardized residuals are zero. This causes plot.lm to break when calculating 'ylim' for any of the plots of standardized residuals. Example:

Dear all, the application of the function rstandard() in the base package to a glm object does not produce residuals standardized to have variance one: the reason is that the deviance residuals are divided by the dispersion estimate and not by the square root of the estimate for the dispersion. Should the function not be changed to produce residuals with a variance about 1? R 1.8.1 on

I am trying to find studentized bootstrap intervals for the skewness of a data set. I have tried the following (nerve.dat is a set of data containing observations on one variable) (using Windows XP): x <- scan("e:/Flashbackup2009/Nonparametrics/nerve.dat") n <- length(x) library(e1071) skewness(x) library(boot) sampleskew <- function(x,d) {return(skewness(x[d]))} bb <-

Hi Peter and others, If it helps, I wrote a small function glm.scoretest() for the statmod package on CRAN to compute score tests from glm fits. The score test for adding a covariate, or any set of covariates, can be extracted very neatly from the standard glm output, although you probably already know that. Regards Gordon --------------------------------------------- Professor Gordon K

Dear colleagues, hi. I am using the function gls in R-package nlme and I want to compute either the leverage values or the studentized residuals. Can any body tell me how can do that? With regards ___________________________________________________________ 围蟻畏蟽喂渭慰蟺慰喂蔚委蟍[elided Yahoo spam]] 螔伪蟻蔚胃萎魏伪蟿蔚 蟿伪 蔚谓慰蠂位畏蟿喂魏维 渭畏 未喂伪胃苇蟿蔚喂 蟿畏谓 魏伪位蠉蟿蔚蟻畏 未蠀谓伪蟿萎 蟺蟻慰蟽蟿伪蟽委伪 魏伪蟿维 蟿蠅谓 蔚谓慰蠂位畏蟿喂魏蠋谓 渭畏谓蠀渭维蟿蠅谓

Well, it is clear enough that the problem is in interpreting the documentation. However, when you claim you tested something, and found it inconsistent with tables, it would be advisable to back it up with examples! The description in the help files and in the sources is admittedly confusing. The original paper has this, rather more clear, description in the abstract: "We consider the

We cannot help you understand what you are doing if you do not show us what you are doing. Here are some discussions about how to communicate questions about R [1][2][3]. [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example [2] http://adv-r.had.co.nz/Reproducibility.html [3] https://cran.r-project.org/web/packages/reprex/index.html -- Sent from my phone.

Dear all, I wanted to compare Bonferroni vs TukeyHSD correction over a range of groups and group sizes, and wanted to use the function qtukey. In the help documentation it says qtukey(p, nmeans, df, nranges = 1, lower.tail = TRUE, log.p = FALSE) Arguments q vector of quantiles. p vector of probabilities. nmeans sample size for range (same for each group). df degrees of freedom for s (see

Hi, I am newbie user of ggplot and would like some assistance in implementing time series plots. I'd like to know how the tsdiag plot can be made in ggplot? Thanks Harsh Singhal Decisions Systems, Mu Sigma Inc.

Dear list members, I want to save residuals above or less than a certain value to an R object. I have performed a multiple linear regression, and now I want to find out which cases have a residual of above + 2.5 and ? 2.5. Below I provide the R commands I have used. Reg<-lm(a~b+c+d+e+f) # perform multiple regression with a as the dependent variable. Residuals<-residuals(reg) # store

Greetings to r-help land. I've run into some program crashes and I've traced them back to methods() behavior after the package gdata is loaded. I provide now a minimal re-producible example. This seems bugish to me. How about you? dat <- data.frame(x = rnorm(100), y = rnorm(100)) lm1 <- lm(y ~ x, data = dat) methods(class = "lm") ## OK so far library(gdata)