[R] about the 95%CI around the median...

2008-08-03 Thread Fernando Marmolejo Ramos
Dear people I've learnt that by using the "boxplot.stats" command in the "grDevices" library I can get the 5-number summaries of a boxplot, plus other important information, like the confidence interval around the median. I'm interested in knowing the actual formula to used in that package to cal

[R] Decomposing tests of interaction terms in mixed-effects models

2008-08-03 Thread Andrew Robinson
Dear R colleagues, a friend and I are trying to develop a modest workflow for the problem of decomposing tests of higher-order terms into interpretable sets of tests of lower order terms with conditioning. For example, if the interaction between A (3 levels) and C (2 levels) is significant, it ma

[R] Sweave and ggplot2

2008-08-03 Thread Sorn . Norng
Hi all, I've been trying to run Sweave with R code embedded - using the ggplot2 package and in particular the qplot command. There appears to be a problem in Sweave not picking up that qplot is a function. Has anybody else tried to use qplot in Sweave and have you been successful? Any help woul

Re: [R] Changing values

2008-08-03 Thread J Dougherty
On Sunday 03 August 2008 01:54:43 pm Andrew Ramsey wrote: > Hello-- > > I am a relatively new user to R and I cannot find the information I > need. Please help. > > I have a very large data set with values including letters, numbers, > and symbols (sometimes within the same vector value [ie X9-]. >

Re: [R] Determining model parameters

2008-08-03 Thread Ben Bolker
~ The beta distribution only applies to data that are bounded between 0 and 1 (in some cases strictly bounded, i.e. 0 wrote: |> charter.net> writes: |> |>> If I have some data that based on the historgram and other plots it "looks" |> like a beta distribution. Is there |>> a function or functi

Re: [R] Bubble plots

2008-08-03 Thread Michael Bibo
Cody Hamilton Edwards.com> writes: > > Is there a way to create a 'bubble plot' in R? > > For example, if we define the following data frame containing the level of y observed for 5 patients at three > time points: > > time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5)) > y<-c('a','b','c

Re: [R] Changing values

2008-08-03 Thread jim holtman
Exactly how is the translation of "X9-" -> 00911 done? Are there unique mappings of character sequences to numbers? How many different ones might there be? Why do you have leading zeros on the result? If they are changed to numeric, then the default printing results in '911'. Can you provide a

[R] Changing values

2008-08-03 Thread Andrew Ramsey
Hello-- I am a relatively new user to R and I cannot find the information I need. Please help. I have a very large data set with values including letters, numbers, and symbols (sometimes within the same vector value [ie X9-]. I've imported the data using read.fwp and it arrives in list fo

Re: [R] Determining model parameters

2008-08-03 Thread Ben Bolker
charter.net> writes: > If I have some data that based on the historgram and other plots it "looks" like a beta distribution. Is there > a function or functions within R to help me determine the model parameters for such a distirbution? library(MASS) ?fitdistr _

[R] Determining model parameters

2008-08-03 Thread rkevinburton
This may be a begining question. If so, please bear with me. If I have some data that based on the historgram and other plots it "looks" like a beta distribution. Is there a function or functions within R to help me determine the model parameters for such a distirbution? Similarily for other "

Re: [R] Convert date to decimal days

2008-08-03 Thread cls59
You're right, the problem was with the print function not displaying enough precision. Thanks! -- View this message in context: http://www.nabble.com/Convert-date-to-decimal-days-tp18800462p18802117.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Convert date to decimal days

2008-08-03 Thread Prof Brian Ripley
I think you will find that the 'catastrophic loss of significance' is in the printing at the default level. Try print(x, digits=15). On Sun, 3 Aug 2008, cls59 wrote: Hello all, I have a quick question about formatting date strings.. Datetime strings, methinks. If you do just want days, u

Re: [R] convert for loop into apply()

2008-08-03 Thread Charles C. Berry
If length( levels( a1$cat ) ) is not very large, and the structure suggested in the toy example holds in the actual case (probes in a2 do not 'overlap' and are in order according to a2$st and a2$en, then this might do: unsplit(lapply(levels(a1$cat), function(x) { + tmp1 <- subset( a1, ca

[R] Convert date to decimal days

2008-08-03 Thread cls59
Hello all, I have a quick question about formatting date strings.. I am currently debugging some Matlab code someone else wrote and since it is so bad that I have to go through it line by line I figured that I would just rewrite the thing in R. The code produces plots of wave spectra with decim

Re: [R] Will This Computer Run 64 bit Ubuntu/R?

2008-08-03 Thread Prof Brian Ripley
On Sun, 3 Aug 2008, Dirk Eddelbuettel wrote: On Sun, Aug 03, 2008 at 05:38:01PM +0100, Prof Brian Ripley wrote: The Pemtium D 940 supports EMT64 according to http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it). That is the key: EMT64 is Intel's name for x86_64 and so your box

[R] missing F statistic in anova.gam

2008-08-03 Thread Michael A. Milligan
Hello, I have encountered results which I am not sure how to interpret when using anova.gam to compare 2 different models. For certain tests the results do not include an F- or associated p-statistic. This happens when comparing certain models and not others, and I do not discern a patten exp

Re: [R] access object inside a function

2008-08-03 Thread Duncan Murdoch
On 03/08/2008 12:13 PM, Jonas wrote: Duncan, Patrick, Sorry for not being very clear (and having non-working code in my example). What I'm actually trying to achieve is creating a wrapper around the NelsonSiegel function in the fBonds package in order to (among other things) create my own plot f

Re: [R] Will This Computer Run 64 bit Ubuntu/R?

2008-08-03 Thread Dirk Eddelbuettel
On Sun, Aug 03, 2008 at 05:38:01PM +0100, Prof Brian Ripley wrote: > The Pemtium D 940 supports EMT64 according to > http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it). > > That is the key: EMT64 is Intel's name for x86_64 and so your box should > run x86_64 Linux which I belie

Re: [R] Will This Computer Run 64 bit Ubuntu/R?

2008-08-03 Thread Prof Brian Ripley
The Pemtium D 940 supports EMT64 according to http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it). That is the key: EMT64 is Intel's name for x86_64 and so your box should run x86_64 Linux which I believe Debian/Ubuntu refer to as 'amd64' (which understandably Intel does not f

Re: [R] Bubble plots

2008-08-03 Thread Franz Mueter
I believe ?symbols is what you are looking for. Of course, you need to convert categories to numeric values for plotting: For example, adding some actual data to your data frame: > D <- cbind(D, P = runif(15)) > symbols(as.numeric(D$time), as.numeric(D$y), circles=D$P, inches=0.2, ann=F, xaxt="n

Re: [R] (sans objet)

2008-08-03 Thread Dirk Eddelbuettel
On Sun, Aug 03, 2008 at 06:19:54PM +0200, pir2.jv wrote: > > - Peut-on poser des questions en français? Pas vraiment. R-help est une liste organisee en anglais. > - Très simple. Je fais un programme "truc.r" (sous emacs) > Ce programme réalise des graphiques. > Je voudrais les enregistrer au fur

Re: [R] drop1() seems to give unexpected results compare to anova()

2008-08-03 Thread Prof Brian Ripley
On Sun, 3 Aug 2008, Thomas P C Chu wrote: Thanks to Mr Dalgaard for his advice and everyone else who has contributed. Inclusion of an error term at the end of sim.set$y = ... line did cure my problems with drop1() and step(). I suppose it is my own inexperience in carrying out simulations cau

[R] (sans objet)

2008-08-03 Thread pir2.jv
- Peut-on poser des questions en français? - Très simple. Je fais un programme "truc.r" (sous emacs) Ce programme réalise des graphiques. Je voudrais les enregistrer au fur et à mesure -- sinon, quand la session est finie, ces graphiques sont perdus-- J'ai essayé: pdf("truc.pdf") Mais le fichi

Re: [R] access object inside a function

2008-08-03 Thread Jonas
Duncan, Patrick, Sorry for not being very clear (and having non-working code in my example). What I'm actually trying to achieve is creating a wrapper around the NelsonSiegel function in the fBonds package in order to (among other things) create my own plot function. However, all of the objects i

Re: [R] access object inside a function

2008-08-03 Thread Patrick Burns
Duncan Murdoch wrote: On 03/08/2008 11:29 AM, Jonas wrote: hi, my apologies if this has been covered numerous times before and it's only my lack of search and/or R skills that is stopping me from finding the solution. i'm trying to access an object defined and created inside a function (the ob

[R] gmaps, projection, inset

2008-08-03 Thread John P. Burkett
Running R version 2.6.1 under Gentoo Linux, I'm trying to produce a thematic map of the USA using the gmaps package. The result thus far has two problems from my point of view. First, the projection (Miller cylindrical?) elongates southern states and flattens northern ones unattractively. I'd

Re: [R] running strucchange?

2008-08-03 Thread Burgs
Thanks to all for your replies. It is appreciated. Prof Brian Ripley wrote: > > On Sat, 2 Aug 2008, Burgs wrote: > >> >> Yes, I'm using Windows. I'm getting an error when I use the source() >> command, it seems not to like the path to the .R file I'm putting in, >> although it is the correc

Re: [R] access object inside a function

2008-08-03 Thread Duncan Murdoch
On 03/08/2008 11:29 AM, Jonas wrote: hi, my apologies if this has been covered numerous times before and it's only my lack of search and/or R skills that is stopping me from finding the solution. i'm trying to access an object defined and created inside a function (the object is not returned by

[R] access object inside a function

2008-08-03 Thread Jonas
hi, my apologies if this has been covered numerous times before and it's only my lack of search and/or R skills that is stopping me from finding the solution. i'm trying to access an object defined and created inside a function (the object is not returned by the function), but i can't seem to get

[R] Installing R on Apache Web Server

2008-08-03 Thread Ajay ohri
Hi List, I have run R from Windows Terminal server and it works fine. I now need to install it on apache Web server / Amazon EC 2 server to test R (with a GUI) for computing on a cloud (using a framework I wrote about).Intent is to benchmark R 's superior algorithms on very big files > 1-2 gb (aga

[R] Will This Computer Run 64 bit Ubuntu/R?

2008-08-03 Thread Tom La Bone
After doing some reading about 64-bit systems and software I am still somewhat uncertain about some things. I have a Dell Dimension XPS 400 with a dual core Intel Pentium D 940 (3.2 GHz) and 4 Gb of memory. I currently dual boot the system with XP Professional and Ubuntu 8.04 (32 bit). If I simply

[R] neighbors of an element of a matrix

2008-08-03 Thread rostam shahname
Hi R users, I wonder if there is any function which would render the value of the neighbors of the given element [i,j] of a matrix. Thanks, Rostam [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] losing row.names in matrix operations

2008-08-03 Thread Gabor Grothendieck
Constructing a zoo object is well covered in ?zoo and in the three vignettes: vignette("zoo") vignette("zoo-quickref") vignette("zoo-faq") On Sat, Aug 2, 2008 at 10:27 AM, rcoder <[EMAIL PROTECTED]> wrote: > > Thanks for your reply Gabor. > > As I already have a dates column in my dataframe (in

Re: [R] How do I get an inset graph (i.e. graph within a graph)?

2008-08-03 Thread Gabor Grothendieck
On Sun, Aug 3, 2008 at 3:51 AM, Arthur Roberts <[EMAIL PROTECTED]> wrote: > How do I get an inset graph (i.e. graph within a graph)? Your input is > greatly appreciated. See ?subplot in the TeachingDemos package. __ R-help@r-project.org mailing list ht

Re: [R] How do I get an inset graph (i.e. graph within a graph)?

2008-08-03 Thread Mark Difford
Hi Arthur, This can be done quite easily using the appropriate arguments listed under ?par; and there are other approaches. Ready-made functions exist in several packages. I tend to use ?add.scatter from package ade4. It's a short function, so it's easy to customize it, but it works well straight

Re: [R] drop1() seems to give unexpected results compare to anova()

2008-08-03 Thread Thomas P C Chu
Thanks to Mr Dalgaard for his advice and everyone else who has contributed. Inclusion of an error term at the end of sim.set$y = ... line did cure my problems with drop1() and step(). I suppose it is my own inexperience in carrying out simulations caused such gaffe. Thomas __

Re: [R] Gaps in time series.

2008-08-03 Thread Gabor Grothendieck
If you are using zoo for this you can do: library(zoo) set.seed(1) x <- zoo(rnorm(7), as.numeric(1:7)) y <- zoo(rnorm(4), c(1, 2, 6, 7)) xy <- merge(x,y) xy$y - xy$x # or if you just want the difference at times existing in both y - x On Sat, Aug 2, 2008 at 5:16 PM, <[EMAIL PROTECTED]> wrote:

Re: [R] output components of GAM

2008-08-03 Thread Gavin Simpson
is this gam::gam or mgcv::gam? If the latter, and if I understand what you want, if gam.out contains your GAM model fitted using gam() from mgcv, then: out.term <- predict(gam.out, type = "terms") gives you the contributions of each covariate on the response for the fitted model. read :?predict.g

Re: [R] R command history -- can it be like Matlab's?

2008-08-03 Thread Peter Dalgaard
Not quite true that you can't type anything. What happens (for me) is that you are still in reverse-i-search, so you can get this effect from "^R l ^C d". > (reverse-i-search)`l': ls() (reverse-i-search)`l': ls() > (reverse-i-search)`l': ls() (reverse-i-search)`ld': levels(ftpain3) <- list

Re: [R] R command history -- can it be like Matlab's?

2008-08-03 Thread Gad Abraham
Peter Dalgaard wrote: Gad Abraham wrote: Prof Brian Ripley wrote: A patch to do this was posted on 2007-09-29 by Glenn Davis. Some people not addicted to Matlab find the behaviour very inconvenient and prefer the getline/readline behaviour (triggered by ^R/^S) of Rterm and R on Unixen. Now

Re: [R] boxplot help

2008-08-03 Thread Bernardo Rangel Tura
Em Sex, 2008-08-01 às 14:16 -0700, Rajasekaramya escreveu: > hi > > I have list of matrix of lenggth 61 containg the mean values..I want to make > a boxplot for each of the matrix. > I used a for loop but i cant figure out the way to save in the boxplots > > > all.the.mean > [[1]] > mean > 0.5 >

Re: [R] convert for loop into apply()

2008-08-03 Thread Hans W. Borchers
Also, your request can easily be formulated as an SQL statement, for example utilizing the 'sqldf' package: library(sqldf) a1 <- data.frame(id = 1:6, cat = paste('cat', rep(1:3, c(2,3,1))), st = c(1, 7, 30, 40, 59, 91), en = c(5, 25, 39, 55

Re: [R] Smartest way to evaluate question forms

2008-08-03 Thread vlasto
First of all, thank you for your reply and for the links. I see that my problem might be overly detailed in its description and not very clear in its outlines. In fact the question was: If there is any existing module for R to evaluate data, with parameters being only boolean values. (Which u ge

[R] output components of GAM

2008-08-03 Thread Andreas Winter
I would like to request help with the following: I am trying to use a Generalized Additive Model (gam) to examine the density distribution of fish as a function of latitude and longitude as continuous variables, and year as a categorical variable. The model is written as:   gam.out   <-  gam(Dens

Re: [R] losing row.names in matrix operations

2008-08-03 Thread rcoder
Thanks for your reply Gabor. As I already have a dates column in my dataframe (in column row.names), is it possible to preserve this whilst still making a data set suitable for rollapply()? Thanks, rcoder Gabor Grothendieck wrote: > > If its regular you can convert it to ts or zoo. > If its

[R] problem with nested loop for regression

2008-08-03 Thread rcoder
Hi everyone, I'm experiencing difficulty getting the results I want when I use a nested for loop. I have a data set to which I perform some calculations, and then try to apply a regression over a rolling window. The code runs, but the regression results I am getting (intercept and slope) are simp

[R] [R-pkgs] sensitivity

2008-08-03 Thread IOOSS Bertrand 165542
Hi All, the new release v1.4-0 of the "sensitivity" package is now available on CRAN The "sensitivity" package is devoted to factor screening and global sensitivity analysis of numerical model output. Here are the new features list: New functionalities: * sequential bifurcation method

[R] [R-pkgs] New package: noia

2008-08-03 Thread Arnaud Le Rouzic
Hi the list, A new version (0.92) of my package 'noia' will be available soon on CRAN mirrors, and I think it might be a good opportunity to introduce it shortly to the R community. In summary: 'noia' will be of absolutely no interest for 99.99% of you. The 0.01% remaining are quantitative g

[R] How do I get an inset graph (i.e. graph within a graph)?

2008-08-03 Thread Arthur Roberts
Hi, all, How do I get an inset graph (i.e. graph within a graph)? Your input is greatly appreciated. Best wishes, Art University of Washington Department of Medicinal Chemistry __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list