search for: aiminy

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) m.svm<-svm(Y~P+Aa+As+Cur,data=p5.new) summary(m.svm) plot(m.svm,p5.new,As~Cur) Here is output: > install.packages("e1071") --- Please...

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) m.svm<-svm(Y~P+Aa+As+Cur,data=p5.new) summary(m.svm) plot(m.svm,p5.new,As~Cur) Here is output: > install.packages("e1071") --- Please...

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","ARG","ASN",..: 12 16 7 18 11 10 19 19 19 1 ... $ As :...

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

...ntent type 'application/zip' length 4200608 bytes (4.0 Mb) opened URL downloaded 4.0 Mb package 'fields' successfully unpacked and MD5 sums checked Warning: cannot remove prior installation of package 'fields' The downloaded packages are in E:\Documents and Settings\aiminy\Local Settings\Temp\RtmpwsmfG9\downloaded_packages updating HTML package descriptions > library(fields); Error in library(fields) : there is no package called 'fields' 4.8 update.packages() fails. You may not be able to update a package which is in use: Windows `locks' the package...

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

I am trying someone's example to build a R package, but I get some errors, Can someone help me to figure it out? thanks, * using log directory '/home/aiminy/myrpackages/roots.Rcheck' * using R version 2.6.1 (2007-11-26) * checking for file 'roots/DESCRIPTION' ... OK * this is package 'roots' version '0.1-1' * checking package dependencies ... OK * checking if this is a source package ... WARNING Subdirectory 'roots/src...

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

> 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 masked from p5 ( position 4 ) : Aa As Cur Ms P Y > m.svm<-svm(as.factor(Y)~.,data=p5) > summary(m.sv...

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

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

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

...e de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush dts acpi mmx fxsr sse sse2 ss ht tm syscall lm pni monitor ds_cpl cid cx16 xtpr ts bogomips : 5537.79 clflush size : 64 cache_alignment : 128 address sizes : 36 bits physical, 48 bits virtual power management: [aiminy at omega ~]$ cat /proc/meminfo MemTotal: 2058192 kB MemFree: 311068 kB Buffers: 924 kB Cached: 15216 kB SwapCached: 176776 kB Active: 1543584 kB Inactive: 153280 kB HighTotal: 0 kB HighFree: 0 kB LowTotal: 2058192 kB LowFree:...

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, I am new to R, How to delete a matrix column. Thanks,