similar to: Need help in calculating studentized residuals/leverage values of non-linear model [nls()]

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)

Hello everyone, This is my first post to the mailing list, so I hope I am posting my message the correct way. I am trying to present my dataset in a 3d scatterplot using cloud() in the {lattice} package. I hope to explicitly identify a subset of my observations. The observations in this subset need to satisfy the following critera (values on x,y, and z axes above a certain cutoff value). It

Hi all, I'm trying to apply loess regression to my data and then use the fitted model to get the *standardized/studentized residuals. I understood that for linear regression (lm) there are functions to do that:* * * fit1 = lm(y~x) stdres.fit1 = rstandard(fit1) studres.fit1 = rstudent(fit1) I was wondering if there is an equally simple way to get the standardized/studentized residuals for a

An embedded and charset-unspecified text was scrubbed... Name: nem el?rhet? URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110420/4913654b/attachment.pl>

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 =

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]] 螔伪蟻蔚胃萎魏伪蟿蔚 蟿伪 蔚谓慰蠂位畏蟿喂魏维 渭畏 未喂伪胃苇蟿蔚喂 蟿畏谓 魏伪位蠉蟿蔚蟻畏 未蠀谓伪蟿萎 蟺蟻慰蟽蟿伪蟽委伪 魏伪蟿维 蟿蠅谓 蔚谓慰蠂位畏蟿喂魏蠋谓 渭畏谓蠀渭维蟿蠅谓

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

A colleague was asking me if R does multi-level modelling as opposed to multiple regression. Since I have no knowledge of multi-level modelling (except 5 minutes googling ) I thought that I would as here. Does are offer any multi-level modeling packages? It looked like arm might be one but I was not sure. Thanks

Sorry, I am sure, this must be documented somewhere (but there are that many docs and tutorials to scan for topics..., actually a great thing... but if you are in a hurry..): I want to save a plot as .ps (or .eps): > postscript("plot1.eps", horizontal=FALSE, onefile=FALSE,height=8,width=8,pointsize=10) > plot(hpfit$fit,rstudent(hpfit),xlab="Fitted

Inline below: 2011/4/21 Jeremy Miles <jeremy.miles at gmail.com> > > Just because it comes from a book does not make it true or correct. Amen! > Books are subject to considerably less peer review than journal > articles. Yes, but ... Peer review among journals is uneven, especially for those from private for-profit publishers. And even for top flight journals, dealing with

Hello, I have 4 questions in R and I am looking for 4 answers! the questions are the below: 1) is there a function in R for the Least Significant Difference method for multiple comparisons? 2) Is there a function in R that gives me the studentized deleted residuals and the leverage points? 3) is there a way that I could create a link in a program of mine, so as to select an option when my program

Hi, Im trying to compile kernel 3.0.4 inside Centos 5 DomU The steps I do are Make bzImage Make modules Make modules_install Depmod -a mkinitrd /boot/initrd.img-3.0.4 3.0.4 but them im getting this -bash-3.2# mkinitrd /boot/initrd.img-3.0.4 3.0.4 ls: /etc/modprobe.d/*.conf: No such file or directory No module ehci-hcd found for kernel 3.0.4, aborting. If I exclude the

Hey everyone, I have a question about vif-common.sh. I run multiple bridges attached on dummy interfaces, which allow me to put guests in seperate subnets (routed through the dom0). As you might expect I already have quite extensive iptables scripts to accomidate this kind of routing. I was just hoping someone on this list can confirm, that I understand what the iptables lines in vif-common.sh

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;

I am using the Unix version of R (version 1.7.0), installed via fink on a G4 Macintosh. I recently upgraded from version 1.6.0 and found that the "car" library now has a problem: ---Begin transcript--- >library(car) Attaching package 'car': The following object(s) are masked from package:base : dfbeta dfbeta.lm dfbetas dfbetas.lm hatvalues hatvalues.lm

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

Hi, For the past day, I have been researching about making dovecot to support mailfront (http://untroubled.org/mailfront/) imapfront-auth natively and able to use mailfront imapfront-auth with dovecot imap v1.2.16 without starting any dovecot related services and no /etc/dovecot.conf file is required. More information about this can be read at:

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 <-

Hey hope this opne is a "make senser". Has anyone tried the "Subject line" combo?, I''ve compiled a dozen src.rpm kernel 3.0 and 3.1 packages with no luck to try native upstream Dom0 kernel support with no success, just a kernel oops and constant reboots. I followed Fedora''s wiki on building a custom kernel, even used F16''s .config which also was a

Hey hope this opne is a "make senser". Has anyone tried the "Subject line" combo?, I''ve compiled a dozen src.rpm kernel 3.0 and 3.1 packages with no luck to try native upstream Dom0 kernel support with no success, just a kernel oops and constant reboots. I followed Fedora''s wiki on building a custom kernel, even used F16''s .config which also was a