Re: [R] time function behavior for ts class objects

Try not reducing options(digit=)! I get (with the default value 7) time(xts) Jan Feb Mar Apr May Jun Jul Aug 1973 1973.000 1973.083 1973.167 1973.250 1973.333 1973.417 1973.500 1973.583 1974 1974.000 1974.083 1974.167 1974.250 1974.333 1974.417 1974

Re: [R] How to write loop

Alex Roy [Sun, Jun 14, 2009 at 06:43:52AM CEST]: > Hi Jiim, Thanks . I want to do the following: > > 1. each time I need to drop one column, say first column 1 from matrix X. > 2 then take out row 1 of the remainning matrix and that row becomes > response (y) > 3. do lasso regression on remain

[R] How to capture output from a subprocess?

## open a pipe a <- pipe("mplayer -slave -quiet /media/wind/Music/a.mp3","w") write("get_file_name",a) ## send get_file_name command to player. ## ANS_FILENAME='a.mp3' write("quit",a) ## exit mplayer Can I capture the output from mplayer into a R object? In this case, it is "ANS_FILENAME='a.mp3'".

Re: [R] How to write loop

Hi Jiim, Thanks . I want to do the following: 1. each time I need to drop one column, say first column 1 from matrix X. 2 then take out row 1 of the remainning matrix and that row becomes response (y) 3. do lasso regression on remaining X to y. 4. store the coefficients Similarly, in next run

Re: [R] triangular plot

Also triplot in TeachingDemos -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of David Carlson > Sent: Friday, J

Re: [R] Gtk objects disappears

I would recommend just running gtkMain(), so that GTK+ blocks R. Then you need your GUI to call gtkMainQuit() when it's time to kill R. On Fri, Jun 12, 2009 at 6:59 AM, Olivier Nuñez wrote: > Dear John, > > I have a question. > > When I run a RGtk code in my terminal (without using the R GUI) >

[R] time function behavior for ts class objects

Hi all- I am trying to use the time function for ts class objects and do not understand the return value. I want to use it to set up a time trend in arima fits. It does not seem to return a correct linear sequence that matches the underlying time series. I am running: R version 2.8.1 (2008

Re: [R] color in pairs

Try this: panel.cor <- function(x, y, ...) { r <- abs(cor(x, y)) u <- par('usr') names(u) <- c("xleft", "xright", "ybottom", "ytop") bgcolor <- ifelse(r > 0.5, 'red', 'blue') do.call(rect, c(col = bgcolor, as.list(u))) points(x, y) } pairs(USJudgeRatings, panel=panel

[R] color in pairs

Dear list Is there any way that I can color the background (not the plotted points) of each box in "pairs" ? I'd like to color them according to the correlation coefficiency of each data pair. thanks. N __ R-help@r-project.org mailing list https://s

Re: [R] External signal in ODE written in C (using deSolve and approx1?)

Dear Thomas Thank you for your excellent help. I agree method 3 would be the best one, and this is the one I have been trying to get working (like you described). So far my skills with pointers, what is returned etc in pure C code is limited, so my attempts have failed so far. If it is possible to

Re: [R] Approximation to next integer

thanks. Solved! 2009/6/14 Henrique Dallazuanna > See ?round > > On Sat, Jun 13, 2009 at 8:53 PM, Carletto Rossi wrote: > >> Hi, i have to approximate a decimal number to the next integer >> >> e.g >> >> 2.76 ---> 3 >> >> is there a function in R that make this operation?? >> >> thanks >> >>

Re: [R] Approximation to next integer

See ?round On Sat, Jun 13, 2009 at 8:53 PM, Carletto Rossi wrote: > Hi, i have to approximate a decimal number to the next integer > > e.g > > 2.76 ---> 3 > > is there a function in R that make this operation?? > > thanks > >[[alternative HTML version deleted]] > > __

[R] Approximation to next integer

Hi, i have to approximate a decimal number to the next integer e.g 2.76 ---> 3 is there a function in R that make this operation?? thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

Re: [R] Fast JSON <-> R converter?

On Fri, Jun 12, 2009 at 7:16 PM, Duncan Temple Lang wrote: > > It is not so much that rjson is implemented in R that makes it slow, just > that it does not use vectorized operations. > > > The package RJSONIO > > http://www.omegahat.org/RJSONIO > Great, I'll check it out. Thanks! kynn

Re: [R] Insignificant variable improves AIC (multinom)?

Oops. In my previous email I meant to say the following: In the AIC approach, you include a new variable or delete an existing variable when the change in the "log-likelihood" value is 2 or more. Ravi. Ravi Varadhan, Ph.D. A

Re: [R] Insignificant variable improves AIC (multinom)?

Hi Werner, AICs of nested models are compared on additive scale, not on multiplicative scale. So, you have to think about how much the AIC is decreased when you add the new variable, not the factor by which it is reduced. If you are doing a stepwise selection based on AIC, then the p-value a

Re: [R] R and Sound

Mike McGonagle wrote: Hello, I recently have started to learn about R, and am interested in using it to generate data that would then get used to produce audio. Basically, I am "envisioning" using R to create the data, based on probability functions (and fractal data), which then gets output

Re: [R] BLAST search from R

Nigel Birney wrote: > > Hi, > > I have FASTA sequences and would like to do BLAST search of all against > all. Can I do this in R? > > thanks, > Nigel > I would suggest asking this question on the Bioconductor list: http://www.bioconductor.org/docs/mailList.html -- View this message in

[R] BLAST search from R

Hi, I have FASTA sequences and would like to do BLAST search of all against all. Can I do this in R? thanks, Nigel -- View this message in context: http://www.nabble.com/BLAST-search-from-R-tp24011761p24011761.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Insignificant variable improves AIC (multinom)?

> >Hi, > > > >I am trying to specify a multinomial logit model using the multinom function > from the nnet package. Now I add another independent variable and it halves > the > AIC as given by summary(multinom()). But when I call Anova(multinom()) from > the > car package, it tells me that t

Re: [R] coupled ODE population model

Hi, Justin: In addition to the example below, Thomas Petzoldt also has two vignettes included with his "simecol" package. These might help you more easily use his package. Have you looked for what you want using the "RSiteSearch.function" in the "RSiteSearch" package? Conside

[R] How to write loop

Dear all, I want to do the following process as a loop ( to run automatically with dimension of X, here 50). How can I do that? Your cooments will be highly appreciable. Alex *# Code:* library(lars) library(chemometrics) X<-matrix(rnorm(2500),ncol=50) dim(X) # [1] 50 50 X1<-X[,

Re: [R] Order data frame

Dear Marcio, Take a look at http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:sort HTH, Jorge On Fri, Jun 12, 2009 at 2:54 PM, MarcioRibeiro wrote: > > Hi listers, > I have a data frame with 5 variables... And I wan

Re: [R] Order data frame

I asked Google for "order data.frame" and got some good answers Uwe Ligges MarcioRibeiro wrote: Hi listers, I have a data frame with 5 variables... And I want reorder the data frame according a specific variable X... I checked the sort and order functions... In my case I think that the OR

Re: [R] Hmisc summarize() with level "" in by variable

Sorry about the bug, which is now fixed. You can get the fix by entering source('http://biostat.mc.vanderbilt.edu/cgi-bin/viewvc.cgi/*checkout*/Hmisc/trunk/R/summary.formula.s?rev=661') until we update the package. Frank Michael Erickson wrote: I was using summarize() in a data set in which

Re: [R] Programmatically copying a graphic to the clipboard

Possible using RGtk2. Just need to get it into a GdkPixbuf. See gtkClipboardSetImage(). To get a graphic as a pixbuf, you can use cairoDevice (drawing to a GdkPixmap and copying over), or output the graphic to a temporary file and read it back in. Michael On Fri, Jun 12, 2009 at 5:42 AM, Hadley W

[R] Fitting Mixture of Non-Central Student's t Distributions

Dear all, I am attempting to model some one-dimensional data using a mixture model of non-central Student's t distributions. However, I haven't been able to find any R package that provides this functionality. Could there be a way to "manipulate" the EM algorithms from the mixdist or mixtools pac

Re: [R] formatting numbers along axes as percents for perspective plot

Further option would be scatterplot3d in {scatterplot3d} which does provide facility for sending a vector of axis labels to the x,y, and z axes. David On Jun 12, 2009, at 10:48 PM, Stephen Samaha wrote: Hi, that's what I did for the time being. Ideally, I'd like to have the % sign after t

Re: [R] 'Errors' with Ubuntu

Unbuntu != R ==> wrong list Uwe Ligges Len Vir wrote: Good Day To You All, I have installed Ubuntu 9.04 on a Laptop running Vista Business. When I try to start my little machine, GRUB keeps telling me there is an error: Either: Error 25, Error 21, or Error 5. Once in a while, it is 'E

[R] spatstat: help with calculating number of clusters, cluster size

Hi, I am new to spatial statistics and have a few basic queries. I am analysing the spatial pattern of trees on a large plot. I wanted to know if there is there some way of obtaining estimates of the number of clusters and cluster-radius (in an iterative fashion) after fitting a spatial point proc

Re: [R] External signal in ODE written in C (using deSolve and approx1?)

Glenn Woodart wrote: Dear list The deSolve package allows you to specify the model code in C or Fortran. Thanks to the excellent vignette this works fine. However I have not yet managed to use forcing functions in C code. In pure R code this works very well with approxfun() specified outside th

Re: [R] coupled ODE population model

Justin Frank wrote: I'm fairly new to R, and I'm trying to write out a population model that satisfies the following; the system consists of s species, i= 1, 2,...,s network of interactions between species is specified by a (s x s) real matrix, C[i,j] x[i] being the relative population of the

Re: [R] Insignificant variable improves AIC (multinom)?

Werner Wernersen wrote >Hi, > >I am trying to specify a multinomial logit model using the multinom function >from the nnet package. Now I add another independent variable and it halves >the AIC as given by summary(multinom()). But when I call Anova(multinom()) >from the car package, it tells me

[R] Griddy Gibbs

Hello, I am trying to code "Griddy gibbs" sampling just with the altorithm which were considered by Ritter and Tanner(1992) and R to implement Bayesian GARCH with griddy gibbs. It is a project for final exam of Bayesian statistics. If you have source code for "Griddy Gibbs"

[R] Hmisc summarize() with level "" in by variable

I was using summarize() in a data set in which one of the levels of the by variable was "". The summary statistic was consistently off by one level and the "" level was not in the output data frame. I tried to report it as a bug, but I could not log into the Hmisc bug reporting website to do so.

[R] 'Errors' with Ubuntu

Good Day To You All, I have installed Ubuntu 9.04 on a Laptop running Vista Business. When I try to start my little machine, GRUB keeps telling me there is an error: Either: Error 25, Error 21, or Error 5. Once in a while, it is 'Error 22'. And then there is no way to go any further: End of th