search for: aimin

...ts have "P","AA",index","delta" columns I want to do histogram for variable "delta" in each data set, and label the histogram by mean, standard deviation and number of observation for "delta" column. Do you know how to do these in R? thanks Aimin Yan [[alternative HTML version deleted]]

I try to type this in my R-winEdt. but I got these. Do you know? Aimin > install.packages('http://rh-mirror.linux.iastate.edu/CRAN/bin/windows/contrib/2.4/plotrix_2.1-6.zip') Warning in download.packages(pkgs, destdir = tmpd, available = available, : no package 'http://rh-mirror.linux.iastate.edu/CRAN/bin/windows/contrib/2.4/plotrix_2.1-6...

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

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 You can now run LaTeX on '...

......................... aa have 19 levels, and there are different number of observation for each levels. I want to pick 75% of observations of each levels randomly to generate a training set, and 25% of observation of each levels to generate a testing set. Does anyone know to do this? Thanks Aimin Yan

...<- lmer(Y ~1+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1)) in addtion, I try these two models, it seems they are same. what is the difference between these two model. Is this a cell means model? m00 <- glm(sc ~aa-1,data = p5) m000 <- glm(sc ~1+aa-1,data = p5) thanks, Aimin Yan

...611 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -1.4693363 -0.8816572 -0.5038361 0.9541089 2.0939872 Number of Observations: 1030 Number of Groups: p aa %in% p 5 92 AIC,BIC,LogLik is NA. Does anybody know why? thanks, Aimin Yan [[alternative HTML version deleted]]

...knncat.train))/sum(tb.knncat.train),2) > cat("NNB prediction accuracy for training",nnb.accuracy.train,"\n") NNB prediction accuracy for training 0.67 I thought I should to get : 100%-36.88%=63.12%=0.63. But in my calculation, I get 0.67. Is my calculation wrong? thanks, Aimin

I am trying to use scatterplot function in CAR like the following: scatterplot(X~Y) I want to label X points and Y ponits using the different color. Any idea for this? Aimin

...17 20718709708 20970185038 [82] 21220574642 21469868628 21718058379 21965136456 22211096520 22455933247 22699642261 22942220059 23183663956 [91] 23423972019 23663143017 23901176371 24138072104 24373830797 24608453551 24841941947 25074298011 25305524181 [100] 25535623276 25764598466 thanks, Aimin Yan

...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 <- nnet(omega~aa_three+bas+bcu+aa_ss, data=training, size=2, linout=TRUE) # weights: 57 initial value 89153525.582093 iter 10 value 15036439.951888 iter 20 value 15010796.121891 iter 30 value 15000761.804392 iter 40 value 14955839.294531 iter 50 value 14934746....

Hello, I am new to R, How to delete a matrix column. Thanks,

if I want to put fig1plot to the left, figYPplot to the right figYAaplot on the bottom. How to modify the following cod to do these? \begin{figure}[H] \centering \begin{minipage}[t]{0.5\textwidth} \centering <<label=fig1plot,fig=TRUE,echo=FALSE>>= <<fig1plot>> @ %\caption{Caption 1} \end{minipage} \begin{minipage}[t]{0.5\textwidth} \centering

how to produce a Row Reduced Echelon Form for a matrix in R? Aimin Yan

I want to integrate R code to my WinEdt file. Can someone tell me how to do this? I copy R code to my WinEdt file, but it does't work. Aimin Yan

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 1 1 1 1 1 1 1 1 1 ... $ Aa : Factor w/ 19 levels "ALA","AR...

...63560 103.67100 inward 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 outward 4 154l_aa CYS 57.20 49.85 -0.0549589 72.97930 outward 5 154l_aa TYR 28.87 31.75 0.0526457 96.11660 inward 6 154l_aa ASN 31.14 31.09 0.0632711 55.65980 outward does anyone know how to these? Aimin > head(p1982) p aa as ms cur sc 1 154l_aa ARG 152.04 108.83 -0.1020140 92.10410 2 154l_aa THR 15.86 28.32 0.2563560 103.67100 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 4 154l_aa CYS 57.20 49.85 -0.0549589 72.97930 5 154l_aa TYR 28.87 31.75 0.052...

I try a example http://omega.psi.iastate.edu/bootstrapFile.html it doesn't give me output. I don't why. For another example http://omega.psi.iastate.edu/trivial.html it seems works except it doesn't display figure. Does anyone know how to figure it out? Aimin

...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 can give me warning like "x is used before, overwrite it or not?" when I use x again? Thanks, Aimin

I use tune.svm to tune gamma and cost for my training dataset. I use PC, it runs very slowly. Does anyone know how to make it faster? Aimin