Re: [R] Help with SVM package Kernlab

2009-12-24 Thread Steve Lianoglou
Hi, On Fri, Dec 25, 2009 at 12:49 AM, Vishal Thapar wrote: > Hi Steve, > > Thank you so much for the reply. The response to your queries are: > What do these commands return over your data? > > 1. is(train500) > -->"data.frame" "list"   "oldClass"   "mpinput"    "vector" > 2. is(train500$clas

Re: [R] Help with SVM package Kernlab

2009-12-24 Thread Vishal Thapar
Hi David, Thanks for your reply. The package is Kernlab. I agree with you when you mention that R didn't like the structure of the testSeq object but I can't figure out why. Here is the output for "dput(testSeq) and str(testSeq). Any insight that you can see would be really helpful to me. Thanks

Re: [R] Help with SVM package Kernlab

2009-12-24 Thread Vishal Thapar
Hi Steve, Thank you so much for the reply. The response to your queries are: What do these commands return over your data? 1. is(train500) -->"data.frame" "list" "oldClass" "mpinput""vector" 2. is(train500$class) --> "NULL" "OptionalFunction" "output" 3. is(train500[1,5])

Re: [R] help in merging

2009-12-24 Thread utkarsh . singhal
Thanks, but I kind of new this way already and it doesn't seem an optimal thing to do. What I was looking for is to pass some argument in 'merge' itself which doesn't change the ordering of 'x'. Or, more than that, I am interested in knowing that why is it changing the o

Re: [R] Help with SVM package Kernlab

2009-12-24 Thread David Winsemius
On Dec 24, 2009, at 11:42 PM, Vishal Thapar wrote: Hi useR's, I am resending this request since I got no response for my last post and I am new to the list so pardon me if I am violating the protocol. I am trying to use the "Kernlab" package for training and prediction using SVM's. I am

Re: [R] Help with SVM package Kernlab

2009-12-24 Thread Steve Lianoglou
Hi, Comments in line: On Thu, Dec 24, 2009 at 11:42 PM, Vishal Thapar wrote: > Hi useR's, > > I am resending this request since I got no response for my last post and I > am new to the list so pardon me if I am violating the protocol. > > I am trying to use the "Kernlab" package for training and

[R] problem with .Rhistory: please ignore

2009-12-24 Thread Erin Hodgess
There was a "permissions" problem on .Rhistory All is well now. Sorry, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodg...@gmail.com __ R-help@r-project.org

[R] problem with .Rhistory

2009-12-24 Thread Erin Hodgess
Dear R People: I just updated to R-2.10.1 on an Ubuntu Karmic Koala and am getting the following with saving .Rhistory e...@erin-desktop:~$ R R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTE

[R] Help with SVM package Kernlab

2009-12-24 Thread Vishal Thapar
Hi useR's, I am resending this request since I got no response for my last post and I am new to the list so pardon me if I am violating the protocol. I am trying to use the "Kernlab" package for training and prediction using SVM's. I am getting the following error when I am trying to use the pred

Re: [R] checking package errors-'require' calls not declared and no visible global function definition

2009-12-24 Thread Gabor Grothendieck
Don't list c in the NAMESPACE file or else you can use the .Rbuildignore file to completely eliminate any files in your sources tree from the built source. On Thu, Dec 24, 2009 at 10:31 PM, rusers.sh wrote: >  Thanks. >  BTW, if i have three functions, a,b and c. I only want the two functions > "

Re: [R] checking package errors-'require' calls not declared and no visible global function definition

2009-12-24 Thread rusers.sh
Thanks. BTW, if i have three functions, a,b and c. I only want the two functions "a" and "b" to be exported for use and "c" is only used by myself because it is not very stable. How to set this? Thanks a lot. 2009/12/24 Gabor Grothendieck > The first one is a warning and the second is a no

Re: [R] checking package errors-'require' calls not declared and no visible global function definition

2009-12-24 Thread Gabor Grothendieck
The first one is a warning and the second is a note. You did not get any errors. On Thu, Dec 24, 2009 at 9:40 PM, rusers.sh wrote: > Hi, >  I got two errors after i checking my package. The errors are as follows, > Error1: > * checking for unstated dependencies in R code ... WARNING > 'library'

[R] checking package errors-'require' calls not declared and no visible global function definition

2009-12-24 Thread rusers.sh
Hi, I got two errors after i checking my package. The errors are as follows, Error1: * checking for unstated dependencies in R code ... WARNING 'library' or 'require' calls not declared from: boot np Error2: * checking R code for possible problems ... NOTE myfunction: no visible global functio

Re: [R] Question to use R plot GO pie chart

2009-12-24 Thread Martin Morgan
Waverley @ Palo Alto wrote: > Hi, > > I have a list of IPI gene IDs. I want to find out whether there is a > package which can map the gene ontology to these IPIs, and plot the > pie chart to demonstrate the molecular function distributions. > > The input is like the following gene IPI IDs: > IP

Re: [R] How to separate a data set by its factors

2009-12-24 Thread Carl Witthoft
G Quote: " I am climbing up a steep R learning curve, and so would appreciate some help." Please help stamp out technical illiteracy: a steep learning curve is a GOOD thing. Think about what a "learning curve" represents. Time is the x-axis and "knowledge" or 'ability' is the y

Re: [R] An unprofessional message

2009-12-24 Thread L.A.
Well, a beginner just learning here, so >ToYou<-"MERRY CHRISTMAS" >print(ToYou) L.A. -- View this message in context: http://n4.nabble.com/An-unprofessional-message-tp978419p978670.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] An unprofessional message

2009-12-24 Thread Peter Ehlers
Oh my, a useR for 1.5 years and still sending HTML. Tsk, tsk. -Peter Ehlers Maithili Shiva wrote: Dear R helpers, � I understand that this is absolutely unprofessional on my part and this group doesn't entertain such things. I have been associted with this group since�last 1 and half years a

[R] Question to use R plot GO pie chart

2009-12-24 Thread Waverley @ Palo Alto
Hi, I have a list of IPI gene IDs. I want to find out whether there is a package which can map the gene ontology to these IPIs, and plot the pie chart to demonstrate the molecular function distributions. The input is like the following gene IPI IDs: IPI:IPI8860.1|SWISS-PROT:Q9BXJ4-1|TREMBL:Q

Re: [R] An unprofessional message

2009-12-24 Thread Ivan Gregoretti
Great animation, really amusing. Merry Christmas gone Psychedelic Christmas. Ivan Ivan Gregoretti, PhD National Institute of Diabetes and Digestive and Kidney Diseases National Institutes of Health 5 Memorial Dr, Building 5, Room 205. Bethesda, MD 20892. USA. On Thu, Dec 24, 2009 at 5:55 PM,

Re: [R] An unprofessional message

2009-12-24 Thread Yihui Xie
I have wrapped it into a Flash animation, so we can really see it now: http://yihui.name/en/2009/12/merry-christmas-using-r/ Regards, Yihui -- Yihui Xie Phone: 515-294-6609 Web: http://yihui.name Department of Statistics, Iowa State University 3211 Snedecor Hall, Ames, IA On Thu, Dec 24, 2009

[R] How to dbReadTable() only a limited number of rows? (RMySQL)

2009-12-24 Thread Peng Yu
I only want to load a limited number of rows by dbReadTable(). I don't see an option in the help. Is there an option to do so? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] How to separate a data set by its factors

2009-12-24 Thread David Winsemius
On Dec 24, 2009, at 3:24 PM, James Rome wrote: I have a large data set of airport data and wish to analyze it by hour and day of the week. hour and day of the week are factors. I can do something such as: histogram(~() | , type="count", breaks=60) which displays the data the way I want it in p

[R] lapack error on Redhat Linux Fedora 11

2009-12-24 Thread Hongbin Zhang
Hi, The following errors occurs when I run the nlme's example: -- Error in chol.default((value + t(value))/2) : lapack routines cannot be loaded In addition: Warning message: In chol.default((value + t(value))/2) : unable to lo

Re: [R] How to separate a data set by its factors

2009-12-24 Thread Steve Lianoglou
Hi, On Thu, Dec 24, 2009 at 3:24 PM, James Rome wrote: > I think my question boils down to "how do you replace a whole data set > by its factored subsets in all of the usual R commands? I think the answer to your question is: I'm not sure that there's a way to do that in "all of the usual R com

[R] How to separate a data set by its factors

2009-12-24 Thread James Rome
I have a large data set of airport data and wish to analyze it by hour and day of the week. hour and day of the week are factors. I can do something such as: histogram(~(Arrival.Val) | DAY*Hour, type="count", breaks=60) which displays the data the way I want it in principle, but the plots are too

Re: [R] help in merging

2009-12-24 Thread milton ruser
Hi there, You can add a order column on x or y and after use this field to order z. Like z<-z[order(z$orderfield),] To generate a order on x or y you can do something like x$xorder<-1:nrow(x) cheers milton On Thu, Dec 24, 2009 at 2:26 PM, wrote: > > Hi All, > I want to "merge" two datase

[R] help in merging

2009-12-24 Thread utkarsh . singhal
Hi All, I want to "merge" two datasets by column "ID" and I don't want the result to be sorted by "ID". I am doing the following: > z = merge(x, y, by = "ID", sort=F) The result is not sorted by "ID". But (as oppose to what I expected) it is not even in the original order of eith

[R] Multiple CHOLMOD errors when attempting poisson glmm

2009-12-24 Thread postava-davig.m
Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod

Re: [R] An unprofessional message

2009-12-24 Thread Yihui Xie
As a useR, you'd better present this message in a more professional way :) Let's see Christmas shining in R: n = length(speed <- runif(angle <- runif(x <- strsplit("MERRY CHRISTMAS", "")[[1]], 0, 360), 0, 15)) x11(10, 3) par(mar = rep(0, 4), bg = "black") for (j in 1:1000) { angle = angle + sp

Re: [R] Running 32 bit R in terminal on Mac OS

2009-12-24 Thread David Winsemius
On Dec 24, 2009, at 1:33 PM, David Winsemius wrote: Try this in an R session: if(.Machine$sizeof.pointer == 4){print("32 bit")} else {print("64 bit")} Dropped an "e" somehow... if(.Machine$sizeof.pointer == 4){print("32 bit")} else {print("64 bit")} # for me it's [1] "64 bit" On De

Re: [R] Running 32 bit R in terminal on Mac OS

2009-12-24 Thread David Winsemius
Try this in an R session: if(.Machine$sizeof.pointer == 4){print("32 bit")} else {print("64 bit")} # for me it's [1] "64 bit" On Dec 24, 2009, at 1:04 PM, saleem1000 wrote: Hello, I am using Mac OS on R. When I start R thru the application menu What's the "application menu"? How di

Re: [R] as.Date question

2009-12-24 Thread Gabor Grothendieck
zzz1 is POSIXct so looking at: > as.Date.POSIXct function (x, ...) { z <- floor(unclass(x)/86400) attr(z, "tzone") <- NULL structure(z, class = "Date") } we see as.Date.POSIXct takes the POSIXct object, zzz1, and converts it to Date relative to GMT. There is no time zone argument on

[R] Running 32 bit R in terminal on Mac OS

2009-12-24 Thread saleem1000
Hello, I am using Mac OS on R. When I start R thru the application menu I have the option of running both the 32 bit version as well as the 64 bit version. When I type in R in terminal I believe it starts up the 64 bit version. How can I start the 32 bit version instead?? Saleem -- View this

Re: [R] how to do multiple responses in a linear regression

2009-12-24 Thread Benilton Carvalho
you need to be more clear on your question... what is it (exactly) that you want? Is it the following? y1 ~ x1 + ... + xm y2 ~ x1 + ... + xm ... yn ~ x1 + ... + xm ? if so: lm(cbind(y1, y2, ..., yn) ~ x1+x2+...+xm) b On Dec 24, 2009, at 3:33 PM, Hao Cen wrote: > Hi, > > I have multiple res

Re: [R] how to do multiple responses in a linear regression

2009-12-24 Thread Jorge Ivan Velez
Dear Jeff, Take a look at the following example. It may get you started. # Some data set.seed(123) x1 <- rnorm(100) x2 <- runif(100) x3 <- rpois(100, 2) er <- rnorm(100) y1 <- 3 + 2*x1 + .5*x2 + x3 + er y2 <- 1 + .5*x1 + 3*x2 + .2*x3 + er # model m <- lm(cbind(y1, y2) ~ x1 + x2 + x3) summary(m)

[R] how to do multiple responses in a linear regression

2009-12-24 Thread Hao Cen
Hi, I have multiple responses y1, y2, .., yn, and would like to do linear regression for each of them with x1, x2, ..., xm. Instead of doing regression n times, it it possible to do it all at once? I tried lm(y1+y2 ~ x1 + x2 + x3) and lm added y1 y2 and then did the regression. thanks Jeff ___

Re: [R] as.Date question

2009-12-24 Thread Linlin Yan
I am afraid that although in same literally, they are indeed different functions: as.Date.POSIXct and as.Date.POSIXlt. But I am not sure why they are designed like this, which causes the confusion as you mentioned. On Thu, Dec 24, 2009 at 11:02 PM, MAL wrote: > Mark, not sure that's the answer. >

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
may be ?table ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask wrote: > Dear list > > I have a response variable coded 0/1 i.e. a binary response. There are > 20,000 individual responses that I would like to aggregate into numbers of > each category (i.e. 0/1) by > group called dn (350 different g

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
?aggregate ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask wrote: > Dear list > > I have a response variable coded 0/1 i.e. a binary response. There are > 20,000 individual responses that I would like to aggregate into numbers of > each category (i.e. 0/1) by > group called dn (350 different grou

[R] aggregate binary response data

2009-12-24 Thread Graham Leask
Dear list I have a response variable coded 0/1 i.e. a binary response. There are 20,000 individual responses that I would like to aggregate into numbers of each category (i.e. 0/1) by group called dn (350 different groups) and by month mth (there are several hundred responses per month. W

Re: [R] Column naming issues using read.table

2009-12-24 Thread John Kane
I'm not exactly an expert so this is not likely a good way to do it but if the actual variable names are constant across the files why not just read in the data and assign the names later? see skip in ?read.table. x <- read.table("d:/junk1.txt", skip=2) should read in the data. You can cl

Re: [R] as.Date question

2009-12-24 Thread MAL
Mark, not sure that's the answer. Usually one has x=y --> f(x)=f(y) which doesn't seem to hold here (put x=zzz1, y=zzz2, f=as.Date()). Or do I overlook something? - Original Message - From: "Marek Janad" To: Sent: Thursday, December 24, 2009 00:08 Subject: Re: [R] as.Date question

Re: [R] Two Easy questions

2009-12-24 Thread John Kane
?text for the first question ? plot(xx$Eight, xx$Punctuation) text(xx$Eight, xx$Punctuation, xx$Name) You will need to play around with the values for Punctuation to get the names to be beside the dots --- On Thu, 12/24/09, Jose Narillos de Santos wrote: > From: Jose Narillos de Santos >

Re: [R] Two Easy questions

2009-12-24 Thread Cedrick W. Johnson
Jose Narillos de Santos wrote: Second, I have read that with R we can dowload directly stock market historical information from google or yahoo. Can anyone guide me how to do it or more easily where I can find information (a web link or something similar) Thanks for all¡¡¡ Check out getSymbol

[R] Two Easy questions

2009-12-24 Thread Jose Narillos de Santos
Sorry all for these two easy questions: First, I have a matrix with trhee columns: A= Name Eight Puntuation Pepe 1,85 10 Paco 1,7 7 Pablo 1,82 6 I want to scatter the two columns (I could use pairs or other functions...) but the only doubt is that I cannot find the way to add in the as

[R] An unprofessional message

2009-12-24 Thread Maithili Shiva
Dear R helpers,   I understand that this is absolutely unprofessional on my part and this group doesn't entertain such things. I have been associted with this group since last 1 and half years and have been immensely benefited by the noble service rendred by many R helpers.   So I take this oppo

Re: [R] prcomp : plotting only explanatory axis arrows

2009-12-24 Thread Prof Brian Ripley
First, this is about biplot, not prcomp. Second, you seem to want to get a single-variable plot out of a biplot, which contradicts the 'bi' and hence I would not expect there to be a simple way to do this. The simplest thing to do would be to edit biplot.default via biplot.default <- stats::

Re: [R] What is ".Machine$double.eps"?

2009-12-24 Thread S Ellison
See ?.Machine >>> "Saji Ren" 24/12/2009 06:07:27 >>> Hello,everyone: I met this notation when I read the original code of function "quantile".There is one sentence as below: >eps <- 100 * .Machine$double.eps when I input ".Machine$double.eps" in R, it returns "[1] 2.220446e-16". Can anyone sh

[R] Error with Package "Kernlab" for SVM prediction

2009-12-24 Thread Vishal Thapar
Hi All, I am trying to use the "Kernlab" package for training and prediction using SVM's. I am getting the following error when I am trying to use the predict function: > predictSvm = predict(modelforSVM, testSeq); Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") : contrasts can be ap

Re: [R] Help with makeClusters for Snow

2009-12-24 Thread Ishwor Gurung
2009/12/24 SVM : > > Hi, > > Would it help if I was to install Ubuntu, and try loading it there? I've > seen quite a few comments about running rmpi on Linux, through google > search. > > Do you know if I have to install any specific packages/libraries before > trying to install RMpi? Yes. You nee