search for: rstudent

...w_chla=1/chla/chla; run; ODS listing close; %macro sort_event(cut_off,last=0); /*proc glm data=sort_dataset; class year; model logchla=year cos1 sin1 cos2 sin2 cos3 sin3 cos4 sin4 /solution; by station; where bloom=0; output out=chla_res predicted=pred student=studres cookd=cookd rstudent=rstudent u95=u95; lsmeans year / at (cos1 sin1 cos2 sin2 cos3 sin3 cos4 sin4)=(0 0 0 0 0 0 0 0); ODS output ParameterEstimates=parmest LSmeans=lsmeans; run;*/ proc glm data=sort_dataset; class year month; model chla=/solution; by station; weight w_chla; where bloom=0; output...

...hat 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 influence influence.glm influence.lm rstudent rstudent.glm rstudent.lm >x <- 1:10 >y <- 2*x+3+2*rnorm(10) >model<-lm(y~x) >qq.plot(model) Bus error ---End transcript--- No plot appears. After the bus error, R no longer operates normally, and I have to kill the terminal window in order to quit. i just upgraded from R v...

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

Hi, I am trying to figure out how to get the correlation coefficient for a QQ plot (residual plot). So to be more precise, I am creating the plot like this: qq.plot(rstudent(regrname), main = rformula, col=1) But want to also access (or compute) the correlation coefficient for that plot. Thanks, Tom [[alternative HTML version deleted]]

...ectively. 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 = function(dataset,dependent,independent){ dataset$predicted = predict(lm(dependent~independent)) dataset$stdres = rstudent(lm(dependent~independent)) m = 1 for(i in 1:length(dataset$stdres)){ dataset$outlier_counter[i] = if(dataset$stdres[i] >= 3 | dataset$stdres[i] <= -3) {m} else{0} } j = length(which(dataset$outlier_counter >= 1)) while(j>=1){ print(dataset[which(dataset$o...

Command I am using for the plot: plot(Raw[][Plate==101]~well[][Plate==101], xlab="Well", ylab="Raw", main="Plate 101") I only want to label points on the graph where well equals B13, G13, K14 and N14 with the name of the well. Thank you for your help. -- View this message in context:

...ed 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 Response",ylab="Studentized Residuals") unfortunately, I cannot find the axes and labels in the .ps anymore. where are they gone? how to fix it? thanks a lot christoph -- Christoph Lehmann Phone: ++41 31 930 93 83 Department o...

...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) residuals(mod)/(summary(mod)$sigma*sqrt(1-lm.influence(mod)$hat)) -- View this message in context: http://r.789695.n4.nabble.com/studentized-and-standarized-residuals-tp3732997p3732997.html Sent from the R help mailing list archive at Nabble.com.

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 loess model? BTW my apologies if there is something here that I'm missing. All the best, * * *Oliver * [[alternative HTML version deleted]]

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 17 rows, 7 of them are "new rows", where all the entries are NA's, even its rownames. Because of this I'm not sure of the dimension of a[rstu >3, ] (Do...

I just started using R last week. I have a dataset with 3 columns - Plate, Well and Raw I need to make a simple plot(Well~Raw) but only when Plate = 101 Thanks for your help. -- View this message in context: http://r.789695.n4.nabble.com/Simple-conditional-plot-tp3654300p3654300.html Sent from the R help mailing list archive at Nabble.com.

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>

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

....lm* formula.lm* hatvalues.lm [17] influence.lm* kappa.lm labels.lm* logLik.lm* [21] model.frame.lm model.matrix.lm nobs.lm* plot.lm [25] predict.lm print.lm proj.lm* qr.lm* [29] residuals.lm rstandard.lm rstudent.lm simulate.lm* [33] summary.lm variable.names.lm* vcov.lm* Non-visible functions are asterisked > > library(gdata) gdata: read.xls support for 'XLS' (Excel 97-2004) files ENABLED. gdata: read.xls support for 'XLSX' (Excel 2007+) files ENABLED. Attaching p...

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

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

..., and it seems that it returns elements sigma and coefficients that are only based on the first variable (first column of the residual matrix wt.res) and give wrong results for other variables. This will influence functions dfbeta.lm (coefficients), dfbetas.lm (coefficients, sigma), dffits (sigma), rstudent.lm (sigma) and covratio (sigma). lm.influence finds these elements in compiled code and this is harder to fix. MASS (the book & the package) avoid using compiled code in their (univariate) studentized residuals, and instead use a clever short-cut. In addition to these, there are a couple of ot...

Hi, P-values in Durbin-Watson test obtained through the use of functions available in packages "lmtest" and "car" are different. The difference is quite significant. function "dwtest" in "lmtest" is much faster than "burbinwatson" in "car". Actually, you can take a nap while the latter trying to calculated Durbin-Watson test. My question

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

Am I asking for too much: for any object that a stat proc returns ( y <- lm( y~x) , etc ) ) , is there a super convenient function like give_all_extractors( y ) that lists all extractor functions , the datatype returned , and a text descriptor field ("pairwisepval" "lsmean" etc) That would just be so convenient. What are my options for querying an object so that I can