Re: [R] Help with makeClusters for Snow

2009/12/24 SVM sza...@myway.com: 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.

[R] Error with Package Kernlab for SVM prediction

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 applied

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

See ?.Machine Saji Ren saji@gmail.com 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

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

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 -

[R] An unprofessional message

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

[R] Two Easy questions

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

Re: [R] Two Easy questions

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

Re: [R] Two Easy questions

?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 narillosdesan...@gmail.com wrote: From:

Re: [R] as.Date question

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 marek.ja...@gmail.com To: r-help@r-project.org Sent: Thursday, December 24, 2009

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

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

[R] aggregate binary response data

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.

Re: [R] aggregate binary response data

?aggregate ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask grahamle...@btinternet.comwrote: 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

Re: [R] aggregate binary response data

may be ?table ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask grahamle...@btinternet.comwrote: 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

Re: [R] as.Date question

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 diver...@univecom.ch wrote: Mark, not sure

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

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] how to do multiple responses in a linear regression

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

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

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

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

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] as.Date question

zzz1 is POSIXct so looking at: as.Date.POSIXct function (x, ...) { z - floor(unclass(x)/86400) attr(z, tzone) - NULL structure(z, class = Date) } environment: namespace:base we see as.Date.POSIXct takes the POSIXct object, zzz1, and converts it to Date relative to GMT. There is no

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

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

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

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 Dec 24, 2009,

Re: [R] An unprofessional message

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

[R] Multiple CHOLMOD errors when attempting poisson glmm

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

[R] help in merging

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 either x or y.

Re: [R] help in merging

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, utkarsh.sing...@global-analytics.comwrote: Hi All,

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

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] How to separate a data set by its factors

Hi, On Thu, Dec 24, 2009 at 3:24 PM, James Rome jamesr...@gmail.com 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

[R] lapack error on Redhat Linux Fedora 11

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

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

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

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

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

Re: [R] An unprofessional message

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 xieyi...@gmail.com Phone: 515-294-6609 Web: http://yihui.name Department of Statistics, Iowa State University 3211 Snedecor Hall, Ames, IA On

Re: [R] An unprofessional message

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

[R] Question to use R plot GO pie chart

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:

Re: [R] An unprofessional message

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

Re: [R] An unprofessional message

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] How to separate a data set by its factors

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

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

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:

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

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

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

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 rusers...@gmail.com wrote: Hi,  I got two errors after i checking my package. The errors are as follows, Error1: * checking for unstated dependencies in R code ...

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

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 ggrothendi...@gmail.com The first one is a warning and

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

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 rusers...@gmail.com wrote:  Thanks.  BTW, if i have three functions, a,b and c. I only want the

[R] Help with SVM package Kernlab

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

[R] problem with .Rhistory

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

[R] problem with .Rhistory: please ignore

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

Re: [R] Help with SVM package Kernlab

Hi, Comments in line: On Thu, Dec 24, 2009 at 11:42 PM, Vishal Thapar vishaltha...@gmail.com 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

Re: [R] Help with SVM package Kernlab

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

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

Re: [R] Help with SVM package Kernlab

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 mpinputvector 2. is(train500$class) -- NULL OptionalFunction output 3. is(train500[1,5]) -- factor

Re: [R] Help with SVM package Kernlab

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.

Re: [R] Help with SVM package Kernlab

Hi, On Fri, Dec 25, 2009 at 12:49 AM, Vishal Thapar vishaltha...@gmail.com 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.