[R] Announcing new features for the crantastic.org community portal

We’re very pleased to announce an update to http://crantastic.org/, the community site for R packages. Notable new features include: * Easier login - you can use your existing account from Google, Twitter, OpenID, Yahoo!, Flickr, or MyOpenID. * Timelines with recent changes * Task View integration

[R] apply on xts

Hi, I do not understand why after I called apply on a function that returns an xts (getIdvAdjSeries) it returns a matrix whose columns are just numeric value of time series in xts instead of a list of xts objects. Basically, I called the following: apply(matrix(tickers,ncol=1),1,FUN=getDivAdjSer

[R] Can hist plot % of count instead of raw count?

I do not want to use histogram because I have more control over base graphs that R provide. I'm wondering if hist function provides a parameter that can set y axis as % of count instead of raw count? If not, I'd like to know how I can add text or line on histogram function? I don't like the idea o

[R] grid.polygon() + color gradient

Hi, I wonder whether there is a way to generate a polygon (a triangle in my case) with color gradient using grid.polygon() in package grid? I tried something like library(grid) grid.polygon(x=c(0, 0.5, 1), y=c(0.5, 1, 0.5), gp=gpar(col=NA, fill=colorRampPalette(c("green", "lightgray"),

[R] Using by() and stacking back sub-data frames to one data frame

Dear all, I have a code where I subset a data frame to match entries within levels of an factor (actually, the full script uses three difference factors do do that). I'm very happy with the precision with which I can work with R, but since I loop over factor levels, and the data frame is big, the

[R] Lose of decimal when using write.table to text file

Good Day, I have the following: >Date<-c("08/05/08","08/06/08","08/07/08") >Weight<-c(209.4,211.8,210.0) >planned.meal<-cbind(Date,Weight) >planned.meal DateWeight 1 08/05/08, 209.4 2 08/06/08, 211.8 3 08/07/08, 210.0 >write.table(planned.meal, file="plannedMeal1.txt", + quote=FA

Re: [R] Can hist plot % of count instead of raw count?

On Jun 24, 2009, at 9:29 PM, R_help Help wrote: I do not want to use histogram because I have more control over base graphs that R provide. I'm wondering if hist function provides a parameter that can set y axis as % of count instead of raw count? So you want 100 times the proportions t

Re: [R] apply on xts

On Jun 24, 2009, at 9:26 PM, R_help Help wrote: Hi, I do not understand why after I called apply on a function that returns an xts (getIdvAdjSeries) it returns a matrix whose columns are just numeric value of time series in xts instead of a list of xts objects. Basically, I called the fo

Re: [R] tail probability

On Jun 24, 2009, at 8:01 PM, Roslina Zakaria wrote: Hi r-users, If I have 10 data generated from multi skew-t distribution, how do I find the tail probability of alpha=.01? You have given us the tail probability, so I assume you are asking for the the value of x above which the tail

Re: [R] Using by() and stacking back sub-data frames to one data frame

try do.call(rbind, yourByList) hth, Kingsford Jones On Wed, Jun 24, 2009 at 9:34 PM, Stephan Lindner wrote: > Dear all, > > > I have a code where I subset a data frame to match entries within > levels of an factor (actually, the full script uses three difference > factors do do that). I'm very

Re: [R] curvedarrow (some graphics problem)

Hi there, This is now the code % library(grid) vp <- viewport( x = unit(0, "npc"), y = unit(0, "npc"), just = c("left", "bottom"), xscale = c(-1, 1) , yscale = c(-1, 1)) vp2 <- viewport( # probably not needed but I had trouble placing the xaxis x = unit(0,"npc"), y = unit(0.5, "npc"

Re: [R] The gradient of a multivariate normal density with respect toits parameters

Ravi Varadhan: > You may want to look at the paper by Dwyer on "Some applications of matrix > derivatives in multivariate analysis" (JASA 1967), especially the Table 2 > on p. 617. Thanks! That’s a *very* useful paper. Note there also is a small corrigenda available: Corrigenda: Some Applicatio

[R] [R-pkgs] RImageJ initial release

Hello, I've uploaded the initial release (0.0-142) of the RImageJ package to CRAN. The package allows R to use the facilities of the java based image processing and analysis platform ImageJ. See also the announcement on my blog, with a very simple example showing how to use the package: http:

Re: [R] Model fitting with GAM and "by" term

On Tue, 2009-06-23 at 14:18 -0400, Paul Simonin wrote: > Hello R Users, > I have a question regarding fitting a model with GAM{mgcv}. I have data > from several predictor (X) variables I wish to use to develop a model to > predict one Y variable. I am working with ecological data, so have data

Re: [R] gradually switching regression

Hi, Although Bacon&Watts (1971) assume a transition, the aim is to estimate the breakpoint where the linear relationship changes. The start- and end-point of the transition phase are not parameters to estimate ; it is a trick to estimate the model. As a possible alternative, you could have a

Re: [R] Customize axis labels in xyplot

I found indeed that there were no attachements yesterday and corrected the mistake. They are in my previous post. Thanks. -- View this message in context: http://www.nabble.com/Customize-axis-labels-in-xyplot-tp24126788p24180049.html Sent from the R help mailing list archive at Nabble.com. __

[R] Lempel—Ziv complexity

Hi I'm lookoing for a R-plugin to compute the Lempel-Ziv complexity. Does anyone have still heard about this algorithme? Regards -- View this message in context: http://www.nabble.com/Lempel%E2%80%94Ziv-complexity-tp24181807p24181807.html Sent from the R help mailing list archive at Nabble.com

[R] vector size

Hi, I have a data of size 981.1MB(1028707715) and I intend to calculate the length of the data using tapply function in R. I was able to read the data into R but when I tried to use the factor function, I had an error message "Error: cannot allocate vector of size 2.0 Gb". Can anybody tell me what

[R] CRAN: 64 bit OS X Leopard build

Hello, I am after a stable 64 bit binary of R for OS X Leopard (i.e. 2.8). There seems to be the siggestion that thery should be available from CRAN: "leopardBinaries of universal (32-bit and 64-bit) package builds for Mac OS X 10.5 or higher" But when I follow the link there is only a co

Re: [R] vector size

ogbos okike wrote: Hi, I have a data of size 981.1MB(1028707715) and I intend to calculate the length of the data using tapply function in R. I was able to read the data into R but when I tried to use the factor function, I had an error message "Error: cannot allocate vector of size 2.0 Gb". Ca

Re: [R] implementing Maximum Likelihood with distrMod when only the PDF is known

Dear Guillaume, retrospectively, I'm not sure if the decision to have spezial initialize methods is the optimal way to go. In distrMod and our other packages on robust statistics we don't introduce special initialize methods, but use so-called generating functions. This approach has the advant

[R] Found the stable 64bit OS X leopard builds

Hi, I have found the stable 64bit OS X leopard builds here: http://r.research.att.com/ Dan -- ** Daniel Brewer, Ph.D. Institute of Cancer Research Molecular Carcinogenesis Email: daniel.bre...@icr.ac.uk ***

[R] nlme package - unbalanced data and Croissant (2008)

Dear Anthony: please please! I didn't say 'nlme' "[does] not produce sensible / accurate results", nor did I ever mean it. The sentence you quote is an unfortunate mistake on my part, and like most unfortunate mistakes it first went unnoticed and then ended up in a most visible part of the pa

[R] lattice wireframe within a loop ???

Hi, I have the following problem. Calling wireframe within a loop results into an empty window(s) #generate some data temp = expand.grid(A = 1:3,B = 1:3) temp = cbind(temp, y1 = rnorm(9)) temp = cbind(temp, y2 = runif(9)) #plot y1 and y2 in two different windows for(i in 1:2) { wireframe(y1

Re: [R] values in graph points

Hi Romildo, thigmophobe.labels in the plotrix package might do what you want: thigmophobe.labels(x,y,labels=y,...) I have just had a bug report, so if you have any NA or NaN values in your x and y coordinates, wait for version 2.6-2 Jim __ R-help@r

Re: [R] CRAN: 64 bit OS X Leopard build

Daniel Brewer wrote: Hello, I am after a stable 64 bit binary of R for OS X Leopard (i.e. 2.8). There seems to be the siggestion that thery should be available from CRAN: "leopard Binaries of universal (32-bit and 64-bit) package builds for Mac OS X 10.5 or higher" It says "package build"

Re: [R] lattice wireframe within a loop ???

You need to print() lattice plots inside a loop or a function. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwaliteitszorg / Sect

[R] vertical label on x-axis

Hi, have anybody a hint how it's possible plot the labels in the x-axis vertical to get more space. dotplot( value ~ abbreviate(attribute,minlength=3),data=tabstat10,ylim=c(0,1),col="blue",ylab="Value",type="p",par.settings=my.theme) Many thanks, Christian __

Re: [R] CRAN: 64 bit OS X Leopard build

This thread really belongs on R-sig-mac! I am answering here since if people are going to ask here, some general information would be appropriate, but PLEASE follow up only on R-sig-mac. The CRAN Mac OS X page is a little confusing. The main CRAN distribution of R is for Mac OS X 10.4 and la

Re: [R] vertical label on x-axis

I get it : scales=list(x=list(rot=90)) Cheers, Christian Hi, have anybody a hint how it's possible plot the labels in the x-axis vertical to get more space. dotplot( value ~ abbreviate(attribute,minlength=3),data=tabstat10,ylim=c(0,1),col="blue",ylab="Value",type="p",par.settings=my.theme)

[R] subsetting data frame using a vector of column names / values

Hello I have a data frame d with columns "var1", "var2", "var3" Then I have two vectors: columns <- c("var2", "var3") values <- c(0, 1) Is there a compact way to subset the data frame using these two vectors and get the result equivalent to: select (d, var2==0 & var3==1) ? Thank You Blazej

Re: [R] lattice wireframe within a loop ??? - solved

Thank you. ONKELINX, Thierry schrieb: You need to print() lattice plots inside a loop or a function. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometr

[R] lmList

Hi R-helpers: I have a function f<-y~a+b.x and my data is an unbalanced panel data (a dataframe with columns: firm,year,y,x). I need to run that regression by year. I started to subset the data by year and run the lm () function, but I asked R help how to perform this regressions at once. One

[R] gnls : Rho

Hello list: How to extract the value of "Rho" from a gnls() object. I am using gnls() function similar to res <- gnls(y~SSmicmen(),correlation=corCompSymm(form~1|b),data=dat) Thanks in advance, Mahbub. -- Mahbub Latif School of Mathematical Sciences Queen Mary, University of London United Ki

Re: [R] curvedarrow (some graphics problem)

On Jun 24, 2009, at 1:21 AM, casperyc wrote: [quote]> grid.text("t", x=0.3, y=unit(-1, "line"), vp=vp2) Error in valid.units(units) : Invalid unit grid.text("s", 0.8, unit(-1, "line"), gp=gpar(col="red"), vp=vp2) Error in valid.units(units) : Invalid unit[/quote] Hi, I have loaded the [gr

[R] bwimage

Dear Sir/Madam, I am using R-Image in an MSc statistics dissertation and have found, in various EBImage manuals a functions - bwimage, which I would like to use, however I can't find it within my version of R. Could you tell me if the function has changed or what I need to do to be able to use i

[R] lda.cv help

i am a new user of R working on "Cost sensitivity for Classical classification functions" i have codes for cross validation for both linear and quadratic classifier but i cant not very cost of missclassification. please any body with help should help me thanks         ARIYO,OLUDARE SAMUEL Statis

[R] bwlabel

I am using R-Image in an MSc statistics dissertation and have found, in various EBImage manuals a functions - bwlabel, which I would like to use, however I can't find it within my version of R. Could you tell me if the function has changed or what I need to do to be able to use it. -- View this

[R] I meant bwlabel

Sorry I meant bwlabel on this message Alice --- On Wed, 24/6/09, alicedug...@btopenworld.com wrote: From: alicedug...@btopenworld.com Subject: bwimage To: r-help@r-project.org Date: Wednesday, 24 June, 2009, 12:34 PM Dear Sir/Madam, I am using R-Image in an MSc statistics dissertation

Re: [R] Using a variable for mysql/database Query

Hello, Thanks! Join is not possible, as the database is simply way to large. When we exactly now, what we are looking for maybe, but not yet… By now, I figured the paste command - not familiar with R…, however, it only seems to work for one - the last - entry I am using. Unfortunately, the help

Re: [R] subsetting data frame using a vector of column names / values

Try this: subset(d, eval(parse(text = paste(paste(columns, values, sep = "=="), collapse = " & " On Wed, Jun 24, 2009 at 9:16 AM, Blazej Krzeminski wrote: > Hello > > I have a data frame d with columns "var1", "var2", "var3" > > Then I have two vectors: > columns <- c("var2", "var3") > valu

[R] Test for Temporal autocorrelation in GPS obtained points

Hi, I have data collected from PTT-GPS devices tagged to birds. These data typically are delivery in time slots, comprising the day-light period for example. I'm interested in testing the temporal correlation of observations (points) and eventually exclude locations falling below a given corre

Re: [R] bwimage

Dear Alice, EBImage is a general-purpose image processing package for R, able to read, write, process and analyze images, but is part of the Bioconductor project. Please post your future questions on the BioC mailing list (see http://www.bioconductor.org/docs/mailList.html). There is no func

[R] Odp: subsetting data frame using a vector of column names / values

Hi r-help-boun...@r-project.org napsal dne 24.06.2009 14:16:15: > Hello > > I have a data frame d with columns "var1", "var2", "var3" > > Then I have two vectors: > columns <- c("var2", "var3") > values <- c(0, 1) > > > Is there a compact way to subset the data frame > using these two vectors

[R] parallel rotated strips with color gradient

Hi, I want to produce two parallel rotated strips with color gradient. So far, the sample strip is something produced by this: pushViewport(viewport(x = unit(0.638, "npc"), y =unit(0.386, "npc"), width=.62, height=0.006, angle=137.2)) grid.rect(y=100:1/100, just="top", gp=gpar(col=NA

[R] label sorting x-axis

Thanks, how is it possible that the x-axis labels getting the row sorting instead the alphabetically one? Could i just printed only every n (i.e. 3rd) label? Christian df <- data.frame(value=runif(26),le=letters) df <- df[order(df[,1]),] dotplot(value ~ le,data=df) Maybe you want to have

Re: [R] label sorting x-axis

Try this: dotplot(value ~ factor(le, levels = le),data=df) On Wed, Jun 24, 2009 at 10:24 AM, Christian Schulz wrote: > Thanks, how is it possible that the x-axis labels getting the row sorting > instead the alphabetically one? > Could i just printed only every n (i.e. 3rd) label? > > Christi

Re: [R] curvedarrow (some graphics problem)

You're right, I meant to write "lines"not "line". The strange thing is, although "line" isn't listed in ?unit, it doesn't return an error on my machine, > unit(-1, "line") [1] -1line > unit(-1, "lines") [1] -1lines Reading from http://svn.r-project.org/R/trunk/src/library/grid/src/unit.c "l

[R] Boxplots: side-by-side

Dear R-sians.. I am trying to plot boxplots with side-by-side option.. I tried some of the posted suggestions and could not make it work due to unequal sizes of categories... e.g. weekly measured water depth values are categorized into 5 levels based on their values such measurement is again cate

Re: [R] How to import timestamps from emails into R

One last tiny problem: How do I add months to the scale? It currently just has years http://school.thomaslevine.org/mywall.png Thanks again Tom On Sat, Jun 20, 2009 at 12:14 PM, Thomas Levine wrote: > I wasn't really thinking that far ahead; plot tries to do something, so I > figured I'd try th

Re: [R] curvedarrow (some graphics problem)

I get: > unit(-1, "line") Error in valid.units(units) : Invalid unit > unit(-1, "lines") [1] -1lines I have an older version of R (2.8.1) and of grid (2.8.1). The source for grid 2.8.1 does not have pseudonyms. http://svn.r-project.org/R/branches/R-2-8-branch/src/library/grid/src/unit.c -- Da

Re: [R] driver file

Do you mean something like source? Type ?source for help --- On Tue, 6/23/09, Derek Lacoursiere wrote: > From: Derek Lacoursiere > Subject: [R] driver file > To: r-help@r-project.org > Received: Tuesday, June 23, 2009, 3:05 PM > > Hi, > > How can I, from a single "driver" file, source othe

Re: [R] How to import timestamps from emails into R

See the example in plot.zoo labelled Fancy X Axis: library(zoo) example(plot.zoo) ?plot.zoo On Wed, Jun 24, 2009 at 10:11 AM, Thomas Levine wrote: > One last tiny problem: How do I add months to the scale? It currently just > has years > http://school.thomaslevine.org/mywall.png > > Thanks again

[R] Second LondonR Group meeting

Hello, Mango is pleased to announce that the next LondonR Group meeting will be held on 21st July. At this meeting we plan to have a feedback session from people who attended useR! in Rennes. To help us plan this, if you are attending at Rennes next month and would be interested in giving some

[R] rgdal package -- Error in installing package

R-help, I'm trying to install the rgdal package under Ubuntu but I get the following warning: > install.packages("rgdal") Warning in install.packages("rgdal", configure.args = "--with-gdal-modules=/usr/local/lib") : argument 'lib' is missing: using '/home/luisridaocruz/R/x86_64-pc-linux-gnu-

Re: [R] Boxplots: side-by-side

reproducible code. On Wed, Jun 24, 2009 at 8:53 AM, Santosh wrote: > Dear R-sians.. > > I am trying to plot boxplots with side-by-side option.. I tried some of the > posted suggestions and could not make it work due to unequal sizes of > categories... > > e.g. > weekly measured water depth values

[R] order by decerasing 1st variable and increasing 2nd variable

Hello, I have a data.frame which I would like to sort with the primary key decreasing while the secondry key is increasing e.g. x <- data.frame(One=c(1,1,1,2,2,3,4,5),Two=c(2,3,1,2,3,3,3,3)) I would like to order it so it looks like this: One Two 8 5 3 7 4 3 6 3 3 4 2 2 5 2

[R] Qualified parameters in SOAP body using .SOAP

Hello, I am trying to reach a web service using the SOAP package. I succeeded calling the web service, but not sending parameters to it. After much research and tries, I think I found that the problem lies in the namespace including the parameters in the SOAP body. In short, my question is:

Re: [R] order by decerasing 1st variable and increasing 2nd variable

How about this: > x[order(x$One,-x$Two,decreasing=T),] One Two 8 5 3 7 4 3 6 3 3 4 2 2 5 2 3 3 1 1 1 1 2 2 1 3 2009/6/24 Daniel Brewer : > Hello, > > I have a data.frame which I would like to sort with the primary key > decreasing while the secondry key is increa

Re: [R] rgdal package -- Error in installing package

$ sudo R > install.packages(rgdal) worked for me on my Ubuntu machine but, I already had GDAL library installed with $ sudo apt-get install libgdal1-1.5.0 -- B. H. Braswell (Rob) University of New Hampshire On Jun 24, 2009, at 9:24 AM, Luis Ridao Cruz wrote: R-help, I'm trying to install

[R] [Classification] lifting score in R

Hi all, Could anybody give me some pointers to Cross Validation using Lifting Score as error function, as commonly used in data-mining and classification field in marketing and e-commerce research? Thanks! __ R-help@r-project.org mailing list https://s

Re: [R] Boxplots: side-by-side

Try this Rich tmp <- data.frame( y=rnorm(100), category=rep(factor(letters[1:5]),each=20), level=rep(factor(0:1), length=100)) tmp table(tmp[,2:3]) tmp$y[with(tmp, category=="a" & level=0)] <- NA tmp$y[with(tmp, category=="a" & level==0)] <- NA tmp$y[with(tmp, category=="e" & level==1)] <- NA t

[R] how to undo automatic loading of packages?

I wanted to try out package distrMod, so I did > install.packages('distrMod') > library(distrMod) and played around, saved and quit. Now whenever I start up in this directory, I get distr and lots of other stuff loaded and lots of messages. How do I keep it from automatically loading, other tha

[R] building R packages

Hello, I tried to build an own R package on Windows XP but get an error which I can't solve. I called my package "pack". It works to creat the file "pack" with package.skeleton(). But when I try to compile it with Rcmd build --binary pack I obtain the following error: IO::Seekable::seek missi

Re: [R] The gradient of a multivariate normal density with respecttoits parameters

Karl, I have written a function to compute the gradient of a multivariate normal density (derivatives at a given point, with respect to the parameters of the density). It is in the atttached file. I would appreciate any feedback. There are couple of interesting points: 1. Note is that there

Re: [R] subsetting data frame using a vector of column names / values

This "one liner" works great! Thanks (for all replies) > > Try this: > > subset(d, eval(parse(text = paste(paste(columns, values, sep = "=="), > collapse = " & " > >> >> Hello >> >> I have a data frame d with columns "var1", "var2", "var3" >> >> Then I have two vectors: >> columns <- c("var2

Re: [R] how to undo automatic loading of packages?

On Wed, 24 Jun 2009, dav...@rhotrading.com wrote: I wanted to try out package distrMod, so I did install.packages('distrMod') library(distrMod) and played around, saved and quit. Now whenever I start up in this directory, I get distr and lots of other stuff loaded and lots of messages. How

Re: [R] [Classification] lifting score in R

Maybe the packages caret,RWeka and ROCR are usefuel starting points. Cheers, Christian Hi all, Could anybody give me some pointers to Cross Validation using Lifting Score as error function, as commonly used in data-mining and classification field in marketing and e-commerce research? Thanks!

[R] Random Forest Variable Importance Interpretation

Hi I am trying to explore the use of random forests for regression to identify the important environmental/microclimate variables involved in predicting the abundance of a species in different habitats, there are approx 40 variable and between 200 and 500 data points depending on the dataset. I

Re: [R] [Classification] lifting score in R

Michael, a lift chart for evaluating binary scoring classifiers, as I understand it, plots... lift score: P(Yhat = + | Y = +)/P(Yhat = +) against rate of rate of positive predictions: P(Yhat = +). ...across the continuum of possible cutoffs. If you want to do this, here is how you would do this

[R] "by" question

Hello all I have a big data frame and I regularly want to break it down into subsets, calculate some new data, and add it back to the data frame. At the moment my technique seems a bit ugly and embarrassing. Something like: result <- by(mydata, mydata$some_factor, function (x) { # do something

Re: [R] how to undo automatic loading of packages?

That was unintentional: the mail server crashed and apparently sent the just-started reply. It is nothing to do with 'Autoloads'. What is most likely if that you have saved a reference to a namespace such as 'distrMod' in your .RData file. If that it is the case, when .RData is loaded it wil

Re: [R] "by" question

How about something like this: > x id data 1 1 0.7773207 2 3 0.9606180 3 2 0.4346595 4 3 0.7125147 5 2 0.344 6 2 0.3253522 7 2 0.7570871 8 3 0.2026923 9 3 0.7111212 10 2 0.1216919 > # compute running sum for each ID > x$run <- ave(x$data, x$id, FUN=cumsum) > x id

Re: [R] how to undo automatic loading of packages?

I removed all of the objects that were from the examples I tried, saved, and quit. R starts up cleanly now. Yeah! That was annoying behavior, and I would recommend that the authors change it, please. Thank you for your to-the-point advice, Professor Ripley! -- David ps - Sorry for the lengthy dis

Re: [R] [Classification] lifting score in R

Thanks a lot! Yet is there a way to incorporate the lifting score into Cross Validation, not just a plot? Thanks again! On Wed, Jun 24, 2009 at 9:07 AM, Tobias Sing wrote: > Michael, > > a lift chart for evaluating binary scoring classifiers, as I > understand it, plots... > > lift score: P(Yhat

[R] How to avoid ifelse statement converting factor to character

Hi R-helpers, Please see the below R output. The problem is that after running the ifelse statement, data$SOCIAL_STATUS is converted from a factor to a character. Is there some way I can avoid this conversion? Thanks in advance, Mark Na > str(data) 'data.frame': 2100 obs. of 11 variables: $ D

Re: [R] lmList

Cecilia, this is impossible to tell from the information you provide. The simulated data below... ##START year=rep(1:10,10) year.fe=rep(rgamma(10,10,2),10) firm.fe=rep(rnorm(10,0,3),each=10) x=rnorm(100,0,3) e=rnorm(100,0,3) y=10+firm.fe+year.fe+x+e newdata=data.frame(y,x,e,firm,year,firm.fe,year.

[R] Problem installing RODBC package under UBUNTU

R-help, I'm now on the process of installing RODBC on UBUNTU but an error occured. I copy-paste the output: > install.packages("RODBC") Warning in install.packages("RODBC") : argument 'lib' is missing: using '/home/luisridaocruz/R/x86_64-pc-linux-gnu-library/2.8' trying URL 'http://cran.ii.uib

[R] change the height or scale of the y axis

Hallo, All, I have a question about changing the height or scale of the y axis. When I use following two R codes, I can get two plots. Please look at the y axes, the number of indices (x1, x2, …) on the y axis in the first plot is smaller than that in the second plot, and hence the space between

Re: [R] How to avoid ifelse statement converting factor to character

You can work with levels of factor: levels(DF$SOCIAL_STATUS)[DF$MALE > 4 & DF$SOCIAL_STATUS == "B"] <- "C" On Wed, Jun 24, 2009 at 1:34 PM, Mark Na wrote: > Hi R-helpers, > > Please see the below R output. > The problem is that after running the ifelse statement, data$SOCIAL_STATUS > is convert

Re: [R] "by" question

That seems to work. I should add that to make "ave" work like "by" one can do: mydata$newvar <- ave(1:nrow(mydata), mydata$some_factor, FUN= function (x) { x <- ds[x,] # ... etc... }) Thanks! David [[alternative HTML version deleted]] __ R-h

Re: [R] How to avoid ifelse statement converting factor to character

Explicitly convert it to a factor... data$SOCIAL_STATUS<-factor(ifelse(data$SOCIAL_STATUS=="B" & data$MALE>4, "C", data$SOCIAL_STATUS)) However, note that this would, in general, change the levels attribute to the levels actually present in the converted vector. If you wish to (and it makes sense

Re: [R] Problem installing RODBC package under UBUNTU

Luis, On 24 June 2009 at 16:35, Luis Ridao Cruz wrote: | R-help, | | I'm now on the process of installing RODBC on UBUNTU | but an error occured. I copy-paste the output: You can simply do $ sudo apt-get install r-cran-rodbc As for the error: | > install.packages("RODBC") [...] | co

[R] Coxph frailty model counting process error X matrix deemed singular

Hello, I am currently trying to simulate data and analyze it using the frailty option in the coxph function. I am working with recurrent event data, using counting process notation. Occasionally, (about 1 in every 100 simulations) I get the following warning: Error in coxph(Surv(start, end,

Re: [R] Qualified parameters in SOAP body using .SOAP

Hi Olivier Olivier Cailloux wrote: Hello, I am trying to reach a web service using the SOAP package. I succeeded calling the web service, but not sending parameters to it. After much research and tries, I think I found that the problem lies in the namespace including the parameters in the S

Re: [R] I meant bwlabel

aliceduggan wrote: > > Sorry I meant bwlabel on this message > [...] > I am using R-Image in an MSc statistics dissertation and have found, in > various EBImage manuals a functions - bwimage, which I would like to use, > however I can't find it within my version of R. > Could you tell me if the

[R] Goodness of fit test / pseudo r^2 measure for Zero Inflated Model

Hi I have been using a Zero-Inflated negative binomial model fitted using the pscl zeroinfl command but I would like to extract a goodness of fit measure are there any suitable pseudo R^2 measures available for this type of analysis to try and assess the amount of variation in the data explain

[R] How to recursively build a binary tree with R script?

Greetings! Can someone provide a simple script for a R function that recursively builds a binary tree. I am most familiar with C and pass by reference, but I think R is like Fortran and pass by value. - Mike Beddo __ R-help@r-project.org mailing list

[R] agrep behavior

Dear list - I am a bit puzzled by the behavior of agrep: The following command finds a match: agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.6) But if I _increase_ the maximum distance to .9 it fails: agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.dist

[R] subplot()

Hi, I'm trying to get a color-gradient triangle out of a large plot. Wonder can subplot() generate a triangle? So far, I'm only getting rectangles. Great appreciation for any reply!! Kexin __ R-help@r-project.org mailing list https://stat.ethz.ch/mail

[R] Reading 4-byte float into R from binary file

A simple question, I hope I have a binary file format that I want to slurp into R using readBin. There are a couple of fields in the file that are 4-byte floats. Since R has no "float" data type (4-byte), I am wondering how to go about this task. I could simply do all the work in C, but the

Re: [R] "by" question

You might also want to look at the plyr package, http://had.co.nz/plyr. In particular, ddply + transform makes these tasks very easy. library(plyr) ddply(mtcars, "cyl", transform, pos = seq_along(cyl), mpg_avg = mean(mpg)) Hadley On Wed, Jun 24, 2009 at 11:48 AM, David Hugh-Jones wrote: > That

Re: [R] How to avoid ifelse statement converting factor to character

On Wed, Jun 24, 2009 at 12:34 PM, Mark Na wrote: > The problem is that after running the ifelse statement, data$SOCIAL_STATUS > is converted from a factor to a character. > Is there some way I can avoid this conversion? I'm afraid that ifelse has very bizarre semantics when the yes and no argument

Re: [R] Reading 4-byte float into R from binary file

On 6/24/2009 2:41 PM, Sean Davis wrote: A simple question, I hope I have a binary file format that I want to slurp into R using readBin. There are a couple of fields in the file that are 4-byte floats. Since R has no "float" data type (4-byte), I am wondering how to go about this task. I

Re: [R] agrep behavior

Unable to reproduce: > agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.6) [1] 1 > > agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.89) [1] 1 > agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.9) [1] 1 > agrep("Sta

Re: [R] How to recursively build a binary tree with R script?

On 6/24/2009 2:02 PM, Mike Beddo wrote: Greetings! Can someone provide a simple script for a R function that recursively builds a binary tree. I am most familiar with C and pass by reference, but I think R is like Fortran and pass by value. Here's a tree with a random depth: makeTree <- fu

[R] Why can't I use ADF.test?

Greetings! I'm trying to use R to test the unit root for a univariate data. By this link: http://rss.acs.unt.edu/Rdoc/library/uroot/html/ADF.test.html it tells me that I can use the function ADF.test(). However, when I tried this in R, I got this message: "Error: could not find function "ADF.te

Re: [R] Why can't I use ADF.test?

Since you provided no code, the following is just a guess, Try: install.packages("uroot") library(uroot) Then try your analysis again. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of DongHongwei Sent: Wednesday, June 24, 2009 3:

Re: [R] agrep behavior

Thanks David, for trying to replicate. May be it is a 2.9 problem? I should have attached session info: > agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.9) integer(0) > agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=.7) [1] 1 sessionInfo() R vers

Re: [R] subplot()

The subplot function (at least the one in the TeachingDemos package and the copy in Hmisc) only does rectangles. Clipping regions are also limited to rectangles currently. There are external tools (gimp, imagemagick) that may me of use or triangles and other shapes can be approximated by enoug

  1   2   >