similar to: please help me for svm plot question

I run the following code, all other is ok, but plot(m.svm,p5.new,As~Cur) is not ok Anyone know why? install.packages("e1071") library(e1071) library(MASS) p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv") p5.new<-subset(p5,select=-Ms) p5.new$Y<-factor(p5.new$Y) levels(p5.new$Y) <- list(Out=c(1), In=c(0)) attach(p5.new)

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

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

> install.packages("e1071") Warning: package 'e1071' is in use and will not be installed > library(e1071) > library(MASS) > p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv") > attach(p5) The following object(s) are masked from p5 ( position 3 ) : Aa As Cur Ms P Y The following object(s) are

I used packages 'fields". when I load workspace again, then do installation again. I get the following error. I read R for Windows FAQ (4.8 see below), but I still can't figure it out. Does anyone can help me out? Thanks, jlm > install.packages("fields",repos="http://streaming.stat.iastate.edu/CRAN/ ",); trying URL '

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

I have a set data={A,B,C,D,E,F,G} I want to choose 2 letter from 8 letters, i.e. generate the combination set for choose 2 letters from 8 letters. I want to get the liking: combination set={AB,AC,AD,....} Does anyone konw how to do in R. thanks, Aimin

de Jong, S. (1993) SIMPLS: an alternative approach to partial least squares regression. Chemometrics and Intelligent Laboratory Systems, 18, 251?263 Thanks

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

Hello, I am new to Ubuntu (Linux in general) and I am trying to install R: so far, I am stuck. I have followed the instructions at the CRAN website with little luck. I read an exchange about the same issue posted before, but did not give me enough information to solve things by myself, hence, this email. Based on the former exchange I read, I am trying to include enough information to see if

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.zip' at the

I have a data set like this: if I want to less than 200000 obs from this data set. How can I do these? > str(p1982) 'data.frame': 465979 obs. of 6 variables: $ p : Factor w/ 1982 levels "154l_aa","1A0P_aa",..: 1 1 1 1 1 1 1 1 1 1 ... $ aa : Factor w/ 19 levels "ALA","ARG","ASN",..: 2 16 4 5 18 3 19 3 2 9 ... $ as : num 152.0

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

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

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

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.")" =

Hello fellow Samba folks, I am attempting to mount a cifs share on a RHEL 5 box using mount.cifs. The server is another RHEL 5 box. Both boxes are joined to the same Kerberos realm (AD). I kinit to get my Kerberos tickets. This is the mount command I'm using: mount.cifs //rhel5.server.iastate.edu/benvon ./mnt -o user=benvon,sec=krb5 This results in a password prompt, then a

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