similar to: standardized residuals

Hello, I use nnet to do prediction for a continuous variable. after that, I calculate correlation coefficient between predicted value and real observation. I run my code(see following) several time, but I get different correlation coefficient each time. Anyone know why? In addition, How to calculate prediction accuracy for prediction of continuous variable? Aimin thanks, > m.nn.omega

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

I have a dataset like this: p aa index x y z sdx sdy sdz delta as ms cur sc 1 821p MET 1 -5.09688 32.8830 -5.857620 1.478200 1.73998 0.825778 13.7883 126.91 92.37 -0.1320180 111.0990 2 821p THR 2 -4.07357 28.6881 -4.838430 0.597674 1.37860 1.165780 13.7207 64.09 50.72 -0.0977129 98.5319 3 821p GLU 3 -5.86733 30.4759

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

I am trying to fit a logistic regression model for this data set. Firstly, I want to plot P(Y=1) vs As and P(Y=1) vs Aa. Does any body know how to do these in R. Thanks, Aimin > p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv") > str(p5) 'data.frame': 1030 obs. of 6 variables: $ P : Factor w/ 5 levels "821p","8ABP",..: 1

Hi all, I am a real novice to R. :) I am struggling with a problem for generating colored dendrogram. I have searched the R list and complied/collected a R code which can generated a colored dendrogram based on the rainbow color and 4x4 similarity matrix (say matrix:m). In this dendrogram, each leaf is colored differently. But, I do not want the leaf colored on a random basis. I want to assign

I have a data(ABC) like this: x y A 3 4 A 1 3 B 2 6 B 4 8 C 5 4 C 6 7 I want to split this data into A: x y A 3 4 A 1 3 B B 2 6 B 4 8 C C 5 4 C 6 7 anyone knows how to do that? thanks, Aimin Yan

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

I have a question in R. In directory H:/Delta_angle I have 19 files like this: ALA.delta ASN.delta ASP.delta CYS.delta GLN.delta GLU.delta HIS.delta ILE.delta LEU.delta LYS.delta MET.delta PHE.delta PRO.delta SER.delta THR.delta TRP.delta TYR.delta VAL.delta I want to read these files to 19 data sets in R. All these data sets have "P","AA",index","delta"

Dear R-users, I need to standardize residuals using standard deviation. Is the 'stdres' the proper function? Beside other methods (for standardization and normalization), are there some approaches how to standardize using standard deviation? Here is the sample of my residuals: Zres040 Zres0820 Sres040 Sres0820 PCres040 -2.101740 -2.0682900 1.6328500 0.5046730

I try Sweave and get Sweave-test-1.tex but hot to run LaTeX on 'Sweave-test-1.tex'? I am using WinEdt. thanks, Aimin > Sweave(testfile) Writing to file Sweave-test-1.tex Processing code chunks ... 1 : print term verbatim 2 : term hide 3 : echo print term verbatim 4 : term verbatim 5 : echo term verbatim 6 : echo term verbatim eps pdf 7 : echo term verbatim eps pdf

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

I have a data set like this I want to assign "outward" to Y if sc <90 and assign "inward" to Y if sc>=90. then cbind(p1982,Y) to get like these p aa as ms cur sc Y 1 154l_aa ARG 152.04 108.83 -0.1020140 92.10410 inward 2 154l_aa THR 15.86 28.32 0.2563560 103.67100 inward 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 outward 4 154l_aa CYS 57.20 49.85

Dear R list, I have a question in R, it could be very simple, but I don't know how to do it? for example: I assign 6 to x in beginning of of my R script code > x<-6 ...... After many line code, I forget using x variable before, I use x again, and do assignment like this > x<-45 > x [1] 45 then value 6 of previous x is replaced by 45. I am wondering if there are some way R

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:

Use getwd(), I can get current directory name. Now I want to get current workspace name, Does anyone know how to do that? here is information about my R > sessionInfo() R version 2.4.0 (2006-10-03) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

How to add mean,sd, number of observation in each panel for lattice histogram? Aimin

I draw dotplot using the following code: sd.dotplot<-dotplot(data.47.nmr$om_sd ~ as.factor(data.47.nmr$position) |as.factor(data.47.nmr$pr), data = data.47.nmr,layout=c(1,1), xlab="Position", xlim=range(data.47.nmr$position),ylab="Sd", main="Changes of omega angle in different positions", scales = list(x = list(rot = 45))) However I get same xlim for each

In order to start matlab server in R , I using the following commands getwd() setwd("D:\R_matlab") install.packages("R.oo") install.packages("R.matlab") install.packages("R.utils") library(R.matlab) Matlab$startServer() a minimized MATLAB Command Window come out, but I can't make this window become larger. Does anyone know why? I type

Does anyone know how to solve this question about R.matlab? I am in windowsXP, my matlab is matlab 7.0.0 19920(R14) thanks, Aimin > matlab <- Matlab(host="localhost", port=9998) > if (!open(matlab)) throw("Matlab server is not running: waited 30 seconds.") Error in list("throw("Matlab server is not running: waited 30 seconds.")" =