[R] Keeping only one date for each patient ID and disease SITE

Dear List Subscribers, I am working on the following problem and was wondering if there is some command or set of commands to solve it: Thank you in advance, Eric 1. The dataset Cancer0 may have multiple dates of treatment (DATE) for each patient (ID) with a given disease (SITE). C

Re: [R] Summing data based on certain conditions

This should do the trick: data$date <- as.Date(data$date, '%m/%d/%Y') data$month <- format(data$date, '%Y-%m') by(data$rammday, data$month, sum) Hope that helps, Jason -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Stephan Kol

Re: [R] memory error

Janet Choate-2 wrote: > > Thanx for clarification on stating my problem, Charlie. > > I am attempting to merge to files, i.e.: > hi39 = merge(comb[,c("hillID","geo")], hi.h39, by=c("hillID")) > > if this is relevant or helps to explain: > the file 'comb' is 3 columns and 1127 rows > the file '

Re: [R] Sharing levels across multiple factor vectors

Sorry for spamming. I swear I had worked on that problem a long time before posting. But I just figured it out: I have to change the values, which are represented as integers, not strings. So the following code will do it: df <- data.frame ( a = factor( c( "bob", "alice", "bob" ) ),

[R] Sharing levels across multiple factor vectors

Hi, I've got a data frame with multiple factor columns, but they should share the same set of labels, such as this tiny example: df <- data.frame ( a = factor( c( "bob", "alice", "bob" ) ), b = factor( c( "kenny", "alice", "alice" ) ) ); In my data, though, the strings are enorm

Re: [R] trying to understand lme() results

Hi: On Wed, Mar 31, 2010 at 2:31 PM, array chip wrote: > Hi, I have very simple balanced randomized block design where I total have > 48 observations of a measure of weights of a product, the product was > manufactured at 4 sites, so each site has 12 observations. I want to use > lme() from nlm

Re: [R] memory error

Thanx for clarification on stating my problem, Charlie. I am attempting to merge to files, i.e.: hi39 = merge(comb[,c("hillID","geo")], hi.h39, by=c("hillID")) if this is relevant or helps to explain: the file 'comb' is 3 columns and 1127 rows the file 'hi.h39' is 5 columns and 19758564 rows i s

[R] model set up question

I need to compare gene expression differences between multiple line pairs of alcohol preferring and non-preferring rat lines. I have 5 such line pairs, 3 are unrelated but two were derived independently from the same parent stock. For each line, there are 10 samples. I'll be testing multiple genes,

Re: [R] pdf files in loops

On Wed, 31 Mar 2010 22:06:48 -0400 James Rome wrote: > print() did not help, and I get strange messages about mode(onefile) > not being changed: Either: R> pic <- histogram(...) R> print(pic) or R> print(histogram(...)) > I would actually like all the histograms in one pdf file too... Then

Re: [R] pdf files in loops

On 3/31/2010 10:01 PM, Berwin A Turlach wrote: G'day James, On Wed, 31 Mar 2010 21:44:31 -0400 James Rome wrote: > I need to make a bunch of PDF files of histograms. > [...] > What am I doing wrong? > http://cran.ms.unimelb.edu.au/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-

Re: [R] pdf files in loops

On 3/31/2010 10:01 PM, Berwin A Turlach wrote: G'day James, On Wed, 31 Mar 2010 21:44:31 -0400 James Rome wrote: > I need to make a bunch of PDF files of histograms. > [...] > What am I doing wrong? > http://cran.ms.unimelb.edu.au/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-

Re: [R] pdf files in loops

G'day James, On Wed, 31 Mar 2010 21:44:31 -0400 James Rome wrote: > I need to make a bunch of PDF files of histograms. [...] > What am I doing wrong? http://cran.ms.unimelb.edu.au/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f HTH. Cheers, Berwin =

[R] pdf files in loops

I need to make a bunch of PDF files of histograms. I tried gatelist = unique(mdf$ArrivalGate) for( gate in gatelist) { outfile = paste("../", airport, "/", airport, "taxiHistogram", gate, ".pdf", sep="") pdf(file = outfile, width = 10, height=8, par(lwd=1)) title=paste("Taxi tim

[R] predicted time length differs from survfit.coxph:

Hello All, Does anyone know why length(fit1$time) < length(fit2$n) in survfit.coxph output? Why is the predicted time length is not the same as the number of samples (n)? I tried: example(survfit.coxph). Thanks, parmee > fit2$n [1] 241 > fit2$time [1] 031326061 152 15

Re: [R] memory error

Janet Choate-2 wrote: > > Hi R community, > i have what appears to be a memory allocation problem: > > R(51150) malloc: *** mmap(size=158068736) failed (error code=12) > *** error: can't allocate region > *** set a breakpoint in malloc_error_break to debug > > can anyone tell me how to increas

Re: [R] plotting a function

Take a look at the examples in ?curve. As suggested by Erik, in the future please read the posting guide so you can get a more accurate response. Regards, Francisco Dr. Francisco J. Zagmutt Vose Consulting 1643 Spruce St., Boulder Boulder, CO, 80302 USA www.voseconsulting.com Jin wrote: He

Re: [R] Code is too slow: mean-centering variables in a data framebysubgroup

Dimitri, You might try applying ave() to each column. E.g., use f2 <- function(frame) { for(i in 2:ncol(frame)) { frame[,i] <- ave(frame[,i], frame[,1], FUN=function(x)x/mean(x,na.rm=TRUE)) } frame } Note that this returns a data.frame and retains the grouping column (the first)

Re: [R] Bar plots with bars made of stacked text

On 04/01/2010 05:09 AM, Thomas Levine wrote: I would like to make bar plots where the bars are composed of text like this: http://www.thomaslevine.com/lowres/text_bars.png Is there a package that will help me with this? Thanks Hi Tom, Suppose you have a data frame like this: zoodat names mo

Re: [R] lattice: how to add points to the plot generated by levelplot()?

You can use levelplot(runif(100)~rep(1:10, each=10)+rep(1:10, times=10)) trellis.focus("panel", 1, 1, highlight=FALSE) lpoints(runif(100,0,10), runif(100,0,10), pch=2, col=2, cex=2) trellis.unfocus() # or levelplot(runif(100)~rep(1:10, each=10)+rep(1:10, times=10), panel=function(...

Re: [R] ODD and EVEN numbers

Thank you guys for responding to my question. Have a great day!:-) -- View this message in context: http://n4.nabble.com/ODD-and-EVEN-numbers-tp1747032p1747411.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mail

[R] barplot and line

Hi r-users,   I have this data below and would like to plot a barplot overlap with a line. This is my data: > hist_50     pdf_obs pdf_gen.50 1  0.00   0.00 2  0.083156   0.125366 3  0.132196   0.158230 4  0.126866   0.149432 5  0.120469   0.127897 6  0.121535   0.104096 7  0.103412   0.0821

Re: [R] strip.custom with strip on left for three conditioning variables

Hi This is untested but I think for this to work you will need to add scales eg scales = list(x = list(relation = "free"), y = list(relation = "free")), also you will need to adjust the axis limits and a few other things. With the axis limits it is easier to start with a small number to get t

Re: [R] Problem comparing hazard ratios

P.S. Here is some code that is more directly relevant. Frank cphadjchisq <- function(adjust, candidates, S) { require(rms) f <- coxphFit(adjust, S, type='right', method='efron') if(f$fail) stop('could not fit a model with just adjust') ll2 <- -2*f$loglik[2] chisq.adjust <

Re: [R] Summing data based on certain conditions

?by may also be helpful. Stephan Steve Murray schrieb: Dear all, I have a dataset of 1073 rows, the first 15 which look as follows: data[1:15,] date year month day rammday thmmday 1 3/8/1988 1988 3 81.430.94 2 3/15/1988 1988 3 152.860.66 3 3/22/1988 1

[R] trying to understand lme() results

Hi, I have very simple balanced randomized block design where I total have 48 observations of a measure of weights of a product, the product was manufactured at 4 sites, so each site has 12 observations. I want to use lme() from nlme package to estimate the standard error of the product weight.

Re: [R] Install package of "BRugs"?

Thanks! I updated R to 2.10.1 and everything is fine now! Hua --- On Wed, 3/31/10, Uwe Ligges wrote: > From: Uwe Ligges > Subject: Re: [R] Install package of "BRugs"? > To: "David Winsemius" > Cc: "Hua Li" , "R-help" > Date: Wednesday, March 31, 2010, 2:53 AM > > > On 30.03.2010 22:01, Dav

Re: [R] Bar plots with bars made of stacked text

Here is a start: mydat <- list( Morning=c('Billy','Katy','Jonny'), Noon=c('Billy','Patty','Suzie','Tom','Vicky'), Evening=c('Amy','Jonny','Sally') ) sx <- rep( seq_along(mydat), sapply(mydat, length) ) sy <- unlist( lapply( mydat, seq_along ) ) plot(sx, sy, xlab='', ylab='', xax

Re: [R] plotting a function

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Jin wrote: Hello Dear, I am trying to plot a function to see a minimum point (actually, using "optim"). For example, 1. y=f(x) 2. x has a range 3. pl

[R] plotting a function

Hello Dear, I am trying to plot a function to see a minimum point (actually, using "optim"). For example, 1. y=f(x) 2. x has a range 3. plot(x,y) to see a point x minimizing y I tried "plot(x,y)", but it made an error. Am I doing a right way? Thanks, Jin -- View this message in context: htt

[R] memory error

Hi R community, i have what appears to be a memory allocation problem: R(51150) malloc: *** mmap(size=158068736) failed (error code=12) *** error: can't allocate region *** set a breakpoint in malloc_error_break to debug can anyone tell me how to increase the memory size for R on mac os X? thank

Re: [R] lattice: how to add points to the plot generated bylevelplot()?

Tip 1: This sort of thing is done with panel functions. ?xyplot and ?panel.xyplot is where to start. Tip 2: grid functions are generally not necessary, as translations of the common ones are part of lattice. See ?panel.points Tip 3: Deepayan's book provides a gentler explanation of these matters

Re: [R] ODD and EVEN numbers

(foo %% 2) == 0 See ?"%%" HTH Stephan tj schrieb: Hi, anyone here who knows how to determine if an integer is "odd" or "even" in R? Thanks. tj __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] strip.custom with strip on left for three conditioning variables

I want to use a strip.custom as with useOuterStrips for three conditioning variables. useOuterStrips restricts this to only two conditioning variables, and I cannot figure out how to write strip.custom properly to do this. library(lattice) mtcars$HP <- equal.count(mtcars$hp) #with two factors x2

Re: [R] Summing data based on certain conditions

?ave or ?tapply Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Steve Murray Sent: Wednesday, March 31, 2010 10:20 AM To: r-help@r-project.org Subject: [R] Summing data based on c

Re: [R] ODD and EVEN numbers

On Mar 31, 2010, at 12:43 PM, tj wrote: Hi, anyone here who knows how to determine if an integer is "odd" or "even" in R? Thanks. > 2 %% 2 == 0 [1] TRUE > 3 %% 2 == 0 [1] FALSE is.even <- function(x){ x %% 2 == 0 } > is.even(2) [1] TRUE -- David Winsemius, MD West Hartford, CT ___

Re: [R] lattice: how to add points to the plot generated by levelplot()?

Czerminski, Ryszard wrote: I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() It's best to include a small example of what you have done so far, and what you want to further do. Actual R code with a small sample dataset is gr

Re: [R] ODD and EVEN numbers

See %% under ?Arithmetic tj wrote: Hi, anyone here who knows how to determine if an integer is "odd" or "even" in R? Thanks. tj __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

On 03/31/2010 11:44 AM, Gabor Grothendieck wrote: > I got an error message *AND* R becomes unusable and had to be restarted. > > Using dependency walker it seems to be complaining about R.dll so I > copied all of my ...\R\bin\R*.dll files to .. > ...\win-library\2.10\Rgraphviz\libs. > > Then it c

[R] Summing data based on certain conditions

Dear all, I have a dataset of 1073 rows, the first 15 which look as follows: > data[1:15,]     date year month day rammday thmmday 1   3/8/1988 1988 3   8    1.43    0.94 2  3/15/1988 1988 3  15    2.86    0.66 3  3/22/1988 1988 3  22    5.06    3.43 4  3/29/1988 1988 3  29  

[R] How to work with the result of JRip?

Hello, I want to work with the result of JRip, but I don't know how to access it. I produce a rule like below: > rule = JRip(Cls ~., data=dd) The result is: JRIP rules: === (C7 <= 1) => Cls=TRUE (384.0/0.0) => Cls=FALSE (1344.0/0.0) \par Number of Rules : 2 I want to work with

[R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Hi Martin, It is really a 'PATH' problem. After adding C:\Program Files\Graphviz2.20\bin to the 'PATH' environment variable, the Rgraphviz package can be loaded without any error messages. Sorry that I ONLY set my 'PATH' environment variable correctly for R but not for Graphviz. Thank y

[R] change mbmdr odds ratio

I am using the package 'mbmdr' based on the model-base multifactor dimensionality reduction method by Calle et al. In the paper, mbmdr is used for a case-control study. My data is not case-control, i.e. I have a lot more controls than cases. Does anyone know if the odds ratio in the fuction mbm

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Hi Gabor, I just used your R code below and the code worked properly without any error messages. Attached is the output graph, which may be the one you expected. library("Rgraphviz") set.seed(123) V <- letters[1:10] M <- 1:4 g1 <- randomGraph(V, M, 0.2) g1 plot(g1) Thanks. Howard On W

[R] ODD and EVEN numbers

Hi, anyone here who knows how to determine if an integer is "odd" or "even" in R? Thanks. tj -- View this message in context: http://n4.nabble.com/ODD-and-EVEN-numbers-tp1747032p1747032.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Problem comparing hazard ratios

Ravi Varadhan wrote: Frank, Is there an article that discusses this idea of bootstrapping the ranks of the likelihood ratio chi-square Statistics to assess relative importance of predictors in time-to-event data (specifically Cox PH model)? Thanks, Ravi. Do require(rms); ?anova.rms and see

[R] lattice: how to add points to the plot generated by levelplot()?

I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() similarly as in basic R graphic sequence: plot(...); points(...), but it does not work. I was reading documentation for lattice, but so far without much success. I would very much a

[R] AHRQ Patient Quality Indicators

Is anyone aware of R code that mimic's AHRQ's SAS code for their Prevention Quality Indicators (PQI)? Don't see it anywhere, but wanted to see if anyone else knew of anything. Many thanks -- View this message in context: http://n4.nabble.com/AHRQ-Patient-Quality-Indicators-tp1747243p1747243.

Re: [R] reshaping data

On Mar 31, 2010, at 2:37 PM, Kim Jung Hwa wrote: Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like to convert it back to original format (as a data frame)... is there an equivalent easy command

Re: [R] reshaping data

Try: as.data.frame(with(x, tapply(Val, list(Var1, Var2), sum))) On Wed, Mar 31, 2010 at 3:37 PM, Kim Jung Hwa wrote: > Thanks Henrique and Stephan for your reply! > > Henrique, I'm planning to do some arthitmetic operations on tranformed > (matrix) data and then would like to convert it back t

Re: [R] Off Topic: teenie weenie numbers -- Was: Precision level

There is a story (apocryphal?) about Fred Hoyle, many years ago, having come to the close of a public lecture about his work in Cosmology (I seem to recall that it was to the British Astronomical Association, a society of amateur astronomers, and therefore knowledgeable). He was taking questions.

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

I got an error message *AND* R becomes unusable and had to be restarted. Using dependency walker it seems to be complaining about R.dll so I copied all of my ...\R\bin\R*.dll files to .. ...\win-library\2.10\Rgraphviz\libs. Then it complained about IESHIMS.dll so I copied \Program Files\Internet

Re: [R] reshaping data

Actually, apart from melt() in reshape package. On Wed, Mar 31, 2010 at 2:37 PM, Kim Jung Hwa wrote: > Thanks Henrique and Stephan for your reply! > > Henrique, I'm planning to do some arthitmetic operations on tranformed > (matrix) data and then would like to convert it back to original format (

Re: [R] reshaping data

Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like to convert it back to original format (as a data frame)... is there an equivalent easy command for this too? Thanks, On Wed, Mar 31, 2010 at 2:26 PM

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Peter, it works!!! Thanks a lot to all of you! Yes I mistook a data.frame for a matrix Now both David's and Dennis' solutions work! I am really a beginner with R, and I think it's a wonderful tool, but sometimes a little difficult to use only learning from the manuals and help files... I ho

Re: [R] reshaping data

Hi Kim, look at the reshape() command with direction="wide". Or at the reshape package. HTH, Stephan Kim Jung Hwa schrieb: Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a matrix),

Re: [R] reshaping data

Try this: xtabs(Val ~ Var1 + Var2, data = x) On Wed, Mar 31, 2010 at 3:23 PM, Kim Jung Hwa wrote: > Hi All, > > Can someone help me reshape following data: > > Var1 Var2 Val > A X 1 > A Y 2 > A Z 3 > B X 4 > B Y 5 > B Z 6 > > to some kind of matrix/tabular format (preferably as a matrix), may b

[R] reshaping data

Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a matrix), may be like Var1 X Y Z A 1 2 3 B 4 5 6 Any help would be greatly appreciated, Kim [[alternative HTML version deleted]]

Re: [R] interpretation of p values for highly correlated logistic analysis

Hi Claus, welcome to the wonderful world of collinearity (or multicollinearity, as some call it)! You have a near linear relationship between some of your predictors, which can (and in your case does) lead to extreme parameter estimates, which in some cases almost cancel out (a coefficient of

[R] Bar plots with bars made of stacked text

I would like to make bar plots where the bars are composed of text like this: http://www.thomaslevine.com/lowres/text_bars.png Is there a package that will help me with this? Thanks Tom __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

Re: [R] library sets: A & EMPTY does not work; gset_intersection(A,EMPTY) works

Ryszard, You've made me take a closer look and now I do think that you've found a bug. After a quick look at the package vignette, I see that the authors have indeed overloaded "&" and so it should work for your example. The problem seems to be the order of the class attribute which is used to c

[R] Comparing times in two time zones

I have data that is collected in two different time zones z$OnDateTime <- as.POSIXct(runway$OnDateTime, tz = "EST5EDT", format="%m/%d/%Y %H:%M") is in Eastern time, and zi$ActualOnLocal <- as.POSIXct(oooi$ActualOnLocal, tz="MST7MDT", format="%m/%d/%Y %H:%M") is in Mountain". I converted the runway

Re: [R] Off Topic: teenie weenie numbers -- Was: Precision level

I spoke with a theoretical physicist and he said he encountered 10^-120. Has something to do with attempts to describe/explain the universe... Dimitri On Fri, Mar 26, 2010 at 2:33 PM, Steve Lianoglou wrote: > On Fri, Mar 26, 2010 at 1:11 PM, Barry Rowlingson > wrote: >> On Fri, Mar 26, 2010 at 4

Re: [R] dot plot

Hi kayj kayj schrieb: Also, I do nto know how to combine two plots into one plot? Have a look at ?par . There is some argument "add" or "new" (i am not sure) that you have to set to TRUE after making the first plot and then you can add a second plot with running plot() once more. thanks a

Re: [R] GLM / large dataset question

Geelman, This appears to be your first post to this list. Welcome to R. Nearly 2 days is quite a long time to wait though, so you are unlikely to get a reply now. Feedback : the question seems quite vague and imprecise. It depends on which R you mean (32bit/64bit) and how much ram you have. It

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

On 03/31/2010 08:18 AM, Gabor Grothendieck wrote: > By the way, just in case you did not read the entire message R crashed I think you mean that you got an error message, not that R became unusable? More below... > when I tried to run the code from the vignette. > > On Wed, Mar 31, 2010 at 10:52

Re: [R] Error "singular gradient matrix at initial parameter estimates" in nls

Corrado: Over parameterization/non-identifiability is not determined by the ratio of the number of data values to the number of parameters: if you try to fit a scatter of a zillion points that lie near a straight line to a model with curvature -- the 4p logistic function, say -- you're over-parame

Re: [R] library sets: A & EMPTY does not work; gset_intersection(A,EMPTY) works

It seems that "A & B" works the same way as "gset_intersection(A,B)" as long as A and B are not empty... see below: fuzzy_logic("Yager", p=2) support <- universe <- c('a','b','c','d') A <- gset(support=support, memberships=c(0.2, 0.2, 0.9, 0), universe=universe) B <- gset(support=support, membersh

Re: [R] position of mismatches in all.equal()?

I've actually found a solution to my problem. In case someone could be interested, the which() function seems to do the trick: > test2$let2 <- factor(test2$let2, levels=letters[1:10]) > which(test1[,3]!=test2[,3]) [1] 7 > which(test1[,1:3]!=test2[,1:3]) [1] 27 Le 3/31/2010 15:13, Ivan Caland

[R] Boundary kernel for kernel density plot?

Hello, Is a kernel called the "boundary" kernel available for kernel density plots? In the help, I can only see the following kernels: "gaussian", "rectangular", "triangular", "epanechnikov", "biweight", "cosine" or "optcosine". The boundary kernel apparently improves the estimate in a neighbo

Re: [R] library sets: A & EMPTY does not work; gset_intersection(A,EMPTY) works

On 2010-03-31 9:30, Peter Ehlers wrote: Unless I'm missing something, I don't see any method in pkg:sets for intersection other than gset_intersection. Whoops, a bit quick on the draw. There are of course also set_intersection and cset_intersection, but not AFAICS any method for `&`. -Peter E

Re: [R] Simplifying particular piece of code

Thanks for help, Gustaf! (Kan man säga "man tackar" oxa?) :-) On Wed, Mar 31, 2010 at 17:34, Gustaf Rydevik wrote: > On Wed, Mar 31, 2010 at 5:11 PM, Sergey Goriatchev wrote: >> but >> >> data <- merge(data,data.list) >> >> works. >> >> Neither data or data.list is a list, so do.call does not w

[R] predict.Arima: warnings from xreg magic

When I run predict.Arima in my code, I get warnings like: Warning message: In cbind(intercept = rep(1, n), xreg) : number of rows of result is not a multiple of vector length (arg 1) I think this is because I'm not running predict.Arima in the same environment that I did the fit, so the d

Re: [R] Simplifying particular piece of code

On Wed, Mar 31, 2010 at 5:11 PM, Sergey Goriatchev wrote: > but > > data <- merge(data,data.list) > > works. > > Neither data or data.list is a list, so do.call does not work. > I am very weak on lists, never used them before > > Best, > Sergey Hej Sergey, Ok; I was wondering if the apply thing

Re: [R] library sets: A & EMPTY does not work; gset_intersection(A,EMPTY) works

Unless I'm missing something, I don't see any method in pkg:sets for intersection other than gset_intersection. So you're using the base R function `&` whose help page tells you that its arguments should be vectors. Yours aren't. -Peter Ehlers On 2010-03-31 8:50, Czerminski, Ryszard wrote: Whe

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

By the way, just in case you did not read the entire message R crashed when I tried to run the code from the vignette. On Wed, Mar 31, 2010 at 10:52 AM, Gabor Grothendieck wrote: > Based on your success I thought I would try again. I am not sure why I > had more success this time but this time I

Re: [R] mcmcglmm starting value example

Apparently not, since this your 3rd unanswered thread to r-help this month about this package. Please read the posting guide and find out where you should send questions about packages. Then you might get an answer. "ping chen" wrote in message news:975148.47160...@web15304.mail.cnb.yahoo.c

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Driss, David is right - you should include data/code. Nevertheless, the error you get suggests that what you call a "matrix" is in fact a data.frame; that's usually a crucial distinction. There are two things you might try: Dennis' way: lm(as.matrix(a) ~ b) David's way: for(i in seq_along(

Re: [R] GEE for a timeseries of count (one cluster)

Contact the authors of those packages ? "miriza" wrote in message news:1269981675252-1745896.p...@n4.nabble.com... > > Hi! > > I was wondering if there were any packages that would allow me to fit a > GEE > to a single timeseries of counts so that I could account for > autocorrelation > in the

Re: [R] zero standard errors with geeglm in geepack

You may not have got an answer because you posted to the wrong place. Its a question about a package. Please read the posting guide. "miriza" wrote in message news:1269886286228-1695430.p...@n4.nabble.com... > > Hi! > > I am using geeglm to fit a Poisson model to a timeseries of count data as >

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Based on your success I thought I would try again. I am not sure why I had more success this time but this time I got this far even though I did not change my path at all or make any system changes from what I had before. I did reinstall Rgraphviz but used the previously installed graphviz. I tr

[R] tm package- remove stowords failling

Hi, I just noticed that by inspecting the matrix term that no all stopwords are removed, does someone know how to fix that? library(tm) data("crude") d<-tm_map(crude, removeWords, stopwords(language='english')) dt<-DocumentTermMatrix(d,control=list(minWordLength=3, minDocFreq=2)) inspect( dt) I

Re: [R] Question about 'logit' and 'mlogit' in Zelig

Abraham, This appears to be your 3rd unanswered post to r-help in March, all 3 have been about the Zelig package. Please read the posting guide and find out the correct place to send questions about packages. Then you might get an answer. HTH Matthew "Mathew, Abraham T" wrote in message n

Re: [R] regular expression help to extract specific strings from text

On Wed, Mar 31, 2010 at 8:20 AM, Tony B wrote: > Dear all, > > Lets say I have the following: > >> x <- c("Eve: Going to try something new today...", "Adam: Hey @Eve, how are >> you finding R? #rstats", "Eve: @Adam, It's awesome, so much better at >> statistics that #Excel ever was! @Cain & @Abl

[R] library sets: A & EMPTY does not work; gset_intersection(A,EMPTY) works

When using generalized sets from "sets" library A & EMPTY does not work, but gset_intersection(A,EMPTY) works: example code below Is it a bug? Best regards, Ryszard library(sets) support <- universe <- c('a','b','c','d') A <- gset(support=support, memberships=c(0.2, 0.2, 0.9, 0), universe=unive

Re: [R] regular expression help to extract specific strings from text

strapply in gsubfn can extract matches based on content which seems to be what you want: library(gsubfn) f <- function(...) sapply(list(...), paste, collapse = ", ") DF <- data.frame(x, Source = strapply(x, "^(\\w+):", c, simplify = f), Mentions = strapply(x, "@(\\w+)", c, simpli

Re: [R] Problem comparing hazard ratios

Frank, Is there an article that discusses this idea of bootstrapping the ranks of the likelihood ratio chi-square Statistics to assess relative importance of predictors in time-to-event data (specifically Cox PH model)? Thanks, Ravi. -Original Message- From: r-help-boun...@r-project.or

Re: [R] Error "grid must have equal distances in each direction"

M Joshi, I don't know but I guess that some might have looked at your previous thread on 14 March (also about the geoR package). You received help and good advice then, but it doesn't appear that you are following it. It appears to be a similar problem this time. Also, this list is the wrong

Re: [R] ggplot2: Adding points to a density plot

Thanks a lot! This got me started! Joh Dennis Murphy wrote: > Hi: > > Try this: > > library(ggplot2) > movmed <- ddply(movies, .(decade), summarise, med = median(rating)) > m + geom_point(data = movmed, aes(x = med), y = 0, size = 2) > > HTH, > Dennis > > On Wed, Mar 31, 2010 at 4:46 AM, Joh

[R] (no subject)

The answer depends on what your goals are. If you are simply trying to put all the hazard ratios on a similar scale, then I often use the HR comparing risk of a subject at the 25th percentile to one at the 75th percentile. This has already been suggested by David W. But this is not perfect. The

Re: [R] Error "singular gradient matrix at initial parameter estimates" in nls

Try the function called `nls.lm' which is contained in the "minpack.lm" package. See if it solves your problem. Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Corrado Sent: Wednesday, March 31, 2010 9:13 AM Cc: r-help@r-

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

On Mar 31, 2010, at 9:13 AM, Driss Agramelal wrote: > > Hello and thank you both for your answers! > > Dennis, I tried to simply run > > lm(a ~ b) > > after re-importing "a" as a matrix, but I get the following error > message: > > Error in model.frame.default(formula = a ~ b, drop.unused.level

[R] Installing JGR and rJava

Hello, I am running: R version 2.10.0 (2009-10-26) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 on a 64-bit RedHat box. To encourage collegues to use R, I am trying to install "JGR", which depends on "rJava", on the above machine. At first I received and er

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Hello and thank you both for your answers! Dennis, I tried to simply run lm(a ~ b) after re-importing "a" as a matrix, but I get the following error message: Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : invalid type (list) for variable 'a' so maybe I have to spe

[R] regular expression help to extract specific strings from text

Dear all, Lets say I have the following: > x <- c("Eve: Going to try something new today...", "Adam: Hey @Eve, how are > you finding R? #rstats", "Eve: @Adam, It's awesome, so much better at > statistics that #Excel ever was! @Cain & @Able disagree though :(", "Adam: > @Eve I'm sure they'll so

Re: [R] ggplot2: Adding points to a density plot

Hi: Try this: library(ggplot2) movmed <- ddply(movies, .(decade), summarise, med = median(rating)) m + geom_point(data = movmed, aes(x = med), y = 0, size = 2) HTH, Dennis On Wed, Mar 31, 2010 at 4:46 AM, Johannes Graumann wrote: > Hi, > > Consider something like > > library(ggplot2) > > movi

Re: [R] Error "singular gradient matrix at initial parameter estimates" in nls

Dear JN, Bert, 1) It is not a perfect fit. I do not think I have ever said that. I said that an external algorithms fits the model without any problems: with ~ 500,000 data points and 19 paramters (ki in the original equation), it fits the model in less than 1 second. The data are not artifici

[R] position of mismatches in all.equal()?

Dear R users, I would like to compare two dataframes, actually their categorical variables (as factors) only (there are 12, from column 1 to 12). The reason I do that is that I got 2 datasets from two different methods and I would like to be sure that each method used the same data (3D images)

Re: [R] Problem comparing hazard ratios

Thank you, gentlemen. I greatly appreciate your help. -- Michal J. Figurski, PhD HUP, Pathology & Laboratory Medicine Biomarker Research Laboratory 3400 Spruce St. 7 Maloney Philadelphia, PA 19104 tel. (215) 662-3413 __ R-help@r-project.org mailing li

Re: [R] Printing the function t.test() in R

On Mar 31, 2010, at 8:35 AM, Luwis Tapiwa Diya wrote: Dear R Users, I have a question, how does one print an r function like t.test in R (to get the whole function and not just a summary of the environment)? For example if I type the following: t.test function (x, ...) UseMethod("t.t

Re: [R] time series datasets

What kind of time series data? Finance? Econometrics? The 'quantmod' package allows for downloading of stock and FRED data online for analysis. I believe PerformanceAnalytics also has some datasets as well. -c On 3/31/10, vibha patel wrote: > Hello, > > I am searching large multivariate time se

  1   2   >