Displaying 20 results from an estimated 20000 matches similar to: "standardized residuals"
2007 Jan 28
2
nnet question
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
2010 Nov 10
1
standardized/studentized residuals with loess
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
2006 Nov 20
4
for help about logistic regression model
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
2011 Mar 16
1
Standardized Pearson residuals (and score tests)
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
2006 Nov 26
1
plot p(Y=1) vs as
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
2006 Dec 26
1
Colored Dendrogram
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
2007 Jan 19
2
split data set
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
2004 Apr 05
0
studentized deleted residuals and NA's
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
2006 Oct 24
2
for help
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"
2005 Jan 18
0
standardised residuals using standard deviation
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
2006 Nov 23
2
Sweave question
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
2011 Mar 14
3
Standardized Pearson residuals
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
2006 Dec 08
1
question for if else
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
2007 Mar 07
1
how to avoid to overwrite object
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
2005 Dec 06
1
standardized residuals (rstandard & plot.lm) (PR#8367)
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:
2007 Mar 04
2
lattice histogram
How to add mean,sd, number of observation in each panel for lattice histogram?
Aimin
2007 Mar 24
1
how to get current workspace name?
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
2007 Apr 24
1
specify xlim for each panel
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
2006 Dec 17
1
Start Matlab server in R 2.4.0
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
2006 Dec 20
1
R.matlab question
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.")" =