similar to: How to Save the residuals of an LM object greater or less than a certin value to an R object?

Residuals are stored as a numeric vector. The R software comes with a document "Introduction to R" that discusses basic math functions and logical operators that can create logical vectors: abs( stdresiduals ) > 2.5 It also discusses indexing using logical vectors: stdresiduals[ abs( stdresiduals ) > 2.5 ] Note that in most cases it is worth going the extra step of making your

Hi all, I want to make a contingency table in R. I want to tabulate two variables, one as the independent and second as the dependent variable. The IV has two categories, namely, birth complications, and no birth complications. The frequency of birth complication category is fifty, and the frequency of no birth complication category is 34. The categories and frequencies of DV follows.

I am trying to perform goodness of fit test using R. I am using this website http://wiener.math.csi.cuny.edu/Statistics/R/simpleR/stat013.html for help. However, I am unable to carry out the test successfully. My code follows. It is taken from the website just mentioned. freq=c(22,21,22,27,22,36) # frequencies obtained after rolling the dice 150 times. prob=c(1,1,1,1,1,1)/6 # specify expected

Hi Peter, the "residuals()" function returns the residuals of a model fitted using the "lm" function. For instances, using the example included in the help of lm: ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <- gl(2, 10, 20, labels = c("Ctl","Trt")) weight <-

Dear list, I have the following how-to-do it in R, questions. Suppose I have ten independent variables, and one dependent variable. I want to find the Pearson correlation of all the IVs with the DV, but not the correlation between the IVs. What I know so far, about R, that I have to type the cor () function ten times, each time requesting for a correlation between one IV and the DV. I was

Hi all, I would like to read the results of the commands type in the terminal window in Microsoft Word. As a blind user my options are somewhat limited and are time consuming if I want to see the results of the commands that I have type earlier. for example if my first two commands were x<-c(1,2,3,4,5) mean(x) and I have typed ten more commands after the first two commands it is not easy

Dear list, I am slightly confused about how I can do the following in R. I want to perform robust multiple regression. I?ve used the Boot function in CAR package to find confidence intervals and standard errors. Inadition to these, I want to find the robust estimates for the F test and r-square. Finally, I would like to know the significance levels of bootstrap results. Below I explain

Hi For first question, maybe I am completely wrong but cor(swiss[,-1], swiss[,1]) should give you what you want in one step. Second question Without an example it is hard to say but maybe aggregate is the way forward. > aggregate(iris[,1:4], list(iris$Species), function (x) c(mean=mean(x), sd=sd(x))) Group.1 Sepal.Length.mean Sepal.Length.sd Sepal.Width.mean Sepal.Width.sd 1

Hi all, I am trying to use type message with sink, like this: sink("all.Rout", type="message") 1+3 sink() readLines(con = "all.Rout") So to get the following output: > 1+3 [1] 4 Obviously this doesn't work. I tried some variations (based on the explanations in the help) but am missing something on how to make it work. Any suggestions? (p.s: I need

After reading more, I understand I didn't formulate my last question correctly, so please allow me to rephrase: What I am looking for is a way to save the R console session output. That is, a command that would combine the results of using: ?sink # And ?savehistory My motivation for this is that doing it will allow someone who is a blind user of R to be able to easily export his results to

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:

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

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

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)

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

Hi all, I've got some questions to post in regard to running asterisk in a production-grade environment, specifically targeting high-density IVR applications. No VoIP involved, just straight PSTN -> * and perhaps the occasional outdials or agent-based predictive dialing. 1) Which user would you run * under? 2) What other security-related issues do you have to resolve? 3) How do you handle

Hi, I''m looking for past efforts to integrate Rails and Needle. I am thinking of doing it myself but would be interested in seeing what others have done before. If there are examples of other IoC frameworks being used with Ruby I''d be interested in hearing about them too. I find IoC to be a helpful concept in general and, having glanced at a few existing Rails app.s, I

Hello, I''m experimenting with a xhtml templating language and have some questions about ActionPack, particularly to do with compilable templates. I''m mostly looking for further information about compilable templates. The templating language currently works, though it isn''t a compilable template handler. I''ve yet to benchmark it, but it seems to be running

This bug is not quite fixed - the example from my original report now = works using R-2.2.1, but plot(Uniform, 6) does not. The bug is due to if (show[6]) { ymx <- max(cook, na.rm =3D TRUE) * 1.025 g <- hatval/(1 - hatval) # Potential division by zero here # plot(g, cook, xlim =3D c(0, max(g)), ylim =3D c(0, ymx),=20 main =3D main, xlab =3D