[R] 'Errors' with Ubuntu

2009-06-13 Thread Len Vir
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

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

2009-06-13 Thread Michael Erickson
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. I

[R] Griddy Gibbs

2009-06-13 Thread 최현희
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,

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

2009-06-13 Thread Peter Flom
Werner Wernersen pensterfuz...@yahoo.de 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

Re: [R] coupled ODE population model

2009-06-13 Thread Thomas Petzoldt
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] External signal in ODE written in C (using deSolve and approx1?)

2009-06-13 Thread Thomas Petzoldt
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

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

2009-06-13 Thread Soumya Prasad
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

Re: [R] 'Errors' with Ubuntu

2009-06-13 Thread Uwe Ligges
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

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

2009-06-13 Thread David Winsemius
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

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

2009-06-13 Thread Balzer Susanne
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

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

2009-06-13 Thread Michael Lawrence
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

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

2009-06-13 Thread Frank E Harrell Jr
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

Re: [R] Order data frame

2009-06-13 Thread Uwe Ligges
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

Re: [R] Order data frame

2009-06-13 Thread Jorge Ivan Velez
Dear Marcio, Take a look at http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:sort 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 mes...@pop.com.br wrote: Hi listers, I have a data frame with 5

[R] How to write loop

2009-06-13 Thread Alex Roy
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

Re: [R] coupled ODE population model

2009-06-13 Thread spencerg
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? Consider,

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

2009-06-13 Thread Werner Wernersen
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 this added

[R] BLAST search from R

2009-06-13 Thread Nigel Birney
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] BLAST search from R

2009-06-13 Thread Ben Bolker
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 context:

Re: [R] R and Sound

2009-06-13 Thread Uwe Ligges
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] Insignificant variable improves AIC (multinom)?

2009-06-13 Thread Ravi Varadhan
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

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

2009-06-13 Thread Ravi Varadhan
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.

Re: [R] Fast JSON - R converter?

2009-06-13 Thread Kynn Jones
On Fri, Jun 12, 2009 at 7:16 PM, Duncan Temple Lang dun...@wald.ucdavis.edu 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.

[R] Approximation to next integer

2009-06-13 Thread Carletto Rossi
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

Re: [R] Approximation to next integer

2009-06-13 Thread Henrique Dallazuanna
See ?round On Sat, Jun 13, 2009 at 8:53 PM, Carletto Rossi nuovo...@gmail.com 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]]

Re: [R] Approximation to next integer

2009-06-13 Thread Carletto Rossi
thanks. Solved! 2009/6/14 Henrique Dallazuanna www...@gmail.com See ?round On Sat, Jun 13, 2009 at 8:53 PM, Carletto Rossi nuovo...@gmail.comwrote: 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??

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

2009-06-13 Thread Glenn Woodart
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

[R] color in pairs

2009-06-13 Thread Naoki Irye
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

Re: [R] color in pairs

2009-06-13 Thread Henrique Dallazuanna
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.cor) On

[R] time function behavior for ts class objects

2009-06-13 Thread Ron Burns
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

Re: [R] Gtk objects disappears

2009-06-13 Thread Michael Lawrence
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 onu...@iberstat.es wrote: Dear John, I have a question. When I run a RGtk code in my terminal (without

Re: [R] triangular plot

2009-06-13 Thread Greg Snow
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, June

Re: [R] How to write loop

2009-06-13 Thread Alex Roy
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

[R] How to capture output from a subprocess?

2009-06-13 Thread Ronggui Huang
## 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'. I try

Re: [R] How to write loop

2009-06-13 Thread Johannes Huesing
Alex Roy alexroy2...@gmail.com [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