Re: [R] SVM regression

2009-12-12 Thread Eleni Christodoulou
Thank you very much! Eleni On Fri, Dec 11, 2009 at 7:19 PM, Steve Lianoglou mailinglist.honey...@gmail.com wrote: Hi Eleni, On Dec 11, 2009, at 12:04 PM, Eleni Christodoulou wrote: Dear R users, I am trying to apply SVM regression for a set of microarray data. I am using the

Re: [R] About zero-inflation poisson model

2009-12-12 Thread Achim Zeileis
On Fri, 11 Dec 2009, Xiongqing Zhang wrote: Hello all, I am Xiongqing Zhang, come from Beijing of China. I know you from the web site: http://finzi.psych.upenn.edu/Rhelp08/2008-February/154627.html. I am not very clear about the R-project software. Look at http://www.R-project.org/

Re: [R] some problems with ram usage and warnings

2009-12-12 Thread Tom Knockinger
Datum: Fri, 11 Dec 2009 21:52:30 -0500 On Dec 11, 2009, at 11:08 AM, Tom Knockinger wrote: Hi, i am new to the R-project but until now i have found solutions for every problem in toturials, R Wikis and this mailing list, but now i have some problems which I can't solve with this

[R] read.csv to read output of system()?

2009-12-12 Thread Marianne Promberger
Dear list, I have a file that is comma delimited but contains some erroneous non-delimiter commas. I would like to replace these commas with semicolons and then read the correct file into R as a data frame. I want to do this from within R, without changing the original data file. My current

[R] cube root on array

2009-12-12 Thread Rodrigo Tsai
Dear R developers, is that right? -27^(1/3) [1] -3 c(-27,27)^(1/3) [1] NaN 3 i'm using sign( c(-27,27) ) * abs( c(-27,27)) ^(1/3) , thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] cube root on array

2009-12-12 Thread Marc Schwartz
On Dec 12, 2009, at 9:31 AM, Rodrigo Tsai wrote: Dear R developers, is that right? -27^(1/3) [1] -3 library(fortunes) fortune(^) Thomas Lumley: The precedence of ^ is higher than that of unary minus. It may be surprising, [...] Hervé Pagès: No, it's not surprising. At least to me...

[R] Antw: Export R output to Word/RTF?

2009-12-12 Thread Frank Bloos
I am using SWord from statconn: http://rcom.univie.ac.at/download.html It allows you to put r-commands into Word (similar to odfWeave). Output is directed to Word including figures and tables. It is still a beta-version but works fine. Frank Bloos Wenjie Lee wenjieleemaill...@gmail.com

Re: [R] read.csv to read output of system()?

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 7:54 AM, Marianne Promberger wrote: Dear list, I have a file that is comma delimited but contains some erroneous non-delimiter commas. I would like to replace these commas with semicolons and then read the correct file into R as a data frame. I want to do this from within

Re: [R] read.csv to read output of system()?

2009-12-12 Thread Marianne Promberger
Thanks for both replies. Let me start by giving a better minimal example, although indeed the regex replacement is not my problem. system(echo \var1,var2,var3\none,two,three\none,this is a comment,with commas.,three\ file.csv) On 12/12/09 11:02, David Winsemius wrote: You need to

Re: [R] read.csv to read output of system()?

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 12:01 PM, Marianne Promberger wrote: Thanks for both replies. Let me start by giving a better minimal example, although indeed the regex replacement is not my problem. system(echo \var1,var2,var3\none,two,three\none,this is a comment,with commas.,three\ file.csv) On

Re: [R] read.csv to read output of system()?

2009-12-12 Thread Marianne Promberger
David Winsemius dwinsem...@comcast.net 12-Dec-09 17:12: txt - gsub(( [A-Za-z]+),,\\1;,readLines(file.csv)) read.csv(textConnection(txt), header=TRUE) var1 var2 var3 1 onetwo three 2 one this is a comment;with commas. three

[R] matched pair proportion test

2009-12-12 Thread annie Zhang
Hi, ALL, Is there any function in R that does the exact test for the matched pair proportions (one sided), which I assume is binomial(b+c, .5). Thanks, Annie [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] simple ts.plot question

2009-12-12 Thread arka sinha
* I have a simple question regarding plots of time series in R. I have to plot conc against time **for each individual and display in the same panel for the in-built dataset Indometh in R. * ***I have six time series, say subject1.ts, subject2.ts, ., subject6.ts. The observations are

[R] simple ts.plot question

2009-12-12 Thread arka sinha
*Respected Sir, I have a simple question regarding plots of time series in R. I have to plot conc against time **for each individual and display in the same panel for the in-built dataset Indometh in R. * *I have six time series, say subject1.ts, subject2.ts, ., subject6.ts. The

Re: [R] by function ??

2009-12-12 Thread L.A.
Thanks for all the help, They all worked, But I'm stuck again. I've tried searching, but I not sure how to word my search as nothing came up. Here is my new hurdle, my data has 7 abservations and my results have 2 answers: Here is my data LEAID ratio 3 6307 0.720 1 6307

[R] simple ts.plot question

2009-12-12 Thread arka
* I have a simple question regarding plots of time series in R. I have to plot conc against time **for each individual and display in the same panel for the in-built dataset Indometh in R. * *I have six time series, say subject1.ts, subject2.ts, ., subject6.ts. The observations are

Re: [R] matched pair proportion test

2009-12-12 Thread Peter Dalgaard
annie Zhang wrote: Hi, ALL, Is there any function in R that does the exact test for the matched pair proportions (one sided), which I assume is binomial(b+c, .5). binom.test should fit nicely, I think. Thanks, Annie [[alternative HTML version deleted]]

[R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-12 Thread Hamid
Dear all, here my problem my be somone can help to solve this. I have tow timeseries from different stock market with different length (diff ca. 4000 ) Now I would to add the missing times of the one of this series with proper time (they are minute data) and set the value to 0 since I need to

[R] rpart - classification and regression trees (CART)

2009-12-12 Thread Katie N
Hi, I had a question regarding the rpart command in R. I used seven continuous predictor variables in the model and the variable called TB122 was chosen for the first split. But in looking at the output, there are 4 variables that improve the predicted membership equally (TB122, TB139, TB144,

Re: [R] by function ??

2009-12-12 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of L.A. Sent: Saturday, December 12, 2009 12:39 PM To: r-help@r-project.org Subject: Re: [R] by function ?? Thanks for all the help, They all worked, But I'm stuck again.

Re: [R] Regularized gamma function/ incomplete gamma function

2009-12-12 Thread Martin Maechler
RV I would be very grateful if you could help me with: RV Given the regularized gamma function Reg=int_0^r RV (x^(k-1)e^(-x))dx/int_0^Inf (x^(k-1)e^(-x))dx ; 0rInf RV (which is eventually the ratio of the Incomplete gamma RV function by the gamma function), and which is

Re: [R] save an object by dynamicly created name

2009-12-12 Thread Henrik Bengtsson
saveObject() and loadObject() are using save() and load() in base. The default is to compress the data when saving, which take some time. Using saveObject(..., compress=FALSE) is probably faster, but takes up more disk space. ...and make sure you don't work towards a file system over a network,

Re: [R] by function ??

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 3:38 PM, L.A. wrote: Thanks for all the help, They all worked, But I'm stuck again. I've tried searching, but I not sure how to word my search as nothing came up. Here is my new hurdle, my data has 7 abservations and my results have 2 answers: Here is my data

[R] simple ts.plot question

2009-12-12 Thread arka sinha
*Respected Sir, I have a simple question regarding plots of time series in R. I have to plot conc against time **for each individual and display in the same panel for the in-built dataset Indometh in R. * ***I have six time series, say subject1.ts, subject2.ts, ., subject6.ts. The

[R] Call STATA code from within R

2009-12-12 Thread tonitogomez
Dear R Users, Do you know whether is there any way of calling STATA from within R (possibly in a similar way that is possible to call WinBUGS from within R using the function Bugs)?? Thanks, Manuel -- View this message in context:

[R] Create sequence given start and end vector

2009-12-12 Thread Kevin Ummel
How can I create the following without the 'for' loop? start=c(1,10,20) end=c(4,15,27) out=c() for (i in 1:length(start)) { out=c(out,start[i]:end[i]) } out [1] 1 2 3 4 10 11 12 13 14 15 20 21 22 23 24 25 26 27 I know there must be an easier (and, hopefully, faster) way.

Re: [R] Create sequence given start and end vector

2009-12-12 Thread Jorge Ivan Velez
Hi Kevin, Here is a suggestion using mapply(): start - c(1,10,20) end - c(4,15,27) do.call(c, mapply( seq, start, end)) See ?mapply and ?do.call for more information. HTH, Jorge On Sat, Dec 12, 2009 at 2:27 PM, Kevin Ummel wrote: How can I create the following without the 'for' loop?

Re: [R] Create sequence given start and end vector

2009-12-12 Thread Joe King
Also what about c(seq(1,4,1),seq(10,15,1),seq(20,27,1)) Joe King j...@joepking.com Never throughout history has a man who lived a life of ease left a name worth remembering. --Theodore Roosevelt -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-12 Thread milton ruser
Hi Philip, I must confess that I not understood what is the problem. Could you clarify it a little bit more? Cheers miltinho brazil=toronto On Sat, Dec 12, 2009 at 12:39 AM, philip robinson robin...@students.wwu.edu wrote: I am trying to graphically represent a large set of data who's

[R] multi-dimensional array with different number of dimensions?

2009-12-12 Thread parkbomee
Hi, Is it possible to assign to an array with different dimensions? That is to say, supposing a three dimensional array, the third dimension of the array has matrices of different sizes? array , , 1 [1] [2] [3] [1] 111 , , 2 [1] [2] [3] [1] 111 [2] 11

[R] Replace NAs in a range of data frame columns

2009-12-12 Thread Michael Scharkow
Dear all, I'm stuck in a seemingly trivial task that I need to perform for many datasets. Basically, I want to replace NA with 0 in a specified range of columns in a dataframe. I know the first and last column to be recoded only by its name. I can select the columns starting like this

Re: [R] Replace NAs in a range of data frame columns

2009-12-12 Thread Patrizio Frederic
hi Michael, the following code should work b - a[match('first',names(a)): match('last',names(a))] b[is.na(b)]-0 a[match('first',names(a)): match('last',names(a))] - b cheers, Patrizio 2009/12/13 Michael Scharkow mich...@underused.org: Dear all, I'm stuck in a seemingly trivial task that I

Re: [R] matched pair proportion test

2009-12-12 Thread annie Zhang
Yes, thanks. It's exactly what I want. Annie On Sat, Dec 12, 2009 at 12:46 PM, Peter Dalgaard p.dalga...@biostat.ku.dkwrote: annie Zhang wrote: Hi, ALL, Is there any function in R that does the exact test for the matched pair proportions (one sided), which I assume is binomial(b+c, .5).

Re: [R] multi-dimensional array with different number of dimensions?

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 5:36 PM, parkbomee wrote: Hi, Is it possible to assign to an array with different dimensions? That is to say, supposing a three dimensional array, the third dimension of the array has matrices of different sizes? Use a list. Or populate with NA's array , , 1

Re: [R] Create sequence given start and end vector

2009-12-12 Thread Charles C. Berry
On Sat, 12 Dec 2009, Jorge Ivan Velez wrote: Hi Kevin, Here is a suggestion using mapply(): start - c(1,10,20) end - c(4,15,27) do.call(c, mapply( seq, start, end)) ...which is what I would usually do. But for heavy duty applications, the IRanges package and function may be worth

Re: [R] Replace NAs in a range of data frame columns

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 6:15 PM, Patrizio Frederic wrote: hi Michael, the following code should work b - a[match('first',names(a)): match('last',names(a))] b[is.na(b)]-0 This might not throw an error: b - apply(a[match('first',names(a)): match('last',names(a))], 1:2,

Re: [R] Have you used RGoogleDocs and RGoogleData?

2009-12-12 Thread Duncan Temple Lang
Farrel Buchinsky wrote: It Works! Thanks a lot! Its great. Thanks for letting me know. Glad that fixed things for you. What were your few minor, but important, changes - in a nutshell. I will not understand unless you describe it as high level issues. Basically, recognizing the type of a

[R] xtabs - missing combination

2009-12-12 Thread Patrick Hausmann
Dear list, I am trying to make a contingency table with xtabs but I am getting a 0 where I expect a 'NA'. Here is a simple example: options(stringsAsFactors = FALSE) rn - LETTERS[1:4] df1 - data.frame(r07 = rep(rn, each=4), r08 = rep(rn, 4), value = 1:16) xtabs(value ~ r07 +

[R] Non-linear Weibull model for aggregated parasite data

2009-12-12 Thread Dan E.
Hi, I am trying to fit a non-linear model for a parasite dataset. Initially, I tried log-transforming the data and conducting a 2-way ANCOVA, and found that the equal variance of populations and normality assumptions were violated. Gaba et al. (2005) suggests that the Weibull Distribution is

Re: [R] confint for glm (general linear model)

2009-12-12 Thread casperyc
for an example, counts - c(18,17,15,20,10,20,25,13,12) outcome - gl(3,1,9); treatment - gl(3,3) glm.D93 - glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asymptotic normality to verify the confidence interval (confint.default(glm.D93))

[R] too large dimension problem

2009-12-12 Thread Moohwan Kim
Dear R family When I run the command below, the error message came up. It seems like the problem is about computer capacity. It would be appreciated if anyone could give me a solution. ### N - 415884 tau - diag(1, N)[c(N, 1:(N - 1)),] Error in array(0, c(n, p)) : 'dim' specifies too

[R] Hmisc filled bands colors

2009-12-12 Thread Joe King
Hi all, i am still using the Hmisc package and I like the filled bands part, is there a way though to have different groups have different color of bands, maybe a lighter version of the color of the line that is used? --- Joe King, M.A. Ph.D. Student

[R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-12 Thread Viju Moses
I'm facing this problem on R GUI version 2.10.0 on Windows Vista. I have not changed Windows settings or R GUI settings much except to change from MDI to SDI. Someone else reported this problem a few months ago: https://stat.ethz.ch/pipermail/r-help/2009-April/195714.html but it wasn't followed

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-12 Thread Jorge Ivan Velez
Hi Viju, Here is a suggestion: R options(width = 80) R 1:120 HTH, Jorge On Sat, Dec 12, 2009 at 10:34 PM, Viju Moses wrote: I'm facing this problem on R GUI version 2.10.0 on Windows Vista. I have not changed Windows settings or R GUI settings much except to change from MDI to SDI.

[R] Easily switch between R GUI and R term on Windows?

2009-12-12 Thread Viju Moses
Is there an easy way to switch between R GUI and R terminal? I'm currently using R GUI 2.10.0 on Windows Vista and would like to see something like the R terminal found in Linux. But I may also want to switch back to GUI. I don't remember seeing the terminal option at the R download site, or

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-12 Thread Viju Moses
Dear Jorge, Thanks. But. (My R console width is set by default at 91, which is fine. Output for R1:120 wraps correctly already. ) The problem is that when I am typing a command, it does not wrap like it does in this email while typing. As I am reach the right border, a horizontal scroll

Re: [R] too large dimension problem

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 8:31 PM, Moohwan Kim wrote: Dear R family When I run the command below, the error message came up. It seems like the problem is about computer capacity. It would be appreciated if anyone could give me a solution. You would want to just store such an object but also to

Re: [R] confint for glm (general linear model)

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 8:19 PM, casperyc wrote: for an example, counts - c(18,17,15,20,10,20,25,13,12) outcome - gl(3,1,9); treatment - gl(3,3) glm.D93 - glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asymptotic normality to verify