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

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

Re: [R] too large dimension problem

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 man

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

Dear Jorge, Thanks. But. (My R console width is set by default at 91, which is fine. Output for R>1: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

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

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 during

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

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 >

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

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

[R] Hmisc filled bands colors

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] too large dimension problem

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

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

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] Non-linear Weibull model for aggregated parasite data

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 be

[R] xtabs - missing combination

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 +

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

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

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

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

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

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 study

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

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

Re: [R] matched pair proportion test

Yes, thanks. It's exactly what I want. Annie On Sat, Dec 12, 2009 at 12:46 PM, Peter Dalgaard wrote: > 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.te

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

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 : > Dear all, > > I'm stuck in a seemingly trivial task that I need to perform f

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

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 a[match

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

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

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

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 wrote: > > I am trying to graphically represent a large set of data who's result is > not > strictly un

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

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] Create sequence given start and end vector

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?

[R] Create sequence given start and end vector

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.

[R] Call STATA code from within R

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: http://n4.nabble.com/Call-STATA-code-from-within-R-tp961885p9

[R] simple ts.plot question

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

Re: [R] by function ??

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

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

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] Regularized gamma function/ incomplete gamma function

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 ; 0 (which is eventually the ratio of the Incomplete gamma RV> function by the gamma function), and which is exactly what

Re: [R] by function ??

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

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

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,

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

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 hav

Re: [R] matched pair proportion test

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] simple ts.plot question

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

Re: [R] by function ??

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

*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. Th

[R] simple ts.plot question

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

[R] matched pair proportion test

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 https://st

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

David Winsemius 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 Wonderful, tha

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

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

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

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 ne

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

On 12/12/09 11:02, David Winsemius wrote: > > 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 cor

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

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

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

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 12.12.2009 00:28 >>> Hi R Exper

Re: [R] cube root on array

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] cube root on array

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 https://stat.ethz.ch/m

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

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 idea

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

> 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

Re: [R] About zero-inflation poisson model

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

Re: [R] SVM regression

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