Re: [R] Custom ticks on x axis when dates are involved

On 11/25/2010 06:27 AM, Monica Pisica wrote: ... Now the graph looks very close to what i want, but i know that my ticks actually are not exactly at 01/01/ as i would like, although i suppose my error is not that much in this instance. However i would really appreciate if i can get the ti

[R] Reg : Using an existing JDBC connection in RJDBC

Hi R-helpers, I have an application in which the data is read from the datasource and stored in a table using a JDBC connection. I need to read the data from those tables using the same JDBC connection object created.I saw that in RJDBC we need to create a new connection object to connect to th

[R] Reg : Using an existing JDBC connection in RJDBC

Hi R-helpers, I have an application in which the data is read from the datasource and stored in a table using a JDBC connection. I need to read the data from those tables using the same JDBC connection object created.I saw that in RJDBC we need to create a new connection object to connect to tha

[R] How to change value of y axis from log relative Hazard to relative Hazard

http://r.789695.n4.nabble.com/file/n3058505/file.csv file.csv Hi, Rusers I have a problem in making a rcspline.plot with a Hmisc package. My data is in the upload attachment. My programme as follows: library(Hmisc) A<-read.csv("file.csv",header=TRUE) attach(A) rcspline.plot(factor,Ti

Re: [R] overlap cdf plots and add colors and etc

Hi Jorge, I tried but still it does not work.  Thank you for your time. From: Jorge Ivan Velez Cc: r-help@r-project.org Sent: Thu, November 25, 2010 4:46:37 PM Subject: Re: [R] overlap cdf plots and add colors and etc Hi Roslina, Try par(mar=c(4,4,2,1.2)

Re: [R] Determining a basal correct count

Phil, I wanted to thank you for this solution. I was working on other projects for the past couple of weeks, but I was finally able to come back to this and get it to work. The key was excluding a couple of columns on the right margin of the data base, the syntax wouldn't process those columns,

Re: [R] overlap cdf plots and add colors and etc

Hi Roslina, Try par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(rnorm(100))) plot(ecdf(rnorm(100)), add = TRUE, col = 2) HTH, Jorge On Thu, Nov 25, 2010 at 12:18 AM, Roslina Zakaria <> wrote: > Hi r-users, > > I would like to overlap 2 ecdf plots. > > I tried this below and

[R] Heirarchial clustering based heatmap with user defined colours

Hi guys is it possible to cluster many different sets of values together by using Heir.clustering and generate a heatmap ? input looks like Tf1_peak1_value 20 Tf2_peak1_value 30 Tf3 Tfn . Tf1_peak2_value 40 Tf2_peak2_value 60

[R] overlap cdf plots and add colors and etc

Hi r-users, I would like to overlap 2 ecdf plots.  I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail.

[R] Lifting Wavelet Transform

Hi R users   Thanks in advance   Is lifting wavelet transform implemented in R? If so, which package or codes  can be used for performing that?    assaed...@yahoo.com   Thanks     [[alternative HTML version deleted]] __ R-help@r-project

[R] Help on running regression by grouping firms

Hi there, I have a huge data set with multiple firms years and other firm characteristics. I want to run a regression on the dependent variable and other explanatory variables and calculate the residual terms by grouping the firms in same year and same industry. What I want to do is to divid

Re: [R] Replacing 1 plot in layout

On Nov 24, 2010, at 8:24 PM, Kohleth Chia wrote: Hi, layout(1:2) # split the screen into 2 regions plot(1:2) # plot in region 1 plot(1:2) # plot in region 2 par(mfg=c(1,1)) # Then I want to update the plot in region 1 You cannot "replace" it with another plot (lattice is set up to doe tha

Re: [R] Time series data

Hi: Perhaps something like this, assuming DATE is a Date object (try str(plcm60) to check) - if not, you need to use as.Date() to convert. jandays <- data.frame(DATE = seq(as.Date('2003-01-01'), by = 'days', length = 23)) merge(jandays, plcm60, by = 'DATE', all.x = TRUE) HTH, Dennis On Wed, Nov

Re: [R] Execute SQL Stored Procedure in R

On Nov 24, 2010, at 9:07 PM, Suphajak Ngamlak wrote: Dear all, I would like R to retrieve a table resulting from execute sql stored procedures. What function can I use? Can RODBC do it. You might try looking at the package RSQLiite.extfuns I tried channel1 <- odbcConnect("ptsecmstqa0

Re: [R] delete-d jackknife

Hi: You could always create a matrix of nCd combinations of size d from n = length(data) and apply your delete-d jackknife to that. Here's a simple version for delete-2: # data vector > u <- rnorm(10) > library(combinat) Attaching package: 'combinat' The following object(s) are masked from 'pack

[R] Execute SQL Stored Procedure in R

Dear all, I would like R to retrieve a table resulting from execute sql stored procedures. What function can I use? Can RODBC do it. I tried channel1 <- odbcConnect("ptsecmstqa01-alpha") query <- paste("execute DB.dbo.usp_test") Data<- sqlQuery(channel1, query) However, it returned blank. Be

[R] file reading problem unique to windows. Err info: Error in file(file, ifelse(append, "a", "w")). cannot open the connection

Dear List I asked this question before, got some tips but can't get it solved. it is clear now that this problem only occurs when run on windows (I tested it on windows XP) nothing wrong at all when run on unix. unfortunately, sometimes I have to run it on windows, I appreciate any suggestion on h

[R] Time series data

I have a stock price dataset a snippet of which is: > plcm60[1:15, c(1,3,4,5,6,7)] DATE BIDLO ASKHIPRC VOL RET 1 1/2/03 9.450 9.79 9.700 1531819 0.018907 2 1/3/03 9.670 9.94 9.940 1582192 0.024742 3 1/6/03 9.830 10.05 9.960 1843298 0.002012 4 1/7/03 9.835

[R] Replacing 1 plot in layout

Hi, layout(1:2) # split the screen into 2 regions plot(1:2) # plot in region 1 plot(1:2) # plot in region 2 par(mfg=c(1,1)) # Then I want to update the plot in region 1 plot(1.5,1.5) But the this last plot overlays with my 1st plot. Q: How do I replace my first plot with my last? NOTE ##

[R] delete-d jackknife

Hi dear all, Can aynone help me about delete-d jackknife usually normal jackknife code for my data is: n <- nrow(data) y <- data$y z <- data$z theta.hat <- mean(y) / mean(z) print (theta.hat) theta.jack <- numeric(n) for (i in 1:n) theta.jack[i] <- mean(y[-i]) / mean(z[-i]) bias <- (n - 1) * (mea

Re: [R] Newbie Help Please

Hi, It is really hard to help much without some example data. If you create a minimal example, we can probably help more. My guess is that the data is not in the proper format to do what you want. Take a look at the results of str(YourData). Whatever you are doing boxplots by (example males an

[R] Newbie Help Please

Hi I am doing a project at school for another company. We are very new to R and are close to a deadline and nowhere near finished. With that said, our data files are csv format, and have 3 columns (date, time, precipitation). When I do a boxplot it displays the graph all weird. We combined date an

Re: [R] Aligning Grid Lines with Date Ticks in Trellis Plot

On 2010-11-24 06:59, Elliot Joel Bernstein wrote: I know this issue has been discussed before, but I hoped the advent of pretty.Date would resolve it. In the following code, the first plot produces misaligned vertical grid lines, while in the second plot they are properly aligned. Is there any

Re: [R] Create new string of same length as entry in dataframe

Simple, but was taking me a long time to google. Thanks Phil, happy Thanksgiving. -- View this message in context: http://r.789695.n4.nabble.com/Create-new-string-of-same-length-as-entry-in-dataframe-tp3058232p3058275.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Create new string of same length as entry in dataframe

Suppose you want to replace using all b's: thestring = 'bb' t2$V3 = substring(thestring,1,t2$V2) t2$V3 [1] "bbb""" "bb" "b" "bb" For the second case, thestring = paste(rep('bad',5),collapse='') thestring [1] "badbadbadbadbad" t2$V3 = substring(the

[R] Create new string of same length as entry in dataframe

I suspect that this is simple, but thanks in advance for any advice... I have a dataframe, t2: V1 V2 aaa 3 4 aa6 a1 aa 2 V2 is the length of the string in V1 using nchar(as.character(t1$V1)) I'd like to create a third column, that cont

Re: [R] Population abundance, change point

On Wed, 2010-11-24 at 13:00 -0500, Jonathan P Daily wrote: > I understand that smoothing splines produce continuous models, however the > end product of SiZer is not a single model that is then used in any > predictive manner. Rather, the end product is a map of potential > changepoint locations

[R] Seeking advice on dynamic linear models with matrix state variable.

  Hello, fellow R users,   I recently need to estimate a dynamic linear model in the following form:   For the measurement equation:   Y_t = F_t * a_t + v_t   where Y_t is the observation. It is a 1 by q row vector for each t. F_t is my forecasting variable. It is a 1 by p row vector. a_t is my st

Re: [R] looking for the RMySQL package for R 2.12.0 under XP

On 25/11/2010 4:55 a.m., Henrique Dallazuanna wrote: You need compile it from source code. An alternative is to use RODBC. There are instructions on my page here: http://www.stat.auckland.ac.nz/~dscott/779/index.php Look in the MySQL section (obviously). You may need to pester your admin guy

Re: [R] Getting user input without hitting enter

Thanks! that's what I need. KC ^_^ On 25/11/2010, at 6:45 AM, Greg Snow wrote: > This is a user interface issue. The standard command line user interfaces > all wait for the user to hit enter before sending the information to the > parser, so this cannot be done using the standa

Re: [R] Problem setting the number of digits in xtable

Thanks Brian, You're right ! They were strings and I thought they were numerics. I just had to force the conversion and Everything is now fine. Thanks for your help. It was really stupid from me, but I couldn't imagine that these colomns were interpreted as strings. Sylvain Le 24 nov. 2010 à

Re: [R] diverting output from nested loops

Lara Poplarski [Wed, Nov 24, 2010 at 06:14:58PM CET]: > for (K in ...){ > > ... create object ... > > topDirName <- as.character(paste("topName", K, sep="")) > topDirMake <- dir.create(topDirName) You may want to insert setwd(topDirName) here. > > for (L in ...) { > >

[R] RJDBC doesn't load Oracle driver?

I am having trouble loading Oracle JDBC driver into RJDBC. Basically, i installed RJDBC and used it to connect to MySql. However, when i try to load JDBC Oracle driver, it gives me an error message > drv<-JDBC("oracle.jdbc.driver.OracleDriver", "/pathToDriver/ojdbc6.jar") Error in .jfindClass(as

Re: [R] Tinn R to R problem

Hello, All you need are explained at Help/Main/User guide/HTML from Tinn-R menu. Basically: 1. Close R 2. Menu R/Configure/ Pemanent (Rprofile.site) - Tinn-R will generate a small script inside of the file Rprofile.site (located on the folder 'etc' where R is istalled) 3. Reboot R 4. It is al

[R] RJDBC doesn't load Oracle driver?

After installing RJDBC on RedHat Linux and using it successfully with mySQL JDBC driver i am trying to use it with Oracle database. I use the JDBC driver that works fine elsewhere (use it in DBVizualizer). However, when i try to load the driver, drv<-JDBC("oracle.jdbc.driver.OracleDriver", "/hom

Re: [R] Problem setting the number of digits in xtable

On 11/22/2010 9:16 AM, wphantomfr wrote: DEar list members, I am currently using Sweave with LaTeX which is great. I can use xtable for formatting outp of tables but I have a problem setting the number of decimals in xtables when used with dataframe. I have found an example on the net ith matr

Re: [R] Custom ticks on x axis when dates are involved

Try this: d <- c("4/6/1984", "9/29/1984", "1/19/1985", "3/27/1986", "10/3/1987", "10/8/1987", "1/28/1988", "12/16/1989", "10/11/1991", "10/5/1992", "11/15/1995", "4/7/1996", "10/3/1997", "2/28/1998", "10/11/2000", "10/30/2001", "2/27/2002", "12/28/2002", "10/20/2003", "10/20/2003", "10/20/2003", "

Re: [R] Tobit model on unbalanced panel

Hi Terry On 24 November 2010 15:34, Terry Therneau wrote: > For tobit regression, also see the last example of help(survreg) in the > survival package. How does the survreg() function of the "survival" package account for the panel structure of the data? Best wishes, Arne -- Arne Henningsen h

Re: [R] Tobit model on unbalanced panel

On 23 November 2010 23:32, Achim Zeileis wrote: > On Tue, 23 Nov 2010, Liang Peng wrote: > >> Appreciate any suggestions regarding how to fit an unbalanced panel data >> to >> a Tobit model using R functions. I am trying to analyze how real estate >> capital expenditures (CapEx) are affected by ma

Re: [R] apply over list of data.frames

Fabulous! I'll start with trying out Reduce, and then explore do.call or converting it to an array if I need to. Thank you Dimitris, Baptiste, and Josh for the very quick help! Best, Tim >>> Dimitris Rizopoulos 11/24/2010 2:43 PM >>> try this: DF.lis <- list( one = data.frame(x = c(1,

[R] Is there an equivalent to predict(..., type="linear") of a Proportional hazard model for a Cox model instead?

Hi all, Is there an equivalent to predict(...,type="linear") of a Proportional hazard model for a Cox model instead? For example, the Figure 13.12 in MASS (p384) is produced by: (aids.ps <- survreg(Surv(survtime + 0.9, status) ~ state + T.categ + pspline(age, df=6), data = Aidsp)) zz <- predi

[R] Hospital ANOVA/ANCOVA problem

Hi everyone, I've only been using R for a week or so now, and am now required to work with some pretty mixed data at my hospital. Basically, we're looking to see if hypertension status (nominal) affects cardiac function (multiple & continuous). In order to do this, we've collected multiple echoc

[R] Change the class of columns in a data frame

Hi. First of all, excuse me if I do any mistakes, but English is not a language I use very often. I have a data frame with numbers. A small part of the data frame is this: nominal ordinal 2 2 2 1 2 1 2 2 So, I want to use the gower distanc

[R] Laplace Approximation

Does anyone have any R code that shows how to do a Laplace Approximation? I know there are a variety of these numerical approximation algorithms and I'm pretty open at this point, I'm just curious how it's approximated in R code. I have seen some functions in packages, but I think they all call

Re: [R] Par() Variables, Timing, and Resizing

I would suggest using the layout function instead of the mfcol and pin with par. The layout function gives you more control of the size of areas to be plotted to. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original

Re: [R] Getting user input without hitting enter

This is a user interface issue. The standard command line user interfaces all wait for the user to hit enter before sending the information to the parser, so this cannot be done using the standard command line interface. You could rewrite the source code, but that is probably overkill. The bet

Re: [R] apply over list of data.frames

Hi Tim, It may not be possible, but if you can use an array, it is super easy: ## convert to array Sue <- array(unlist(x), dim = c(3, 3, 2)) ## use apply() on first two dimensions (collapse across 3rd) apply(X = Sue, MARGIN = c(1, 2), FUN = mean) Cheers, Josh On Wed, Nov 24, 2010 at 11:37 AM,

Re: [R] apply over list of data.frames

try this: DF.lis <- list( one = data.frame(x = c(1,2,3), y = c(1,2,2), z = c(3,1,1)), two = data.frame(x = c(2,5,2), y = c(2,3,1), z = c(4,1,2)) ) Reduce("+", DF.lis) / length(DF.lis) I hope it helps. Best, Dimitris On 11/24/2010 8:37 PM, Tim Howard wrote: R users, This probably in

Re: [R] apply over list of data.frames

Hi, Try this, do.call(`+`, x) / length(x) HTH, baptiste On 24 November 2010 20:37, Tim Howard wrote: > R users, > This probably involves a simple incantation of one of the flavors of apply... > that I can't yet figure out. Consider a list of data frames. I'd like to > apply a function (mean

Re: [R] Cannot update R

Thanks Phil, I may not be as versed in UNIX as you. Can I ask you to show me the BASH command sequence. I did find the original package, I must have downloaded it mannualy. But I do not know where I installed it, how can I find that? Thanks Alon -- View this message in context: http://r.7

Re: [R] Cannot update R

Thank you Josh Updating, you suggest I do it through the GUI of the package manager? unless you know how to do it from the the command line I tried to get rid of rjags but it still says it cannot find it, How do I locate the library? Thanks Alon -- View this message in context: http://r.78

[R] apply over list of data.frames

R users, This probably involves a simple incantation of one of the flavors of apply... that I can't yet figure out. Consider a list of data frames. I'd like to apply a function (mean) across the list and return a dataframe of the same dimensions where each cell represents the mean of that cell a

Re: [R] Change the class of columns in a data frame

you can do it like this: # example data DF <- data.frame(nominal = rep(2, 5), ordinal = c(1,2,2,1,2)) DF$nominal <- factor(DF$nominal) DF$ordinal <- ordered(DF$ordinal) DF str(DF) I hope it helps. Best, Dimitris On 11/24/2010 8:27 PM, Άγνωστος μέσα στο πλήθος wrote: Hi. First of all, excu

[R] Par() Variables, Timing, and Resizing

Hello, R-help, Although this is my first post to the mailing list, I have been a subscriber for a while now, and have found all of the posts to be very informative in my quest to attempt to master R. Thank you for that in advance. This question involves the use of the par() variables pin, oma,

[R] Change the class of columns in a data frame

Hi. First of all, excuse me if I do any mistakes, but English is not a language I use very often. I have a data frame with numbers. A small part of the data frame is this: nominal ordinal 2 2 2 1 2 1 2 2 So, I want to use the gower distanc

[R] Custom ticks on x axis when dates are involved

Hi, I have a set of irregular time series and i want to produce a simple plot, with dates on x axis and attribute value of y axis. This is simple enough but my x axis is divided automatically by ticks every 5 years. I would like to have a tick every year at January 1st. I am not sure how i ca

Re: [R] Cannot update R

Two comments: http://www.r-project.org/nosvn/R.check/r-release-linux-ix86/rjags-00check.html CRAN checks show rjags is having errors (maybe because of the library issue, I'm not certain) You are using an old unstable version of 2.12.0---the stable release has come out since. It might be worth u

Re: [R] Finding the nearest data in intraday data from two zoo objects

On Wed, Nov 24, 2010 at 11:04 AM, Santosh Srinivas wrote: > Hello Group, > > I have the following options and future data in zoo objects > >> head(optData.z) >                       ExpDt OptTyp Strike TrdPrice TotTrdQty > 2009-01-01 09:55:03 20090129      1   2900 180.        50 > 2009-01-01

Re: [R] Cannot update R

Alon - It says it couldn't find it in /usr/local/lib -- is it in some non-standard location? If so, you could modify the ~/.R/Makevars file to contain PKG_LIBS=/location/of/jagslibrary before trying to build the package. - Phil On Wed, 24 Nov 2

Re: [R] Cannot update R

Thanks Phil, To complicate matters still let me add that I have been using JAGS and rjags for a while, never had a problem , so R can find JAGS and I can interface with JAGS using rjags... but I cannot update it? Your opinion? Maybe I need to update a path to rjags? maybe I can move it to the m

[R] rbind error (maybe a problem with chron package?)

Hello, I have a list of data frames that I would like to combine to a single data frame. I typically do this with: do.call(rbind, list.of.dataframes) This time, I get Error in class(x) <- cl : cannot set class to "array" unless the dimension attribute has length > 0 I'm not sure what's hap

Re: [R] Filled contour plot showing labeled isolines?

On Nov 24, 2010, at 10:19 AM, jt306 wrote: Sorry, I guess I didn't make it clear. I mentioned it in the subject of my post, but not in the body. Is it possible to have a filled contour plot (showing the colors) with the isolines labeled? Then levelplot _was_ (perhaps) the right answe

Re: [R] adding blocks to arrays

On Nov 24, 2010, at 12:08 PM, Cable, Samuel B Civ USAF AFMC AFRL/RVBXI wrote: I have a four-dimensional array. I have a three-dimensional array I want to add to it. My particular problem is conceptually very simple because the new array is just going to go into the end of the 4D array;

[R] Tinn R to R problem

I am trying to send complete R code to R from Tinn R by pressing R send:contiguous button. I get following error: > source(.trPaths[5], echo=TRUE, max.deparse.length=150) Error in source(.trPaths[5], echo = TRUE, max.deparse.length = 150) : object '.trPaths' not found Anyone has any idea? Th

Re: [R] Filled contour plot showing labeled isolines?

Sorry, I guess I didn't make it clear. I mentioned it in the subject of my post, but not in the body. Is it possible to have a filled contour plot (showing the colors) with the isolines labeled? On Wed, Nov 24, 2010 at 5:35 AM, David Winsemius [via R] < ml-node+3057276-1726133352-203...@n4.na

Re: [R] shifting down ylab in a plot

thanks Phillip for the trick. However I am a bit surprised that something as trivial as this is not accessible via par options... If somebody know how to do it with par ... or anything from graphics would be nice to know Cheers -- View this message in context: http://r.789695.n4.nabble.com/shif

[R] Aligning Grid Lines with Date Ticks in Trellis Plot

I know this issue has been discussed before, but I hoped the advent of pretty.Date would resolve it. In the following code, the first plot produces misaligned vertical grid lines, while in the second plot they are properly aligned. Is there any way to get something along the lines of the first c

Re: [R] the first. from SAS in R

Often the purpose of first/last in sas is to facilitate grouping of observations in a sequential algorithm. This purpose is better served in R by using vectorized methods like those in package plyr. Also, note that first/last has different meanings in the context of "by x;" versus "by x notsorted

Re: [R] Cannot update R

I believe R is complaining that it can't find the jags library. You can find out about it here: http://www-fis.iarc.fr/~martyn/software/jags/ There appear to be rpms available. - Phil Spector Statistical Computi

Re: [R] Population abundance, change point

I understand that smoothing splines produce continuous models, however the end product of SiZer is not a single model that is then used in any predictive manner. Rather, the end product is a map of potential changepoint locations along a gradient. Are you suggesting that SiZer would not find a

Re: [R] The nls2 function automatically prints the object!

Ivan, Try this: invisible(capture.outuput( model <- nls2 (...) )) On Wed, Nov 24, 2010 at 12:07 PM, Ivan Allaman wrote: > > Good morning gentlemen! > > When I use the function nls2, and store it in an object, that object is > automatically printed, without the summary or to draw the object. Fo

Re: [R] Running R from SAS

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of ziad.elmou...@tnsglobal.com > Sent: Wednesday, November 24, 2010 3:27 AM > To: r-help@r-project.org > Subject: [R] Running R from SAS > > Hello All, > > I am interested in runni

[R] Cannot update R

Hello I cannot update the package rjags. Although I am a su. Running openSUSE R version : language R version.string R version 2.12.0 Under development (unstable) (20

[R] diverting output from nested loops

Dear List, I have a series of nested loops with the structure shown below, and I am struggling to figure out how to divert output to folders created with dir.create() within the loops. What I need is for the output to end up as topNameK/subNameL/objNameM.pdf; what I get instead is a series of dir

Re: [R] Population abundance, change point

On Wed, 2010-11-24 at 10:59 -0500, Jonathan P Daily wrote: > I agree that SiZer is not the ultimate answer to all changepoint analysis, > but that is why there are so many changepoint detection methods used. I > will clarify, though, that my understanding of SiZer (which may be wrong) > was that

[R] adding blocks to arrays

I have a four-dimensional array. I have a three-dimensional array I want to add to it. My particular problem is conceptually very simple because the new array is just going to go into the end of the 4D array; i.e. will be added into the final index. I think I have seen an R command that can do t

Re: [R] cross validation using e1071:SVM

thank you so much for your help. if i am not wrong then createDataPartition can be used to create stratified random splits of a data set. is there other way to do that? Thank you -- View this message in context: http://r.789695.n4.nabble.com/cross-validation-using-e1071-SVM-tp3055335p3057684.

[R] negative binomial regression, unbalanced panel

I am a student who is doing empirical work for his thesis and trying to switch to R. I am familiar with Stata, and at the moment I am trying to replicate some of my previous work. I have a large unbalanced panel data set, observations for different countries between 1970 and 2007. My dependent

Re: [R] Finding the nearest data in intraday data from two zoo objects

Try data.table with the roll=TRUE argument. Set your keys and then write : futData[optData,roll=TRUE] That is fast and as you can see, short. Works on many millions and even billions of rows in R. Matthew http://datatable.r-forge.r-project.org/ "Santosh Srinivas" wrote in message news:

Re: [R] Population abundance, change point

> To: gavin.simp...@ucl.ac.uk > CC: carus...@gmail.com; marchy...@hotmail.com; r-help@r-project.org; > r-help-boun...@r-project.org > Subject: Re: [R] Population abundance, change point > From: jda...@usgs.gov > Date: Wed, 24 Nov 2010 10:59:17 -0500

Re: [R] Performance tuning tips when working with wide datasets

Richard, Try data.table. See the introduction vignette and the presentations e.g. there is a slide showing a join to 183,000,000 observations of daily stock prices in 0.002 seconds. data.table has fast rolling joins (i.e. fast last observation carried forward) too. I see you asked about that on

Re: [R] How to start default browser on R

On Mon, 2010-11-22 at 23:31 -0800, Stephen Liu wrote: > Hi David, > > I see. File: gives the full path to the .html file created/download. Then > the > browser will open that file. Thanks. > > I don't have hwriter package installed. > > A side question, what will be the corresponding comman

Re: [R] looking for the RMySQL package for R 2.12.0 under XP

Hello: That package needs to be compiled from source on the computer you are using. Instructions for doing so can be found at "http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL";. I'm not sure of all the issues, but I believe it needs to be configured to the specific version o

[R] Haver connectivity progress and issues

Dear All, I have been trying to connect to haver to download data and have made some good progress. I am stuck at a point and would appreciate your help. Haver provides a dll called DLXAPI32 using which i can connect to excel/vb/c etc. Now i tried calling it from R and used following code: o

[R] Finding the nearest data in intraday data from two zoo objects

Hello Group, I have the following options and future data in zoo objects > head(optData.z) ExpDt OptTyp Strike TrdPrice TotTrdQty 2009-01-01 09:55:03 20090129 1 2900 180.50 2009-01-01 09:55:31 20090129 1 2900 188.50 2009-01-01 09:55:37

Re: [R] Population abundance, change point

I agree that SiZer is not the ultimate answer to all changepoint analysis, but that is why there are so many changepoint detection methods used. I will clarify, though, that my understanding of SiZer (which may be wrong) was that the smoothing splines are just a vessel for finding the changepoi

Re: [R] looking for the RMySQL package for R 2.12.0 under XP

You need compile it from source code. On Wed, Nov 24, 2010 at 1:48 PM, PtitBleu wrote: > > Hello, > > I convinced our new Ph student to work with R and I managed to convince the > guy in charge of the computer networks to install R 2.12.0 (we haven't the > admin rights). > > I wanted to install

[R] Error in solve.default(V, W) with ProDenICA and the required sample size for stable result

Hi ALL, Recently I confronted with a error in using package ProDenICA, when let khttp://finzi.psych.upenn.edu/R/library/ProDenICA/html/ProDenICA.html), Error in solve.default(V, W) : 'b' must be compatible with 'a' I have examined the code of function ProDenICA(), but could not find the reason.

[R] looking for the RMySQL package for R 2.12.0 under XP

Hello, I convinced our new Ph student to work with R and I managed to convince the guy in charge of the computer networks to install R 2.12.0 (we haven't the admin rights). I wanted to install on his computer the same packages I use (but I have R 2.11) and especially RMySQL. Problem : there is

Re: [R] Conditional Operation in Matrix

Hi, Would that work for you? d <- a-b d[a==0 | b==0] <- 0 The output is as the one expected HTH, Ivan Le 11/24/2010 15:54, Romildo Martins a écrit : Hello, I would like to perform the operation a - b, only for the non-zero values. Thanks a lot! RMB a 0 0 0 0 4 0 0 0 0 0 0 0

Re: [R] problem with parameters

Well since you said that it cannot recognize the value of 'p', this suggests that you have not defined that object. Here is an example of what you might be trying to do (you can just copy and paste the code into R): ## create x matrix of 0s x <- matrix(0,

[R] Conditional Operation in Matrix

Hello, I would like to perform the operation a - b, only for the non-zero values. Thanks a lot! RMB > a 0 0 0 0 4 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 > b

Re: [R] Tobit model on unbalanced panel

For tobit regression, also see the last example of help(survreg) in the survival package. Terry Therneau __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-

[R] R encoding question

Hi, I am using RpgSQL to retrieve data from a PostgreSQL database wich is with encoding UTF8, and I have some Chinese character in one of the columns, unfortunately R can't show it correctly. > df <- dbGetQuery(con, "select * from test") > df ab 1 1 椤惧��娉\xa2 2 2 瑕�� EURO\xa1 I s

[R] nonparametric covariance analysis

Hi there, I want to do a nonparametric covariance analysis and I have tried to use the package "sm" function "sm.ancova" but it didn't work for me because I have more then one covariates (I have 18 covariates and 3 factors). I want to analyse for one factor (who has 13 levels) where the differe

Re: [R] question about constraint minimization

Hi Hans, I really appreciate your help and the detailed reply. Many thanks and have a nice Thanks giving! Best, Hao On Tue, Nov 23, 2010 at 11:46 AM, Hans W Borchers [via R] < ml-node+3055782-1866096000-202...@n4.nabble.com > wrote: > dhacademic gmail.com gmail.com> writes: > > > Hi, > > > >

Re: [R] Population abundance, change point

On Wed, 2010-11-17 at 09:17 -0500, Jonathan P Daily wrote: > Indeed I have looked into various non-standard changepoint analysis > methods. I figured the OP was more interested in traditional methods since > you have to spend less time justifying your methodology. Wavelets are one > potential no

[R] The nls2 function automatically prints the object!

Good morning gentlemen! When I use the function nls2, and store it in an object, that object is automatically printed, without the summary or to draw the object. For example. model <- nls2 (...) Number of iterations to convergence: ... Achieved convergence tolerance: ... Nonlinear regression mo

Re: [R] Need help

Thileepan Sekaran googlemail.com> writes: >I have a affymetrix data and i have performed the Normalize the > data also,Now I want to find the density plot for that data set.Pls do help > me Please (a) post this question to the Bioconductor list, where it will be more appropriate, a

Re: [R] Sub-setting a data frame by partial column names?

Hi, I haven't followed, but I think you want something like: subsetData<- bigData[ , c("Colname1", "Colname2", grep("partialName", colnames(bigData)))] HTH, Ivan Le 11/24/2010 13:51, barah a écrit : Hi Sarah, Thanks a lot for the suggestion. It is working for me. subsetData<- bigData[,grep(

  1   2   >