Re: [R] Alternate to for-loop

Stefan Evert wrote: A couple of remarks on vQ's naive benchmark: f.rep = function(n, m) replicate(n, rnorm(m)) I suppose you meant f.rep = function(n, m) replicate(n, mean(rnorm(m))) which doesn't make a substantial speed difference, though. indeed, thanks; i've already posted a

[R] stopping stepAIC (MASS package)

Dear R-listers, I am building a prediction model starting with many, many variables. I use the 'stepAIC' procedure in the MASS package, and the model building process takes several hours to complete. At the very end, I am occasionally confronted with warnings like this one: Warning messages: 1:

[R] problem with dates in zoo package

i have the following code - assimilating the maximum annual discharge each year ffrom a daily discharge record from year 1989-2005. m - read.table(D:/documents/5 stations/01014000.csv, sep =,) z - zoo(m[,4],as.Date(as.character(m[,3]), %m/%d/%Y)) x - aggregate(z,

[R] Efficient matrix computations

Hi, I am looking for two ways to speed up my computations: 1. Is there a function that efficiently computes the 'sandwich product' of three matrices, say, ZPZ' 2. Is there a function that efficiently computes the determinant of a positive definite symmetric matrix? Thanks, S.A.

[R] difficulty in starting time sequence from 00:00:00

Dear R-experts, Need your help. Dear R- Experts, Seek your help. I created a time sequence using: x[i] -chron(dates, tt, format=c(dates=y-m-d, tt=h:m:s)) first element in the list is displayed as: (09-01-01 00:00:00) Further elements are: (09-01-01 00:01:00) (09-01-01 00:02:00) (09-01-01

Re: [R] How to add direction of time to plot.circular()

Michael Kubovy wrote: Dear r-helpers, I want to show that time is flowing CCW in the following: require(circular) len - 8 labl - as.character(c(0, 1, 1, 1, 0, 0, 1, 0)) r - circular(2*pi* (rep(c(1, 3, 6), each = 200)/len + rnorm(600, 0, 0.025))) r.dens - density(r, bw = 25, adjust = 4,

Re: [R] Create package with Fortran 90 and C code

On Tue, 17 Feb 2009, Nathan S. Watson-Haigh wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 I'm trying to add some Fortran 90 code to an existing package. When I compile and load the file manually like: SHELL R CMD SHLIB file.f90 R dyn.load(file.so) I can use the .Fortran() fine.

[R] Subset Regression Package

Dear all , Is there any subset regression (subset selection regression) package in R other than leaps? Thanks and regards Alex [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Efficient matrix computations

it would be a bit more helpful if we knew more info regarding these matrices, for instance is P diagonal, etc. In any case, you could have a look at crossprod() # and tcorssprod() and, for the determinant maybe prod(eigen(mat, symmetric = TRUE, only.values = FALSE)$values) # or

[R] Ingore this message. I have already solved the problem..

Ingore this message. I have already solved the problem.. Regards, Suresh Suresh_FSFM wrote: Dear R-experts, Need your help. Dear R- Experts, Seek your help. I created a time sequence using: x[i] -chron(dates, tt, format=c(dates=y-m-d, tt=h:m:s)) first element in the list

Re: [R] Efficient matrix computations

sorry, in my previous e-mail it should be tcrossprod() # and prod(eigen(mat, symmetric = TRUE, only.values = TRUE)$values) Best, Dimitris Shimrit Abraham wrote: Hi, I am looking for two ways to speed up my computations: 1. Is there a function that efficiently computes the 'sandwich

Re: [R] Alternate to for-loop

I ran into a similar issue with a simple benchmark the other day, where a plain loop in Lua was faster than vectorised code in R ... hmm, would you be saying that r's vectorised performance is overhyped? or is it just that non-vectorised code in r is slow? What I meant, I guess, was

Re: [R] Efficient matrix computations

Thanks for your suggestions. I'll try to implement what you suggested. Perhaps the following information can help you to think of alternative ways to speed up computations: I am coding the Kalman Filter in R because I have certain requirements that are not provided by the packages regarding

Re: [R] scatterplot and correlation for weird data format

William Simpson wrote: I have data in a format like this: namessexsex viewnum rating rt ahl4f m f 56 -1082246 ahl4f m f 74 85 1444 ahl4f m f 52 151 1595 ahl4f m f

Re: [R] Alternate to for-loop

Stefan Evert wrote: hmm, would you be saying that r's vectorised performance is overhyped? or is it just that non-vectorised code in r is slow? What I meant, I guess, was (apart from a little bit of trolling) that I'd had misconceptions about the speed differences between loops and

Re: [R] problem with dates in zoo package

See ?which.max : library(zoo) z - read.zoo(myfile.dat, sep = ,, format = %m/%d/%Y) f - function(x) time(x)[which.max(x)] ix - tapply(z, floor(as.yearmon(time(z))), f) z[ix] On Tue, Feb 17, 2009 at 3:58 AM, CJ Rubio cjru...@kongju.ac.kr wrote: i have the following code - assimilating the

[R] Combination

Hello, I have a sequence of numbers: seq(1:50) and I would like to have all the possible combinations with this numbers without repeating any combination: 11, 12, 13, ... ,22,23,24,... How can I do it? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de

Re: [R] Combination

Try this: apply(expand.grid(1:50, 1:50), 1, paste, collapse = '') On Tue, Feb 17, 2009 at 8:37 AM, Dani Valverde daniel.valve...@uab.catwrote: Hello, I have a sequence of numbers: seq(1:50) and I would like to have all the possible combinations with this numbers without repeating any

Re: [R] Combination

Hi Dani, see ?combn . combn(1:50,2) gives you all combinations as matrix. you can do sth like apply(combn(1:50,2),2, paste, sep=, collapse=) to get concatenated results. hth. Dani Valverde schrieb: Hello, I have a sequence of numbers: seq(1:50) and I would like to have all the possible

[R] Problem with calculating the residuals in an AR(p)-model

Hi, I'm trying to calculate the residuals in a one-sided AR(p) model: sum (phi_j (X_t-j - EX)) = epsilon My X is an ARMA(1,1)-model, which can be represented as an AR(infty) model. I calculate the order of the model with AIC and the parameters with the Yule-Walker method. For the residuals, I

Re: [R] Using eval in multinom argument

Thanks for your help, I managed to get it working like this: fml-as.formula(paste(y ~, paste(PCnames, collapse=+))) y = class.ind(grp[outgroup]) z1=multinom(formula = fml, data=data.frame(scores)) Daniel Crouch **QUOTE: Forget eval(parse(text = )) See ?as.formula

Re: [R] Comparison of age categories using contrasts

Hi Dylan, Am I trying to use contrast.Design() for something that it was not intended for? ... I think Prof. Harrell's main point had to do with how interactions are handled. You can also get the kind of contrasts that Patrick was interested in via multcomp. If we do this using your

Re: [R] Comparison of age categories using contrasts

On 2/16/2009 10:18 PM, Dylan Beaudette wrote: On Mon, Feb 16, 2009 at 5:28 PM, Patrick Giraudoux patrick.giraud...@univ-fcomte.fr wrote: Greg Snow a écrit : One approach is to create your own contrasts matrix: mycmat - diag(8) mycmat[ row(mycmat) == col(mycmat) + 1 ] - -1 mycmati -

Re: [R] assuming AR(1) residuals in OLS

Thank you for the lightening replies. I tested various corStruct objects (?corClasses) using the nlme package and all work flawlessly. My best regards to all... Constantine Tsardounis __ R-help@r-project.org mailing list

Re: [R] Combination

Eik Vettorazzi wrote: Hi Dani, see ?combn . combn(1:50,2) gives you all combinations as matrix. you can do sth like apply(combn(1:50,2),2, paste, sep=, collapse=) note that you could simplify the above by combn(50, 2, paste, collapse = ) Best, Dimitris to get concatenated results.

Re: [R] Comparison of age categories using contrasts

Hi Dylan, Chuck, contrast(l, a=list(f=levels(d$f)[1:3], x=0), b=list(f=levels(d$f)[4], x=0)) There is a subtlety here that needs to be emphasized. Setting the interacting variable (x) to zero is reasonable in this case, because the mean value of rnorm(n) is zero. However, in the real world

[R] Chromatogram deconvolution and peak matching

Hi, I'm trying to match peaks between chromatographic runs. I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my

Re: [R] stopping stepAIC (MASS package)

Peter Jepsen PJ at DCE.AU.DK writes: [snip] ... I would save valuable time if the procedure would stop automatically when it encounters any anomaly that causes it to warn me. Is this possible? And if so, how? Maybe options(warn=2) will do what you want. Ben Bolker

Re: [R] assuming AR(1) residuals in OLS

Thank you Gabor Grothendieck for your message. I would surely like to say, that if someone wants to assume AR(1) residuals, running the regression y ~ x, could run gls(y~x, correlation = corAR1(0, ~1)) Constantine Tsardounis http://www.costis.name

[R] joining one-to-many

Hello list, I am wondering if a joining one-to-many can be done a little bit easier. I tried merge function but I was not able to do it, so I end up using for and if. Suppose you have a table with locations, each location repeated several times, and some attributes at that location. The

Re: [R] joining one-to-many

On Tue, Feb 17, 2009 at 8:33 AM, Monica Pisica pisican...@hotmail.com wrote: Hello list, I am wondering if a joining one-to-many can be done a little bit easier. I tried merge function but I was not able to do it, so I end up using for and if. Suppose you have a table with locations,

Re: [R] joining one-to-many

on 02/17/2009 08:33 AM Monica Pisica wrote: Hello list, I am wondering if a joining one-to-many can be done a little bit easier. I tried merge function but I was not able to do it, so I end up using for and if. Suppose you have a table with locations, each location repeated several times,

[R] frequency table for multiple variables

Hi r-help! Consider the following data-frame: var1 var2 var3 1 314 2 223 3 223 4 44 NA 5 435 6 223 7 343 How can I get R to convert this into the following? Value 1 2 3 4 5 var1 0 3 2 2 0

Re: [R] joining one-to-many

Try merge(t1, t2) On Tue, Feb 17, 2009 at 9:33 AM, Monica Pisica pisican...@hotmail.com wrote: Hello list, I am wondering if a joining one-to-many can be done a little bit easier. I tried merge function but I was not able to do it, so I end up using for and if. Suppose you have a table

[R] spss-file problem with foreign 0.8-32

Hi, after updating to foreign version 0.8-32, I experienced the following error when I tried to load a SPSS file: Fehler in inherits(x, factor) : objekt cp nicht gefunden Zusätzlich: Warning message: In read.spss(***l.sav, use.value.labels = TRUE, to.data.frame = TRUE) : ***.sav:

[R] Uninstall question

I need to uninstall R 2.7.1 from my Mac. What is the best way to uninstall it? Simply delete the R icon in the Applications folder? Or is it more involved? TIA, Anjan -- = anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine

Re: [R] joining one-to-many

Hi Monica, merge(t1, t2) works on your example. So why don't you use merge? HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie,

Re: [R] Uninstall question

This is on the Mac FAQ: http://cran.cnr.berkeley.edu/bin/macosx/RMacOSX-FAQ.html#How-can-R-for-Mac-OS-X-be-uninstalled_003f HTH, --sundar On Tue, Feb 17, 2009 at 7:17 AM, ANJAN PURKAYASTHA anjan.purkayas...@gmail.com wrote: I need to uninstall R 2.7.1 from my Mac. What is the best way to

Re: [R] joining one-to-many

Ok, I feel properly ashamed. I suppose my real data is a little bit different than my toy data (although i don't know how) because i did try the merge function as simple as merge(t1, t2) and did not work. Maybe a reset of my session will solve my problems and more coffee my confusion.

[R] Survival-Analysis: How to get numerical values from survfit (and not just a plot)?

Hi! I came across R just a few days ago since I was looking for a toolbox for cox-regression. I´ve read Cox Proportional-Hazards Regression for Survival Data Appendix to An R and S-PLUS Companion to Applied Regression from John Fox. As described therein plotting survival-functions works well

Re: [R] Survival-Analysis: How to get numerical values from survfit (and not just a plot)?

Hi, See ?survfit.object if fit is the object you get using survfit, fit$surv will give you the survival probability. Best, arthur Bernhard Reinhardt wrote: Hi! I came across R just a few days ago since I was looking for a toolbox for cox-regression. I´ve read Cox Proportional-Hazards

Re: [R] spss-file problem with foreign 0.8-32

Harry Haupt wrote: Hi, after updating to foreign version 0.8-32, I experienced the following error when I tried to load a SPSS file: Fehler in inherits(x, factor) : objekt cp nicht gefunden Zusätzlich: Warning message: In read.spss(***l.sav, use.value.labels = TRUE, to.data.frame = TRUE) :

[R] R scripts and parameters

A couple of weeks ago I asked how it is possible to run an R script (not a function) passing some parameters. Someone suggested the function commandArgs(). I read the on-line help and found no clarifying example. Therefore I do not know how to use it appropriately. I noticed this function

Re: [R] frequency table for multiple variables

on 02/17/2009 09:06 AM Hans Ekbrand wrote: Hi r-help! Consider the following data-frame: var1 var2 var3 1 314 2 223 3 223 4 44 NA 5 435 6 223 7 343 How can I get R to convert this into the

Re: [R] R crash after fGarch update

Prof Ripley: Many thanks - it did indeed say it cannot find fGarch after I tried your advice - but a completely clean re-install did the trick. John On Tue, Feb 17, 2009 at 1:09 AM, Prof Brian Ripley rip...@stats.ox.ac.ukwrote: Start R with --vanilla, or rename youe saved workspace (.RData).

[R] word wrap using Hershey fonts

Hi, is it possible to wrap a text using Hershey fonts? \n does not work! Thanks in advance, Martina __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Survival-Analysis: How to get numerical values from survfit (and not just a plot)?

I seriously doubt that a survfit object could only contain that information. I suspect that you are erroneously thinking that what print.survfit offers is the entire story. What does str(survfit(formula, data=dataframe) ) show you? data(aml) aml.mdl - survfit(Surv(time, status) ~ x,

Re: [R] R scripts and parameters

Try this: A - 1 B - 2 C - 3 source(myfile.R) Now the code in myfile can access A, B and C. On Tue, Feb 17, 2009 at 10:55 AM, mau...@alice.it wrote: A couple of weeks ago I asked how it is possible to run an R script (not a function) passing some parameters. Someone suggested the function

Re: [R] spss-file problem with foreign 0.8-32

Peter Dalgaard wrote: Yes, something in the logic appears to have gotten garbled. It's in this part of read,spss: if (is.character(reencode)) { cp - reencode reencode - TRUE } else if (codepage = 500 || codepage = 2000) { attr(rval, codepage) -

Re: [R] R scripts and parameters

On 2/17/2009 10:55 AM, mau...@alice.it wrote: A couple of weeks ago I asked how it is possible to run an R script (not a function) passing some parameters. Someone suggested the function commandArgs(). I read the on-line help and found no clarifying example. Therefore I do not know how to use

[R] ARIMA models

is there some sort of R function which can advise me of the best ARIMA(p,q,r) model to use based on the Schwarz criterion e.g for e.g p=0-5, q =0, r=0-5 or for example p+r 5??? or is this something I will have to write my own code for? Thanks Emma -- View this message in context:

[R] How to connect R and WinBUGS/OpenBUGS/LinBUGS in Linux in Feb. 2009

Hi all, I've managed to get JAGS working on my Ubuntu Hardy Linux with a 32-bit computer and AMD processors using R 2.8.1. JAGS is great. I've read that JAGS is the fastest, but that hasn't been my experience. At any rate, I have more experience with WinBUGS under Windows and would like a

Re: [R] Survival-Analysis: How to get numerical values from survfit (and not just a plot)?

Hi Bernhard, I'm wondering what you will expect to get in dividing two proportional survival curves from a fitted cox model. Anyway, you can provide a newdata object to the survfit function containing any combination of cofactors you are interested in and then use summary, eg: fit - coxph(

Re: [R] Comparison of age categories using contrasts

Hi Dylan, Chuck, Mark Difford wrote: Coming to your question [?] about how to generate the kind of contrasts that Patrick wanted using contrast.Design. Well, it is not that straightforward, though I may have missed something in the documentation to the function. In the past I have

Re: [R] Chromatogram deconvolution and peak matching

Hoi Bart, I think you're right that ALS should be applicable to this problem. Unfortunately in writing I see that there is a bug when the spectra are NOT constrained to nonnegative values (the package has been used to my knowledge only in fitting multiway mass spectra thus far, where this

[R] creating a map

I'm trying to create a fairly basic map using R. What i want to get is the map of the country with circles representing a count of students in each state. What I've done so far is as following - map(state) symbols(data1$count,circles=log(data1$count)*3,fg=col,bg=col,add=T,inches=F) this gives

Re: [R] ARIMA models

see auto.arima in the forecast package. On Tue, Feb 17, 2009 at 10:20 AM, emj83 stp08...@shef.ac.uk wrote: is there some sort of R function which can advise me of the best ARIMA(p,q,r) model to use based on the Schwarz criterion e.g for e.g p=0-5, q =0, r=0-5 or for example p+r 5??? or is

Re: [R] How to connect R and WinBUGS/OpenBUGS/LinBUGS in Linux in Feb. 2009

Paul Heinrich Dietrich wrote: Hi all, I've managed to get JAGS working on my Ubuntu Hardy Linux with a 32-bit computer and AMD processors using R 2.8.1. JAGS is great. I've read that JAGS is the fastest, but that hasn't been my experience. At any rate, I have more experience with WinBUGS

[R] Percentiles/Quantiles with Weighting

Hi All, I am looking at applications of percentiles to time sequenced data. I had just been using the quantile function to get percentiles over various periods, but am more interested in if there is an accepted (and/or R-implemented) method to apply weighting to the data so as to weigh recent

Re: [R] creating a map

Two places that have worked examples leap to mind: --- Sarkar's online accompaniment to his book: http://lmdvr.r-forge.r-project.org/figures/figures.html Thumbing through the hard copy I see Figure 6.5 might of interest. --- Addicted to R's graphics gallery:

[R] merging files with different structures

Hello list, Thanks in advance for any help. I have many (approx 20) files that I have merged. For example d1-read.csv(AlleleReport.csv) d2-read.csv(AlleleReport.csv) m1 - merge(d1, d2, by = c(IND, intersect(colnames(d1), colnames(d2))), all = TRUE) m2 - merge(m1, d3, by = c(IND,

Re: [R] spss-file problem with foreign 0.8-32

It is a issue specific to 0.8-32 and some files (most likely those with some (not all) Windows codepages declared). We are trying to collect together some examples, and will update foreign accordingly later in the week. On Tue, 17 Feb 2009, Harry Haupt wrote: Hi, after updating to foreign

Re: [R] merging files with different structures

Hello list, I am sorry for the previous half post. I accidentily hit send. Thanks again in advance for any help. I have many (approx 20) files that I have merged. Each data set contains rows for individuals and data in 2 - 5 columns (depending upon which data set). The individuals in

Re: [R] Percentiles/Quantiles with Weighting

I do know that Harrell's Quantile function in the Hmisc package will allow quantile estimates from models. Whether it is general enough to extend to time series, I have no experience and cannot say. -- David Winsemius On Feb 17, 2009, at 11:57 AM, Brigid Mooney wrote: Hi All, I am

Re: [R] creating a map

You need to give the symbols function the locations where you want the centers of the circles to be. Some datesets with map information also have centers of the states that you can use, for the USA, there is the state.center dataset that may work for you, or the maptools package function

Re: [R] Overdispersion with binomial distribution

Jessica L Hite/hitejl/O/VCU hitejl at vcu.edu writes: I am attempting to run a glm with a binomial model to analyze proportion data. I have been following Crawley's book closely and am wondering if there is an accepted standard for how much is too much overdispersion? (e.g. change in AIC has

Re: [R] creating a map

Thanks Greg, do you know where i can find the sate.center dataset that you mention? On Tue, Feb 17, 2009 at 12:28 PM, Greg Snow greg.s...@imail.org wrote: You need to give the symbols function the locations where you want the centers of the circles to be. Some datesets with map information

[R] Error with make with R-devel

Hi, I am getting an error compiling the R-devel on a suse architecture 64-bit architecture. The cp attribute is sending 'trusted.lov' and an error. This is a sample of the output: make[3]: Entering directory `/lustre/people/schaffer/R-devel/src/library/base' building package 'base' make[4]:

[R] printing out the summary for lm into a txt file

Hi All, I am trying to run several linear regressions and print out the summay and the anova reslts on the top of each other for each model. Below is a sample progarm that did not work. is it possible to print the anova below the summary of lm in one file? thanks for your help

Re: [R] how to control the overall shape of a figure?

Thanks very much, exactly what I need. Oliver On Feb 16, 10:36 pm, Marc Schwartz marc_schwa...@comcast.net wrote: on 02/16/2009 07:51 PM Oliver wrote: hi, I am a R beginner. One thing I notice is that when do graphing is, if I want to draw two figures in a row such as this:

[R] multiple levels of nesting with trellis plots

Hi guys, I have a tricky problem that I'd appreciate your help with. I have two categorical variables, say varA and varB and an associated frequency Freq for combinations of the levels of varA and varB. This was created with a table() call. I'd now like to make panel plots of the frequency. I

Re: [R] Overdispersion with binomial distribution

On Tue, 17 Feb 2009, Ben Bolker wrote: Jessica L Hite/hitejl/O/VCU hitejl at vcu.edu writes: I am attempting to run a glm with a binomial model to analyze proportion data. I have been following Crawley's book closely and am wondering if there is an accepted standard for how much is too much

Re: [R] printing out the summary for lm into a txt file

did not work. Might that mean errors? Care to share? Running that through my wetware R interpreter, I think I am seeing you ask for creation of models with variable outcomes from a the first 100 columns of data. And then you are specifying a formula that includes a weird mixture of 3

Re: [R] Overdispersion with binomial distribution

Thanks for the clarification. I actually had MASS open to that page while I was composing my reply but forgot to mention it (trying to do too many things at once) ... Ben Bolker Prof Brian Ripley wrote: On Tue, 17 Feb 2009, Ben Bolker wrote: Jessica L Hite/hitejl/O/VCU hitejl at vcu.edu

Re: [R] creating a map

It should be in the datasets package that is automatically loaded with R (at least my copy), try ?state and you should see the help for it and others. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 From: Alina Sheyman

Re: [R] printing out the summary for lm into a txt file

The sink() command stops the sinking, so you send the lm output to the file, then stop the sinking before printing out the anova result. So the simplest thing to try is to put the first sink (with the filename and append=T) before you start the loop, remove all calls to sink within the loop,

Re: [R] Percentiles/Quantiles with Weighting

Thanks for pointing me to the quantreg package as a resource. I was hoping to ask be able to address one quick follow-up question... I get slightly different variants between using the rq funciton with formula = mydata ~ 1 as I would if I ran the same data using the quantile function. Example:

Re: [R] Percentiles/Quantiles with Weighting

url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champaign, IL 61820 On Feb 17, 2009, at 1:58 PM, Brigid Mooney wrote:

Re: [R] printing out the summary for lm into a txt file

Thanks a lot for your help, it worked. Greg Snow-2 wrote: The sink() command stops the sinking, so you send the lm output to the file, then stop the sinking before printing out the anova result. So the simplest thing to try is to put the first sink (with the filename and append=T)

Re: [R] printing out the summary for lm into a txt file

this is the error message that I was getting In sink() ... : no sink to remove i got it to work . thanks for the help David Winsemius wrote: did not work. Might that mean errors? Care to share? Running that through my wetware R interpreter, I think I am seeing you ask for

Re: [R] Subset Regression Package

Alex Roy wrote: Dear all , Is there any subset regression (subset selection regression) package in R other than leaps? RSiteSearch({subset regression}) doesn't turn up much other than special-purpose tools for ARMA models etc.. What does leaps not do that you need

Re: [R] Whitening Time Series

Hi Bob - your suggesting worked out great... Many thanks! Also, thanks everyone for the other suggestions! Bob McCall wrote: Look in the package forecast for the function Arima. It will do what you want. It's different than arima function in the stats package. Bob Pele wrote: Hi R

[R] Cross classified or Multiple membership or Hierarchical (3 level ) logistic models using Umacs

Dear R users, I would like to fit cross classified or multiple membership logistic models or a 3 level hierarchical logistic model using the Umacs package. Can anyone advise me on how to proceed or better point me to examples of how its done. Regards, -- Luwis Diya, Leuven Biostatistics and

[R] using sapply to apply function to some columns of a dataframe

Hello: I would like to sum every x columns of a dataframe for each row. For instance, if x is 10, then for dataframe df, this function will sum the first ten elements together and then the next ten: sapply(list(colnames(df)[1:10], colnames(df)[11:20]),function(x)apply( df[,x], 1, sum)) If the

Re: [R] Subset Regression Package

Alex Roy alexroy2008 at gmail.com writes: Dear all , Is there any subset regression (subset selection regression) package in R other than leaps? Lars and Lasso are other 'subset selection' methods, see the corresponding packages 'lars' and 'lasso2' and its description in

[R] cumsum vs. sum

I recently traced a bug of mine to the fact that cumsum(s)[length(s)] is not always exactly equal to sum(s). For example, x-1/(12:14) sum(x) - cumsum(x)[3] = 2.8e-17 Floating-point addition is of course not exact, and in particular is not associative, so there are various possible

Re: [R] using sapply to apply function to some columns of a dataframe

On 17/02/2009 4:42 PM, mwestp...@worldbank.org wrote: Hello: I would like to sum every x columns of a dataframe for each row. For instance, if x is 10, then for dataframe df, this function will sum the first ten elements together and then the next ten: sapply(list(colnames(df)[1:10],

[R] matrix output

Hi friends, I have questions about printing a pretty big size matrix. As you could see from below, the matrix wasn't showed in R at full size (11X11), but it was cut partly into three smaller matrices (11X4,11X4,11X3). I'm wondering if there is a way to show the whole matrix with dimension

Re: [R] matrix output

On 17/02/2009 5:31 PM, phoebe kong wrote: Hi friends, I have questions about printing a pretty big size matrix. As you could see from below, the matrix wasn't showed in R at full size (11X11), but it was cut partly into three smaller matrices (11X4,11X4,11X3). I'm wondering if there is a way

[R] Python and R

Hello all, I am just wondering if any of you are doing most of your scripting with Python instead of R's programming language and then calling the relevant R functions as needed? And if so, what is your experience with this and what sort of software/library do you use in combination with

Re: [R] Percentiles/Quantiles with Weighting

Here is one kind of weighted quantile function. The basic idea is very simple: wquantile - function( v, w, p ) { v - v[order(v)] w - w[order(v)] v [ which.max( cumsum(w) / sum(w) = p ) ] } With some more error-checking and general clean-up, it looks like this: # Simple weighted

Re: [R] frequency table for multiple variables

On Tue, Feb 17, 2009 at 10:00:40AM -0600, Marc Schwartz wrote: on 02/17/2009 09:06 AM Hans Ekbrand wrote: Hi r-help! Consider the following data-frame: var1 var2 var3 1 314 2 223 3 223 4 44 NA 5 435 6 2

Re: [R] Python and R

Esmail Bonakdarian wrote: I am just wondering if any of you are doing most of your scripting with Python instead of R's programming language and then calling the relevant R functions as needed? No, but if I wanted to do such a thing, I'd look at Sage: http://sagemath.org/ It'll give you

Re: [R] Python and R

2009/2/17 Esmail Bonakdarian esmail...@gmail.com: Hello all, I am just wondering if any of you are doing most of your scripting with Python instead of R's programming language and then calling the relevant R functions as needed? I tend to use R in its native form for data analysis and

Re: [R] Python and R

Hello! On Tue, Feb 17, 2009 at 5:58 PM, Warren Young war...@etr-usa.com wrote: Esmail Bonakdarian wrote: I am just wondering if any of you are doing most of your scripting with Python instead of R's programming language and then calling the relevant R functions as needed? No, but if I

[R] Saving data to a MS Access Database

I have the following dataframe: ad - data.frame(dates, av, sn$SectorName) colnames(ad) - c(Date, Value, Tag) which has data (rows 10 to 20, for example) as follows: DateValue Tag 10 2008-01-16-0.20875ConDisc

[R] ylim in plot.ts?

Hello, I tried to use ylim=c(x,y) in a plot.ts(). This was ignored. How can I achieve it to create such graphics? Ciao, Oliver Bandel __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] :How to insert a column in a data frame

Hi dear list, I wonder if somebody can help me with this. I have a text file with 300 rows and around 30 columns and I need to insert a column that has the number 1 in every row. This new column should be placed between columns 6 and 7. As an example: I would want to insert a column

Re: [R] Python and R

On Tue, Feb 17, 2009 at 6:05 PM, Barry Rowlingson b.rowling...@lancaster.ac.uk wrote: 2009/2/17 Esmail Bonakdarian esmail...@gmail.com: When I need to use the two together, it's easiest with 'rpy'. This lets you call R functions from python, so you can do: from rpy import r r.hist(z)

[R] Normal cdf modified function

I wonder if an R package would have a function that calculates the following. Let Y be a normal multivariate function. For example, let Y have 4 dimensions. I want to calculate P(Y1 Z1, Y2 Z2, Y3 Z3, Y4 Z4). There are R functions to do the calculation if all the inequalities are of the type

  1   2   >