Re: [R] regular expression

2007-04-08 Thread Laurent Rhelp
Gabor Grothendieck a écrit : On 4/7/07, Laurent Rhelp [EMAIL PROTECTED] wrote: Gabor Grothendieck a écrit : I assume is to replace single quotes in double quoted strings with single quoted strings containing double quotes in place of the single quotes that's rignt

Re: [R] regular expression

2007-04-08 Thread Gabor Grothendieck
I can't reproduce that error. For me it works fine. What version of R and gsubfn are you using? I switched to a French locale just to be sure. I tried pasting it into an R session and also sourcing it but both worked fine for me. Sys.getlocale() [1] LC_COLLATE=English_United

[R] Simulation of the Frailty of the Cox PH model

2007-04-08 Thread Mohammad Ehsanul Karim
Dear R-list users, I am trying to do simulation of survival data to enable it to run under frailty option. Below is the function a that I am using. My questions are: 1. How do I modify it to get bigger (hopefully significant) value of Variance of random effect? 2. What changes do I have to make

Re: [R] regular expression

2007-04-08 Thread Laurent Rhelp
Gabor Grothendieck a écrit : I can't reproduce that error. For me it works fine. What version of R and gsubfn are you using? I switched to a French locale just to be sure. I tried pasting it into an R session and also sourcing it but both worked fine for me. Sys.getlocale() [1]

Re: [R] regular expression

2007-04-08 Thread Laurent Rhelp
Gabor Grothendieck a écrit : On 4/7/07, Laurent Rhelp [EMAIL PROTECTED] wrote: Gabor Grothendieck a écrit : I assume is to replace single quotes in double quoted strings with single quoted strings containing double quotes in place of the single quotes that's rignt

[R] Labelling boxplot with fivenumber summary

2007-04-08 Thread Daniel Siddle
Just wanted to say thanks very much. I used Chuck's 2nd idea as I found it the easiest to understand as I'm still finding me feet with R. Just for reference for anyone else, fn[1] and fn[5] actually pasted and gave the values at the whiskers (~1.5 IQR) so replaced them with a max and min

Re: [R] regular expression

2007-04-08 Thread Gabor Grothendieck
Your example is a vector of strings not a string with mulitple lines. To turn it into a string with mulitple lines use paste: gsubfn('[][^]*[]', f, paste(Lines, collapse = \n)) If you want to break up the output to create a vector of strings rather than a single multi-line string see ?strsplit

[R] data encapsulation with classes

2007-04-08 Thread Peter Williams
Hi All, I'm new to R from a C and Octave/Matlab background. I am trying to construct some classes in R to which I want to attach pieces of data. First, is attr(obj, 'member name') - data the accepted way of doing this? Second, having attached member data with attr is there anything wrong

Re: [R] data encapsulation with classes

2007-04-08 Thread Roger Bivand
On Sun, 8 Apr 2007, Peter Williams wrote: Hi All, I'm new to R from a C and Octave/Matlab background. I am trying to construct some classes in R to which I want to attach pieces of data. First, is attr(obj, 'member name') - data the accepted way of doing this? No, it isn't. You seem

Re: [R] data encapsulation with classes

2007-04-08 Thread hadley wickham
On 4/8/07, Roger Bivand [EMAIL PROTECTED] wrote: On Sun, 8 Apr 2007, Peter Williams wrote: Hi All, I'm new to R from a C and Octave/Matlab background. I am trying to construct some classes in R to which I want to attach pieces of data. First, is attr(obj, 'member name') - data the

Re: [R] Labelling boxplot with fivenumber summary

2007-04-08 Thread Ranjan Maitra
You may wish to provide brief commented source code: that would be most helpful for people referring to it in the archives. Many thanks and best wishes, Ranjan On Sun, 8 Apr 2007 11:43:47 + Daniel Siddle [EMAIL PROTECTED] wrote: Just wanted to say thanks very much. I used Chuck's 2nd

Re: [R] compile errors with rgl-0.70.570 on FreeBSD

2007-04-08 Thread Bernardo Rangel Tura
On Wed, 2007-04-04 at 21:12 +0200, Rainer Hurling wrote: Thank you for your work on rgl. Reading in another thread about your new (inofficial) version of rgl (see below) I tried it with R version 2.5.0 alpha (2007-03-31 r40986) under FreeBSD 7.0-CURRENT. Unfortunately I got the following

Re: [R] compile errors with rgl-0.70.570 on FreeBSD

2007-04-08 Thread Rainer Hurling
Bernardo, thank you for your tips. But before version 0.70.570 rgl was running fine on my system. All necessary GL files are installed. Duncan Murdoch is working on a configure script to get rgl also work on Mac OSX. Some of his changes breaks things for FreeBSD. We tested it without posting

[R] Confidence-Intervals.... help...

2007-04-08 Thread Jochen.F
Hi... I have to use R to find out the 90% confidence-interval for the sensitivity and specificity of the following diagnostic test: A particular diagnostic test for multiple sclerosis was conducted on 20 MS patients and 20 healthy subjects, 6 MS patients were classified as healthy and 8 healthy

Re: [R] Confidence-Intervals.... help...

2007-04-08 Thread Sarah Goslee
Hm... sounds like a homework problem to me... Maybe start by figuring out how to do it without R - what's the approach, and how would you calculate it? Then search R help for the possible key words you came up with. Sarah On 4/8/07, Jochen.F [EMAIL PROTECTED] wrote: Hi... I have to use R to

[R] Failure of mcsamp() but not mcmcsamp()

2007-04-08 Thread Michael Kubovy
Daer r-helpers, Can anyone help with the following: I get: mcmcsamp(txt.lmer3f) Markov Chain Monte Carlo (MCMC) output: Start = 1 End = 1 Thinning interval = 1 But: mcsamp(txt.lmer3f) Error in as.bugs.array(sims, program = lmer, n.iter = n.iter, n.burnin = n.burnin, :

Re: [R] Reasons to Use R

2007-04-08 Thread Johann Hibschman
On 4/6/07, Wilfred Zegwaard [EMAIL PROTECTED] wrote: I'm not a programmer, but I have the experience that R is good for processing large datasets, especially in combination with specialised statistics. This I find a little surprising, but maybe it's just a sign that I'm not experienced enough

Re: [R] Reasons to Use R

2007-04-08 Thread Gabor Grothendieck
On 4/8/07, Johann Hibschman [EMAIL PROTECTED] wrote: R's pass-by-value semantics also make it harder than it should be to deal with where it's crucial that you not make a copy of the data frame, for fear of running out of memory. Pass-by-reference would make implementing data transformations

[R] [Fwd: Re: regular expression]

2007-04-08 Thread Laurent Rhelp
---BeginMessage--- Uwe Ligges a écrit : I guess your problem has been solved by last night's discussion with Gabor G.? Uwe Ligges Laurent Rhelp wrote: Uwe Ligges a écrit : Laurent Rhelp wrote: Uwe Ligges a écrit : Laurent Rhelp wrote: Dear R-List, I have a great

[R] How do I back transforme ordinary log-krigged prdiction values?

2007-04-08 Thread Zia Uddin Ahmed
I have a question to everybody. After log10 transfprmation, I have done ordinary kriging in gstat in R? I need to back trnasform the prediction values to orgiginal scale. How do I do this in gstat in R? Thanks Zia -- Zia Uddin Ahmed 915 Brad Field Hall Department of Crop and Soil Cornell

[R] Plot symbols dimensions

2007-04-08 Thread Cressoni, Massimo \(NIH/NHLBI\) [F]
I am writing some code to obtain publication-like plots (like the ones can be obtained with SigmaPlot). I am not able to find the dimensions of a point. Parameters like cex make you able to make it bigger or smaller but I need to know the exact dimension (in relative coordinates : if my plot

Re: [R] Reasons to Use R

2007-04-08 Thread Wilfred Zegwaard
Dear Johann and Gabor, It's what amounts to large datasets. There are hundreds of datasets R can't handle, probably thousands or more. I noticed on my computer (which is nothing more that an average PC) that R breaks down after 250 MB of memory. I also note that SPSS breaks down, Matlab, etc.

[R] Could not fit correct values in discriminant analysis by bruto.

2007-04-08 Thread 川口 修治
Dear R-users, I would like to use bruto function in mda package for flexible discriminant analysis. Then, I tried, for example, following approach. x band1 band2band3 1 -1.206780 -1.448007 -1.084431 2 -0.294938 -0.113222 -0.95

[R] Repeated Measures design using lme

2007-04-08 Thread Scott Norton
Hi, I have what I believe is a repeated-measures dataset that I'm trying to analyze using lme(). This is *not* homework, but an exercise in my trying to self-teach myself repeated-measure ANOVA for other *real* datasets that I have and that are extremely similar to the following design. I'm

[R] R:Maximum likelihood estimation using BHHH and BFGS

2007-04-08 Thread joey repice
Dear R users, I am new to R. I would like to find *maximum likelihood estimators for psi and alpha* based on the following *log likelihood function*, c is consumption data comprising 148 entries: fn-function(c,psi,alpha) { s1-sum(for(i in 1:n){(c[i]-(psi^(-1/alpha)*(lag(c[i],-1^2*

[R] Dealing with large nominal predictor in sem package

2007-04-08 Thread adschai
Hi, I am using tsls function from sem package to estimate a model which includes large number of data. Among its predictors, it has a nominal data which has about 10 possible values. So I expand this parameter into 9-binary-value predictors with the coefficient of base value equals 0. I also