Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread Tal Galili
It would best be there - I agree. But since I had no experience with the site - I didn't know the "rules of use" for it. Thanks for clarifying. Tal -- My contact information: Tal Galili E-mail: tal.gal...@gmail.com Phone number: 972-52-7275845 FaceB

Re: [R] Minimum and blanks

2009-09-30 Thread Peter Ehlers
Chris, What is your definition of the minimum of a *character* vector? -Peter Ehlers Chris Li wrote: Hi all, I want to calculate the minimum of a column which contains blanks. R returns "" as the minimum, which is not I want. Is there a way to overcome it? Thanks in advance. Chris _

[R] aproximate a titration kurve to the measure data.

2009-09-30 Thread awayguy
Halo i'm studying chemistry, today we made an experiment and i have to draw a titration kurve for my mess data. we can do it on a mm paper, or we can also use a programe. people from chemistry recomend "R" last year i studied civil eng. and we used Matlab, as I see, R ist very similar to it, but

Re: [R] re ading and analyzing a word document

2009-09-30 Thread PDXRugger
Considering your instructions: #Define words to find to.find <- c( 'the', 'is', 'are' ,'dr') #Read in the file... file.text <- readLines( 'data/letter.txt' ) #Count number of occurnces of deined word in text line.matches <- unlist( lapply( to.find, grep, x = unlist(file.text[2]) ) ) Result:

[R] Input format for kcca in kernlab

2009-09-30 Thread Guru Ananda
Hi, I'm trying to use kcca function from the kernlab package, and it requires 2 input datasets, both of which need to be matrices containing data index by row. I have two data frame objects which I converted to matrices using data.frame(x). Their dimensions are 2264*10 and 2264*4. I then tried ru

[R] re ading and analyzing a word document

2009-09-30 Thread PDXRugger
file=("LUSDR/letter.doc") Howdy Y'all, So i am looking to read a word document in the following formats(.doc) or any type of accessible word processor software (e.g. text .txt, notepad, etc). Had the ability to search certain words, for instance "banana", "peacock","Weapons" "Mass" "Destructio

Re: [R] Calling R functions into C# or C++

2009-09-30 Thread fayssal
Hi, It is possible to call R functions from .NET by using a .NET wrapper over the R (D)Com server. http://www.nabble.com/file/p25682136/RFromCsharp.zip Attached is a sample C# project that shows how to call R from C#. Best Regards, Fayssal El Moufatich http://www.nabble.com/file/p25682136/RF

[R] Confidence intervals PLS prediction

2009-09-30 Thread lievelaurens
I have switched from The Unscrambler to R for pls regression analysis and have been able to calculate scores, coefficients, RMSEP from a large number of PLS1 and PLS2 models. The ultimate goal is to use these models for predicting unknown samples, which again is straight-forward with the built-in

[R] Minimum and blanks

2009-09-30 Thread Chris Li
Hi all, I want to calculate the minimum of a column which contains blanks. R returns "" as the minimum, which is not I want. Is there a way to overcome it? Thanks in advance. Chris -- View this message in context: http://www.nabble.com/Minimum-and-blanks-tp25692189p25692189.html Sent from th

[R] xyplot key

2009-09-30 Thread Andrewjohnclose
Hi, I'm having trouble matching the symbols/color of the key to match those specifiec in the plot. Here is the code I used: xyplot(GCR+GCT ~ FRAC, data=RWF, type=c("g","p"),cex=1.2,pch=c(22,21), xlab=expression(italic("f")),ylab=expression(italic("S")), key=list(points=list(pch=c(22,21),cex=1.

[R] Missing functions

2009-09-30 Thread Kenny Shen
Hi all, I'm new to R and have been working hard to get familiarized with it. A problem I'm facing now is that having installed some packages (psych, doBy), I can't seem to access the functions even through there was no error messages when I load them using library(). I get an error telling me the

[R] Off-topic: looking for a categorical, NMAR dataset

2009-09-30 Thread Andrew Robinson
Dear Colleagues, apologies for this off-topic posting. A Ph.D. student here at U of Melb. is trying to find a dataset to use to demonstrate a technique that he is developing. He needs a binary response and ideally a categorical predictor, although the latter can of course be induced from a conti

Re: [R] Re gression for levels of a factor/xyplot type="r"

2009-09-30 Thread Jorge Ivan Velez
Hi James, Here is a suggestion: # Data set.seed(123) Factor <- rep(1:4, each = 20) X <- rnorm(80, 25, 4) Y <- 5 + 1.5*Factor*X + rnorm(80) df <- data.frame(X, Y, Factor) # Splitting df sdf <- with(df, split(df, Factor)) # betas and R2 do.call(rbind, lapply(sdf, function(miniDF){

[R] mvmeta function in GAP package

2009-09-30 Thread Yan Li
Hi, I am using mvmeta function in GAP package. Could you please tell me the algorithm applied in this function and the relevant references? Thanks. Lee [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.eth

Re: [R] Re gression for levels of a factor/xyplot type="r"

2009-09-30 Thread Gabor Grothendieck
See lmList in the lme4 package. Each component of the result will be one lm and you can take a summary of each. On Wed, Sep 30, 2009 at 11:10 PM, jimdare wrote: > > Hi, > > I'm sure these are basic problems so I apologise in advance for my > ignorance.  I have a dataset with X, Y, and a Factor w

Re: [R] re ading and analyzing a word document

2009-09-30 Thread cls59
PDXRugger wrote: > > Considering your instructions: > > #Define words to find > to.find <- c( 'the', 'is', 'are' ,'dr') > #Read in the file... > file.text <- readLines( 'data/letter.txt' ) > #Count number of occurnces of deined word in text > line.matches <- unlist( lapply( to.find, grep, x

Re: [R] R package for visualizing/analyzing accelerometry data?

2009-09-30 Thread spencerg
library(sos) a <- ???accelerometry # no matches d3 <- ???'3D' # 411 matches; retrieved 20 pages 400 matches d3 <- ???'3D'(99) # retrieved all 411 matches summary(d3) # first 13 of 144 packages findFn2xls(d3) # write d3.xls containing sheets PackageSum2, findFn, call installPackages(d3) # install

[R] Re gression for levels of a factor/xyplot type="r"

2009-09-30 Thread jimdare
Hi, I'm sure these are basic problems so I apologise in advance for my ignorance. I have a dataset with X, Y, and a Factor with 4 levels. I am trying to figure out how to use lm() to find the r2 value, slope and intercept of an X~Y regression, for each level of the Factor. Is there a way to

Re: [R] re ading and analyzing a word document

2009-09-30 Thread cls59
PDXRugger wrote: > > Howdy Y'all, > > So i am looking to read a word document in the following formats(.doc) or > any type of accessible word processor software (e.g. text .txt, notepad, > etc). Had the ability to search certain words, for instance "banana", > "peacock","Weapons" "Mass" "Dest

Re: [R] matrix exponential

2009-09-30 Thread David Winsemius
On Sep 30, 2009, at 10:14 PM, spencerg wrote: RSiteSearch('exponential of a matrix') produced 982 matches. RSiteSearch('{exponential of a matrix}') produced 13. Yes, I know. David Winsemius wrote: On Sep 30, 2009, at 9:18 PM, Kon Knafelman wrote: Hi Guys, Im trying to fin

Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread Ben Bolker
Tal Galili gmail.com> writes: > > Hi Joris, > > Good luck with your work. > > I know how to set up a wiki. But I wouldn't do it for just 1 page. And also, > if you have little experience with it - I am not sure you would find it easy > to "jump" into it. > I'll have a look around to see what o

Re: [R] matrix exponential

2009-09-30 Thread spencerg
RSiteSearch('exponential of a matrix') produced 982 matches. RSiteSearch('{exponential of a matrix}') produced 13. David Winsemius wrote: On Sep 30, 2009, at 9:18 PM, Kon Knafelman wrote: Hi Guys, Im trying to find the exponential of a matrix. Can someone please help me

Re: [R] matrix exponential

2009-09-30 Thread cls59
Kon Knafelman wrote: > > > Hi Guys, > > Im trying to find the exponential of a matrix. > > Can someone please help me do this? > > Thanks a lot > > There was actually a very recent discussion of this, and related operations. See: http://www.nabble.com/imple

Re: [R] matrix exponential

2009-09-30 Thread David Winsemius
On Sep 30, 2009, at 9:18 PM, Kon Knafelman wrote: Hi Guys, Im trying to find the exponential of a matrix. Can someone please help me do this? Didn't we just do this one? -- David Winsemius, MD Heritage Laboratories West Hartford, CT __ R-help

[R] matrix exponential

2009-09-30 Thread Kon Knafelman
Hi Guys, Im trying to find the exponential of a matrix. Can someone please help me do this? Thanks a lot _ View photos of singles in your area Click Here =1&lage=18&uage=55&cl=14

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: Peter Dalgaard [mailto:p.dalga...@biostat.ku.dk] > Sent: Wednesday, September 30, 2009 3:58 PM > To: Nordlund, Dan (DSHS/RDA) > Cc: Douglas Bates; r help > Subject: Re: [R] Rounding error in seq(...) > <<>> > > > > The first page of this document references: >

Re: [R] (windows xp) start script on startup / double clicking on the script

2009-09-30 Thread Gabor Grothendieck
Try creating a Windows batch file along these lines: setlocal set R_PROFILE_USER=C:\tmp\myscript.R "C:\Program Files\R\R-2.9.x\bin\Rgui.exe" endlocal and double click it. The set line sets it up to run your script and the next line runs R. On Wed, Sep 30, 2009 at 7:12 PM, Martin Batholdy wrot

[R] (windows xp) start script on startup / double clicking on the script

2009-09-30 Thread Martin Batholdy
hi, I am still trying to figure out how it is possible to start an R- script via the R-GUI directly by double-clicking on the script file in windows xp. On Mac OS I have the option "start script in the editor" in the preferences. But on windows XP I always get the "ARGUMENT 'c:\...\...'

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Peter Dalgaard
Nordlund, Dan (DSHS/RDA) wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Douglas Bates Sent: Wednesday, September 30, 2009 3:06 PM To: Peter Dalgaard Cc: r help; Duncan Murdoch Subject: Re: [R] Rounding error in seq(...) On

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Erik Iverson
> In an earlier thread on this theme I believe that someone quoted Brian > Kernighan as saying "10 times 0.1 is hardly ever 1" but I haven't been > able to track down the quote. Can anyone point us to such a quote? > It summarizes the situation succinctly, See FAQ 7.31 :) ___

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Douglas Bates > Sent: Wednesday, September 30, 2009 3:06 PM > To: Peter Dalgaard > Cc: r help; Duncan Murdoch > Subject: Re: [R] Rounding error in seq(...) > > On Wed, Sep 30, 200

Re: [R] Rounding error in seq(...)

2009-09-30 Thread David Winsemius
On Sep 30, 2009, at 6:05 PM, Douglas Bates wrote: On Wed, Sep 30, 2009 at 2:32 PM, Peter Dalgaard> wrote: Martin Batholdy wrote: hum, can you explain that a little more detailed? Perhaps I miss the background knowledge - but it seems just absurd to me. 0.1+0.1+0.1 is 0.3 - there is no ro

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Douglas Bates
On Wed, Sep 30, 2009 at 2:32 PM, Peter Dalgaard wrote: > Martin Batholdy wrote: >> >> hum, >> >> can you explain that a little more detailed? >> Perhaps I miss the background knowledge - but it seems just absurd to me. >> >> 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? >> >> why is

Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread Tal Galili
Hi Joris, Good luck with your work. I know how to set up a wiki. But I wouldn't do it for just 1 page. And also, if you have little experience with it - I am not sure you would find it easy to "jump" into it. I'll have a look around to see what other collaborative tools there are out there. Plea

Re: [R] preformatted and '#' in manual pages

2009-09-30 Thread Gábor Csárdi
On Wed, Sep 30, 2009 at 10:51 PM, Duncan Murdoch wrote: > On 9/29/2009 7:31 AM, Gábor Csárdi wrote: >> >> Uwe, thanks, but this does not help, I still get: >> >> LaTeX errors when creating PDF version. >> This typically indicates Rd problems. >> LaTeX errors found: >> ! You can't use `macro parame

Re: [R] preformatted and '#' in manual pages

2009-09-30 Thread Duncan Murdoch
On 9/29/2009 7:31 AM, Gábor Csárdi wrote: Uwe, thanks, but this does not help, I still get: LaTeX errors when creating PDF version. This typically indicates Rd problems. LaTeX errors found: ! You can't use `macro parameter character #' in vertical mode. ...ike the following: \begin {alltt} ##

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Ted Harding
On 30-Sep-09 19:32:46, Peter Dalgaard wrote: > Martin Batholdy wrote: >> hum, >> >> can you explain that a little more detailed? >> Perhaps I miss the background knowledge - but it seems just absurd to >> me. >> >> 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? >> >> why is >> x <

Re: [R] A point in a vector?

2009-09-30 Thread David Winsemius
On Sep 30, 2009, at 4:10 PM, Corrado wrote: Dear list, I have a strange requirement I have a vector, for example v<- c(0,0,0,0,1,2,4,6,8,8,8,8). I have a value,for example x<- 4.8. I would like to understand in which sub interval of v is x. In this case, v would be in the sub interval

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Duncan Murdoch
On 9/30/2009 3:59 PM, Ista Zahn wrote: For my own edification more than anything (I never took computer science): is a = seq(0.1,0.9,by=0.1) a <- as.character(a) a[3] == "0.3" [1] TRUE safe? No. Someone might be in a locale where the comma is used as the decimal separator. Take a look at

Re: [R] A point in a vector?

2009-09-30 Thread baptiste auguie
Hi, assuming v is sorted, try this, v[ findInterval(x,v)+0:1 ] see ?findInterval and perhaps ?cut HTH, baptiste 2009/9/30 Corrado : > Dear list, > > I have a strange requirement I have a vector, for example v<- > c(0,0,0,0,1,2,4,6,8,8,8,8). I have a value,for example x<- 4.8. > > I wou

Re: [R] A point in a vector?

2009-09-30 Thread Henrique Dallazuanna
Try this: v[findInterval(x, v) + 0:1] On Wed, Sep 30, 2009 at 5:10 PM, Corrado wrote: > Dear list, > > I have a strange requirement I have a vector, for example v<- > c(0,0,0,0,1,2,4,6,8,8,8,8). I have a value,for example x<- 4.8. > > I would like to understand in which sub interval of v is

[R] A point in a vector?

2009-09-30 Thread Corrado
Dear list, I have a strange requirement I have a vector, for example v<- c(0,0,0,0,1,2,4,6,8,8,8,8). I have a value,for example x<- 4.8. I would like to understand in which sub interval of v is x. In this case, v would be in the sub interval [4,6] that is in the subinterval starting from

[R] How can I use event.chart (Hmisc) to work with a single row of data?

2009-09-30 Thread George Chen
Hi Everyone, I would like to use event.chart to make an event chart based on a single row of data. It works fine for 2 or more rows, but when I set subset.r to 1 row or use a one row data.frame as the data, I get a dimensional error. Does anyone know a work around or a solution? Thanks, Geor

Re: [R] rcs fits in design package

2009-09-30 Thread Frank E Harrell Jr
Hayes, Rachel M wrote: Hi all, I have a vector of proportions (post_op_prw) such that >summary(amb$post_op_prw) Min. 1st Qu. MedianMean 3rd Qu.Max.NA's 0. 0. 0. 0.3985 0.9134 0.9962 1. summary(cut2(amb$post_op_prw,0.0001)) [0.

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Ista Zahn
For my own edification more than anything (I never took computer science): is > a = seq(0.1,0.9,by=0.1) > a <- as.character(a) > a[3] == "0.3" [1] TRUE safe? -Ista On Wed, Sep 30, 2009 at 3:46 PM, cls59 wrote: > > > Martin Batholdy wrote: >> >> hum, >> >> can you explain that a little more det

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Douglas Bates > Sent: Wednesday, September 30, 2009 12:42 PM > To: Dieter Menne > Cc: r-help@r-project.org > Subject: Re: [R] Condition to factor (easy to remember) > > On Wed, S

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread Peter Dalgaard
Douglas Bates wrote: On Wed, Sep 30, 2009 at 2:42 PM, Douglas Bates wrote: On Wed, Sep 30, 2009 at 2:43 AM, Dieter Menne wrote: Dear List, creating factors in a given non-default orders is notoriously difficult to explain in a course. Students love the ifelse construct given below most, but

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread hadley wickham
On Wed, Sep 30, 2009 at 2:32 PM, Ista Zahn wrote: > An extremely verbose, but (in my view) easy to understand approach is: > >> data.f <- data; data.f[which(data <= 10)] <- levs[1]; data.f[which(data > >> 10)] <- levs[2]; data.f <- factor(data.f) > All those which()s are unnecessary. And if you

Re: [R] Rounding error in seq(...)

2009-09-30 Thread cls59
Martin Batholdy wrote: > > hum, > > can you explain that a little more detailed? > Perhaps I miss the background knowledge - but it seems just absurd to > me. > > 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? > > Unfortunately this comes as an utter shock to many people wh

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread Douglas Bates
On Wed, Sep 30, 2009 at 2:42 PM, Douglas Bates wrote: > On Wed, Sep 30, 2009 at 2:43 AM, Dieter Menne > wrote: > >> Dear List, > >> creating factors in a given non-default orders is notoriously difficult to >> explain in a course. Students love the ifelse construct given below most, >> but I reme

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Gabor Grothendieck
See: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f On Wed, Sep 30, 2009 at 2:40 PM, Michael Knudsen wrote: > Hi, > > Today I was flabbergasted to see something that looks like a rounding > error in the very basic seq function in R. > >> a = seq

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread Douglas Bates
On Wed, Sep 30, 2009 at 2:43 AM, Dieter Menne wrote: > Dear List, > creating factors in a given non-default orders is notoriously difficult to > explain in a course. Students love the ifelse construct given below most, > but I remember some comment from Martin Mächler (?) that ifelse should be >

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Greg Snow
Read FAQ 7.31 (no you are not the first to notice this). In fact this is an issue with most computer programs, not unique or original to R. The FAQ has a link to a general computer science article that goes into great detail on this issue. Reading that article (or something similar) should be

Re: [R] Rounding error in seq(...)

2009-09-30 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Martin Batholdy > Sent: Wednesday, September 30, 2009 12:18 PM > To: Duncan Murdoch > Cc: r help > Subject: Re: [R] Rounding error in seq(...) > > hum, > > can you explain that

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread Ista Zahn
An extremely verbose, but (in my view) easy to understand approach is: > data.f <- data; data.f[which(data <= 10)] <- levs[1]; data.f[which(data > > 10)] <- levs[2]; data.f <- factor(data.f) -Ista On Wed, Sep 30, 2009 at 8:31 AM, Dieter Menne wrote: > > > > David Winsemius wrote: >> >> >>> # T

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Peter Dalgaard
Martin Batholdy wrote: hum, can you explain that a little more detailed? Perhaps I miss the background knowledge - but it seems just absurd to me. 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? why is x <- 0.1 + 0.1 +0.1 not equal to y <- 0.3 Remember that this is in BINARY ar

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Rolf Turner
On 1/10/2009, at 8:17 AM, Martin Batholdy wrote: hum, can you explain that a little more detailed? Perhaps I miss the background knowledge - but it seems just absurd to me. 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? why is x <- 0.1 + 0.1 +0.1 not equal to y <- 0.3 ? You

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Duncan Murdoch
On 9/30/2009 3:17 PM, Martin Batholdy wrote: hum, can you explain that a little more detailed? Perhaps I miss the background knowledge - but it seems just absurd to me. 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? R uses a binary floating point storage format, which is the

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Michael Knudsen
On Wed, Sep 30, 2009 at 8:40 PM, Michael Knudsen wrote: >> a = seq(0.1,0.9,by=0.1) >> a > [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 >> a[1] == 0.1 > [1] TRUE >> a[2] == 0.2 > [1] TRUE >> a[3] == 0.3 > [1] FALSE A friend of mine just pointed out a possible solution: > a=seq(0.1,0.9,by=0.1) > a = s

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Martin Batholdy
hum, can you explain that a little more detailed? Perhaps I miss the background knowledge - but it seems just absurd to me. 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? why is x <- 0.1 + 0.1 +0.1 not equal to y <- 0.3 ? Am 30.09.2009 um 21:04 schrieb Duncan Murdoch: On

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Duncan Murdoch
On 9/30/2009 2:50 PM, Michael Knudsen wrote: On Wed, Sep 30, 2009 at 8:44 PM, Duncan Murdoch wrote: Why? You asked for an increment of 1 in the second case (which is exactly represented in R), then divided by 10, so you'll get the same as 0.3 gives you. In the seq() case you asked for an inc

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Bert Gunter
Inline Below. Bert Gunter Nonclinical Biostatistics 467-7374 http://devo.gene.com/groups/devo/depts/ncb/home.shtml -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Michael Knudsen Sent: Wednesday, September 30, 2009 11:50 AM T

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Michael Knudsen
On Wed, Sep 30, 2009 at 8:44 PM, Duncan Murdoch wrote: > Why?  You asked for an increment of 1 in the second case (which is exactly > represented in R), then divided by 10, so you'll get the same as 0.3 gives > you.  In the seq() case you asked for an increment of a number close to but > not equa

Re: [R] Rounding error in seq(...)

2009-09-30 Thread Duncan Murdoch
On 9/30/2009 2:40 PM, Michael Knudsen wrote: Hi, Today I was flabbergasted to see something that looks like a rounding error in the very basic seq function in R. a = seq(0.1,0.9,by=0.1) a [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 a[1] == 0.1 [1] TRUE a[2] == 0.2 [1] TRUE a[3] == 0.3 [1] F

[R] Rounding error in seq(...)

2009-09-30 Thread Michael Knudsen
Hi, Today I was flabbergasted to see something that looks like a rounding error in the very basic seq function in R. > a = seq(0.1,0.9,by=0.1) > a [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 > a[1] == 0.1 [1] TRUE > a[2] == 0.2 [1] TRUE > a[3] == 0.3 [1] FALSE It turns out that the alternative > a

[R] change axis labels in persp

2009-09-30 Thread Bond, Stephen
I need to change one of the axis produced by persp by giving it labels for the tickmarks. The dimension has months (a factor) so the default decimals don't look good. The graph will finally become a pdf or emf file to be embedded in a document, so persp3d will not work, I think. Thank you all.

[R] start a script directly under windows

2009-09-30 Thread Martin Batholdy
hi, I try to configure my system so when I click on an .R script, the script is executed in the R GUI. I tried to open .R files with different .exe files that I found in the bin-folder of the R folder (R.exe, Rcmd.exe, Rgui.exe ...). But I always get the message; ARGUMENT 'c:\...\...'

Re: [R] how to visualize gini coefficient in each node in RF?

2009-09-30 Thread Liaw, Andy
You can try to hunt for it in the findbestsplit Fortran subroutine. It uses some thing that's equivalent (but easier to compute), not exactly identical. Breiman uses whatever computational shortcuts he could find in his code. Best, Andy From: Chrysant

Re: [R] bwplot scales in alphabetical order

2009-09-30 Thread Tim Clark
Peter, Thanks, that did it! Tim Tim Clark Department of Zoology University of Hawaii --- On Wed, 9/30/09, Peter Ehlers wrote: > From: Peter Ehlers > Subject: Re: [R] bwplot scales in alphabetical order > To: "Tim Clark" > Cc: r-help@r-project.org > Date: Wednesday, September 30, 2009, 2:

Re: [R] Scaling data

2009-09-30 Thread Henrique Dallazuanna
Try this: as.data.frame(sapply(sapply(DF, as.character, simplify = FALSE), strsplit, ";")) On Wed, Sep 30, 2009 at 1:18 PM, Dry, Jonathan R wrote: > Hello all > > I am manipulating some data and wish to expand/unmerge (i.e. do the opposite > of aggregate) rows in a data matrix based on the valu

[R] error in function AMER

2009-09-30 Thread Marilia Gmail
Dear all, I am trying to reproduce the example in the vignette "Using lme4 to fit Generalized Additive Mixed Models" with my dataset. But... > exemplo <- read.table(exemplo.dat,header=T) > mod <- amer(pasvig ~ -1 + harvf + tp(dias,by=harvf) + (1 | pac), data=exemplo) Erro em if (from == to)

[R] extract clusters from a hclusrt object without cutree

2009-09-30 Thread Larson, TR
Hello, I am generating an hclust() object and want to extract a subset of clusters WITHOUT using cutree(). I see this somehow involves using the $merge matrix to associate $labels into their appropriate clusters where clusters are split below a given $height threshold. Can anyone suggest a

Re: [R] Macro variable substitution

2009-09-30 Thread Greg Snow
There was a rather in depth article on doing macros in R in the R newsletter a couple of years ago. There is also the defmacro function in the gtools package based on that article that can make the job easier for many cases. See fortune(236) before taking the assign route. -- Gregory (Greg) L

Re: [R] Scaling data

2009-09-30 Thread Dry, Jonathan R
Hello all I am manipulating some data and wish to expand/unmerge (i.e. do the opposite of aggregate) rows in a data matrix based on the values in a particular column and a seperator, e.g. Col1Col2 n1;n2 6 ...separating by ";" becomes Col1Col2 num16 num26 Any ideas? Als

Re: [R] Managing random number generating, while using Condor parallel computing

2009-09-30 Thread Tal Galili
Thanks/Sorry Dirk. After some more searches I found: http://cran.r-project.org/web/packages/rsprng/index.html With a nice explenation about it here: http://wiki.r-project.org/rwiki/doku.php?id=packages:cran:rsprng Thanks/Sorry again, Tal -- My co

Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread joris meys
Hi Tal, The Wiki is a great idea! Alas I'm not that familiar with setting up Wikis. Any idea how to practically approach this? I'm currently a bit in a hurry, as I have to give the introduction next monday. I'll have a first, basic introduction ready by the end of the week, but I'll continue to d

Re: [R] Managing random number generating, while using Condor parallel computing

2009-09-30 Thread Dirk Eddelbuettel
On Wed, Sep 30, 2009 at 05:02:15PM +0200, Tal Galili wrote: > Might someone advice me about the following question: > How should I handle the RNG (random number generation) in the running of > parallel instances of R on different machines. > http://lmgtfy.com/?q=site%3Ar-project.org+parallel+RNG

Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread Tal Galili
Hi Joris. This is a wonderlful idea. I never came across such a document, but assembaling it sounds very good. What do you think about organizing the work of collecting (by R members and maybe your students) the error massages and code into some wiki page? (maybe on the R wiki) It could be a gre

[R] Openning multiple (writing) connections to a single file - possible ?

2009-09-30 Thread Tal Galili
Hello members of the R-help group. I am trying to run multiple occurences of R (through the distributed computing of Condor), and have all of them write their output into one source. It is similar (I believe), to running R several time on the same machine, and running on each instance of R the sa

Re: [R] Read header csv file

2009-09-30 Thread Don MacQueen
Lucas, Although I can't answer all your questions, I can give some suggestions. I will assume you know how many years of data you have. Suppose it is 2004 through 2009. In each input file, the columns must have names, so I will assume there is a column named "Month". I am not sure if you want

Re: [R] AsciiDoc and R

2009-09-30 Thread Steve Lianoglou
Hi, On Sep 30, 2009, at 7:57 AM, johannes rara wrote: I would like to learn AsciiDoc. Is there any good examples how to use AsciiDoc with R? I know that there is packages called ascii to do this, but it would be nice to see some examples how AsciiDoc works with R. Is there an AsciiDoc distribut

Re: [R] k-modes (Huang) in package klaR?

2009-09-30 Thread Uwe Ligges
Uwe Ligges wrote: Kim Vanselow wrote: Dear R-Users and Developers, I want to calculate something like k-means clustering, but with ordinal data (Braun-Blanquet) to combine this classification technique with a NMDS-Ordination. I found an algorithm especially developed for categorical data:

[R] Managing random number generating, while using Condor parallel computing

2009-09-30 Thread Tal Galili
Hello all, Recently I started playing with running R scripts on the Condor system in my institute. (For more on this, have a look at: Running Long R Jobs with Condor DAG by Xianhong Xie link: http://cran.r-project.org/doc/Rnews/Rnews_2005-2.pdf ) Might someone advice me about the following quest

[R] Re cursive regression

2009-09-30 Thread manta
Hi there, I'm in desperate need to figure out how to solve this issue. I need to estimate a recursive model for a time series data of asset returns. The dependent variable is the asset return and then I have a set of k variables, a lagged value of the dependent variable (plus an intercept) as regr

Re: [R] Optim(...) estimate of stDev far too low

2009-09-30 Thread Peter Ehlers
Luis, I don't think that you want log() in fn. Try it with pdf=-sum(dnorm(x, media, st, TRUE)) (and note the x, rather than xx). -Peter Ehlers Luis Ridao Cruz wrote: R-help, I'm just trying to find the ML (maximum likelihood) estimates of the mean and standard deviation of a set of observ

Re: [R] k-modes (Huang) in package klaR?

2009-09-30 Thread Uwe Ligges
Kim Vanselow wrote: Dear R-Users and Developers, I want to calculate something like k-means clustering, but with ordinal data (Braun-Blanquet) to combine this classification technique with a NMDS-Ordination. I found an algorithm especially developed for categorical data: k-modes (Huang 1998)

Re: [R] Condition to factor (easy to remember)

2009-09-30 Thread Gabor Grothendieck
1. A common way of doing this is cut: > cut(data, c(-Inf, 10, Inf), lab = levs, right = TRUE) [1] Pre Pre Pre Post Post Levels: Pre Post We don't actually need right=TRUE as its the default but if you omit it then it can be hard to remember whether the right end of intervals are included

Re: [R] How to calculate KMO?

2009-09-30 Thread G. Jay Kerns
Dear Moumita, On Wed, Sep 30, 2009 at 10:03 AM, Moumita Das wrote: > Hi All, > How do i calculate KMO for a dataset? > > *Dataset:-* > >    m1 m2 m3 m4 m5 m6 m7 m8 > 1       2 20 20  2  1  4 14 12 > 2       9 16  3  5  2  5  5 15 > 3      18 18 18 13 17  9  2  4 > 4       7

Re: [R] R/PL : cannot build PL/R

2009-09-30 Thread David Winsemius
Wrong mailing list. Suggest reposting on R-SIG-Mac but only after re- reading the Posting Guide. On Sep 30, 2009, at 10:03 AM, Bunny, lautloscrew.com wrote: Dear all, i am trying to get PL / R running on Mac OS X. PostgreSQL is already running and connects just smoothly to R via DBI / RPos

Re: [R] Change directory to implement same programes

2009-09-30 Thread Paul Hiemstra
Hi, You can use list.files() of Sys.glob to get a listing of all the files in a certain directory, e.g. to get all ".R" files and source them: for(f in Sys.glob("C:/Documents and Settings/lma/*.R")) source(f) cheers, Paul Tammy Ma wrote: How do I source a bunch of files in different directo

[R] k-modes (Huang) in package klaR?

2009-09-30 Thread Kim Vanselow
Dear R-Users and Developers, I want to calculate something like k-means clustering, but with ordinal data (Braun-Blanquet) to combine this classification technique with a NMDS-Ordination. I found an algorithm especially developed for categorical data: k-modes (Huang 1998). Is there a function im

[R] How to calculate KMO?

2009-09-30 Thread Moumita Das
Hi All, How do i calculate KMO for a dataset? *Dataset:-* m1 m2 m3 m4 m5 m6 m7 m8 1 2 20 20 2 1 4 14 12 2 9 16 3 5 2 5 5 15 3 18 18 18 13 17 9 2 4 4 7 7 2 12 2 11 11 11 5 7 8 5 19 5 2 20 18 6 7 4 7 4 7 9 3 3 7

[R] R/PL : cannot build PL/R

2009-09-30 Thread Bunny, lautloscrew.com
Dear all, i am trying to get PL / R running on Mac OS X. PostgreSQL is already running and connects just smoothly to R via DBI / RPostgreSQL . Despite finding a couple of posts of people with the same problem (i.e. error message) i could not get it done for me. I use Mac OS 10.5.7 and Postgr

[R] rcs fits in design package

2009-09-30 Thread Hayes, Rachel M
Hi all, I have a vector of proportions (post_op_prw) such that >summary(amb$post_op_prw) Min. 1st Qu. MedianMean 3rd Qu.Max.NA's 0. 0. 0. 0.3985 0.9134 0.9962 1. > summary(cut2(amb$post_op_prw,0.0001)) [0.,0.0001) [0.0001,0.9962]

[R] Optim(...) estimate of stDev far too low

2009-09-30 Thread Luis Ridao Cruz
R-help, I'm just trying to find the ML (maximum likelihood) estimates of the mean and standard deviation of a set of observations: >xx=c(2.5,3.5,4,6,6.5,7.5) fn<-function(params,x=xx) { media<-params[1] st <-params[2] pdf=-sum(dnorm(log(xx),log(media),st,TRUE)) return(pdf) } optim(c(mu,stdev

Re: [R] How to compile R with command completion?

2009-09-30 Thread Martin Morgan
Peng Yu wrote: > On Tue, Sep 29, 2009 at 10:39 PM, Martin Morgan wrote: >> Peng Yu wrote: >>> On Tue, Sep 29, 2009 at 3:47 PM, Tobias Verbeke >>> wrote: Peng Yu wrote: > I want to compile R with command completion. But I don't find such an > option in configure. Can somebody let

[R] Read header csv file

2009-09-30 Thread Lucas Sevilla García
Hi R community, First of all, I want to thank everybody to share their time solving R questions, You are great. Ok, for my questions, I've been looking for a solutions by myself, in forums but I'm just a little bit desesperate so I hope somebody can help me. I have built a code to read files

[R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread joris meys
Dear all, The background: as I didn't find a manual specifically directed towards debugging R code on a novice and intermediate level, I decided to write a student manual myself. I have some basic examples of errors (the classical interpunction and writing mistakes), but have the feeling that I mi

Re: [R] AsciiDoc and R

2009-09-30 Thread David Winsemius
On Sep 30, 2009, at 7:57 AM, johannes rara wrote: I would like to learn AsciiDoc. Is there any good examples how to use AsciiDoc with R? I know that there is packages called ascii to do this, but it would be nice to see some examples how AsciiDoc works with R. Is there an AsciiDoc distribution

Re: [R] bwplot scales in alphabetical order

2009-09-30 Thread Peter Ehlers
Tim, Add the argument as.table=TRUE to your call: bwplot(y~x|id, horizontal=FALSE, as.table=TRUE) Peter Ehlers Tim Clark wrote: Dear List, I know this has been covered before, but I don't seem to be able to get it right. I am constructing a boxplot in lattice and can't get the scales in t

  1   2   >