Re: [R] Help with makeClusters for Snow

2009-12-23 Thread Patrick Connolly
On Tue, 22-Dec-2009 at 09:38PM -0800, SVM wrote: |> [] |> |> On top of that, I couldn't find the MPI/NWS/PVM libraries for windows. I've |> looked for them all over the internet but I can't seem to find any resources |> for it. I even tried installing them directly from R, (i.e. |> install.

Re: [R] use of lm() and poly()

2009-12-23 Thread Benoit Boulinguiez
I guess, I would somehow feel interested that someone out of my field is trying to use the tool I daily manipulate to widen his knowledge or dig further in his own field, though I might at some point recommend him to get back to the basis for some concrete concepts that he wouldn't understand sha

Re: [R] Help with makeClusters for Snow

2009-12-23 Thread Prof Brian Ripley
See http://cran.r-project.org/bin/windows/contrib/2.10/@ReadMe for what packages are not made available as binary packages for Windows. You can install packages Rmpi and rpvm from the sources (but not easily) if you have the requisite software installed. See e.g. http://www.stats.uwo.ca/facul

[R] String question

2009-12-23 Thread Knut Krueger
Hi to all I need a string like temp <- paste("m1","m2","m3",sep=",") But i must know how many items are in the string,afterwards the other option would be to use a vector temp <- c("m1","m2","m3") No problem to get the count of items but I must get afterwards the string "m1,m2,m3" No problem to

[R] COnfidence intervals for estimates of linear model

2009-12-23 Thread Daniel Brewer
Hello, I would like to calculate the 95% confidence intervals for the estimates of a linear model and I just wanted to check that I am doing it correct. Is it just: Estimate + 1.95996*Std.Error to Estimate - 1.95996*Std.Error or is there another approach that doesn't assume a normal distrbution

Re: [R] String question

2009-12-23 Thread baptiste auguie
Will this do? temp <- paste("m", 1:3, sep="",collapse=",") HTH, baptiste 2009/12/23 Knut Krueger : > Hi to all > > I need a string like > temp <- paste("m1","m2","m3",sep=",") > But i must know how many items are in the string,afterwards > the other option would be to use a vector > temp <- c("

Re: [R] String question

2009-12-23 Thread Knut Krueger
Will this do? temp <- paste("m", 1:3, sep="",collapse=",") Unfortunately not, because I explained the example not detailed enough. The string could have different Items, like November, December, Monday, Tuesday, Daylight and so on Therefore I must count the Items of the string separated b

Re: [R] String question

2009-12-23 Thread Jim Lemon
On 12/23/2009 09:21 PM, Knut Krueger wrote: Hi to all I need a string like temp <- paste("m1","m2","m3",sep=",") But i must know how many items are in the string,afterwards the other option would be to use a vector temp <- c("m1","m2","m3") No problem to get the count of items but I must get aft

Re: [R] String question

2009-12-23 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 23.12.2009 11:46:31: > > > Will this do? > > > > temp <- paste("m", 1:3, sep="",collapse=",") > > > > > Unfortunately not, because I explained the example not detailed enough. > The string could have different Items, like November, December, Monday,

Re: [R] COnfidence intervals for estimates of linear model

2009-12-23 Thread ONKELINX, Thierry
Dear Daniel, Have a look at ?predict.lm The interval argument gives you the information that you need. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverst

Re: [R] String question

2009-12-23 Thread Knut Krueger
Jim Lemon schrieb: Not as easy as I thought it would be, but: mlist<-as.list(paste("m",1:sample(5:10,1),sep="")) do.call("paste",c(mlist,sep=",")) Hi Jim, yes it works :-) temp <- c("November", "December","Monday","Tuesday") length(temp) #getting the length of the vector string1=do.call("p

[R] Mean, median and other moments

2009-12-23 Thread Amelia Livington
Hi! Suppose I have a dataset as follows pd = c(10,7,10,11,7,11,7,6,8,3,12,7,7,10,10) I wish to calculate the mean, standard deviation, median, skewness and kurtosis i.e. regular standard statistical measures. average = mean(pd) stdev    = sd(pd) median = median(pd) skew    = skewness(pd) kurt 

Re: [R] String question

2009-12-23 Thread Ted Harding
On 23-Dec-09 11:08:02, Knut Krueger wrote: > Jim Lemon schrieb: >> Not as easy as I thought it would be, but: >> >> mlist<-as.list(paste("m",1:sample(5:10,1),sep="")) >> do.call("paste",c(mlist,sep=",")) > > Hi Jim, > yes it works :-) > > temp <- c("November", "December","Monday","Tuesday") > le

Re: [R] String question

2009-12-23 Thread baptiste auguie
Isn't paste doing exactly this? temp <- c("November", "December","Monday","Tuesday") paste(temp, collapse=",") # "November,December,Monday,Tuesday" HTH, baptiste 2009/12/23 Ted Harding : > On 23-Dec-09 11:08:02, Knut Krueger wrote: >> Jim Lemon schrieb: >>> Not as easy as I thought it would be

Re: [R] Mean, median and other moments

2009-12-23 Thread Jim Lemon
On 12/23/2009 10:24 PM, Amelia Livington wrote: Hi! Suppose I have a dataset as follows pd = c(10,7,10,11,7,11,7,6,8,3,12,7,7,10,10) I wish to calculate the mean, standard deviation, median, skewness and kurtosis i.e. regular standard statistical measures. average = mean(pd) stdev= sd(pd

Re: [R] String question

2009-12-23 Thread Ted Harding
On 23-Dec-09 11:40:12, baptiste auguie wrote: > Isn't paste doing exactly this? > > temp <- c("November", "December","Monday","Tuesday") > paste(temp, collapse=",") ># "November,December,Monday,Tuesday" > > HTH, > baptiste Yes, spot-on! I got involved in the confusion resulting from the use of "

Re: [R] String question

2009-12-23 Thread Knut Krueger
Hi Baptiste, Isn't paste doing exactly this? yes indeed - surprising temp <- c("November", "December","Monday","Tuesday") paste(temp, collapse=",") paste(temp, sep=",") I tried to use sep :-( Arguments |...| one or more *R* objects, to be converted to character vectors. |sep|

Re: [R] String question

2009-12-23 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 23.12.2009 12:08:02: > Jim Lemon schrieb: > > > > Not as easy as I thought it would be, but: > > > > mlist<-as.list(paste("m",1:sample(5:10,1),sep="")) > > do.call("paste",c(mlist,sep=",")) > > > > Hi Jim, > yes it works :-) > > temp <- c("November",

[R] how to create normal qqplot with the 95% confidence interval

2009-12-23 Thread CJ Rubio
hi everyone! season's greetings! is there any way that i can create a normal qqplot showing, aside from the qqline, the 95% confidence limits? thank you very much.. happy holidays! -- View this message in context: http://n4.nabble.com/how-to-create-normal-qqplot-with-the-95-confidence-interva

Re: [R] Playing with rgl: a Youtube video

2009-12-23 Thread Rainer M Krug
On Tue, Dec 22, 2009 at 8:10 PM, Duncan Murdoch wrote: > On 22/12/2009 12:49 PM, Mark Knecht wrote: > >> On Mon, Dec 21, 2009 at 4:42 AM, Duncan Murdoch >> wrote: >> > I've just posted a demo made with the rgl package to Youtube, visible >> here: >> > http://www.youtube.com/watch?v=prdZWQD7L5c >

Re: [R] Playing with rgl: a Youtube video

2009-12-23 Thread Duncan Murdoch
On 23/12/2009 7:13 AM, Rainer M Krug wrote: On Tue, Dec 22, 2009 at 8:10 PM, Duncan Murdoch > wrote: On 22/12/2009 12:49 PM, Mark Knecht wrote: On Mon, Dec 21, 2009 at 4:42 AM, Duncan Murdoch mailto:murd...@stats.uwo.ca>> wrote: > I'v

Re: [R] Mean, median and other moments

2009-12-23 Thread John Sorkin
An answer to your second question: Random number generated by your computer are not truly random, they are pseduo random meaning that if your computer is it the exact same state at two movements, it will produce the same string of "random" numbers. Set.seed effectively puts your machine in a giv

Re: [R] String question

2009-12-23 Thread Gustaf Rydevik
On Wed, Dec 23, 2009 at 11:21 AM, Knut Krueger wrote: > Hi to all > > I need a string like > temp <- paste("m1","m2","m3",sep=",") > But i must know how many items are in the string,afterwards > the other option would be to use a vector > temp <- c("m1","m2","m3") > No problem to get the count of

Re: [R] how to create normal qqplot with the 95% confidence interval

2009-12-23 Thread brestat
Rubio, Look at library(fBasics) the function qqnormPlot(). Below an example: qqnormPlot(rnorm(100)) Best's Walmes Zeviani, Brasil. CJ Rubio wrote: > > hi everyone! > > season's greetings! > > is there any way that i can create a normal qqplot showing, aside from the > qqline, the 95%

Re: [R] by function ??

2009-12-23 Thread Matthew Dowle
You asked how to 'create confidence intervals around the median'. Since intervals was plural and you asked within the context of 'by' LEAID then that made sense and I guessed you meant 'confidence interval around the median ratio of each group'. I already assumed you had posted a small subset o

Re: [R] Problem with expand.grid

2009-12-23 Thread Martin Maechler
> "KJ" == Keith Jewell > on Tue, 22 Dec 2009 16:02:49 - writes: KJ> Hi All, KJ> This example code KJ> KJ> dDF <- structure(list(y = c(4.75587, 4.8451, 5.04139, 4.85733, 5.20412, KJ> 5.92428, 5.69897, 4.78958, 4, 4), t = c(0, 48, 144, 192, 240,

[R] panel data quantile regression?

2009-12-23 Thread Junyu
Hi All, About the quantile regression of panel data, what are the difference about Koenker’ s method of “Quantile Regression for Longitudinal Data” and added a dummy variable to represent individual effect? For example, if I have a five year and three country panel data set , I added following du

Re: [R] [R-SIG-Finance] Fitting ACD model

2009-12-23 Thread Brian G. Peterson
R_help Help wrote: Hi - I'm wondering if there is any existing package in R that facilitates ACD (Autoregressive Conditional Duration) model. I googled. Apparently there were packages dynamo and fACD. However, they do not seem to be on any CRAN. I am wondering if anyone could point me to an exist

Re: [R] Cohen's kappa, unequal score ranges

2009-12-23 Thread Julia Myatt
Hi Jim, Thanks for your help, I tried that function which did result in an output. The only issue now is my summary looks like this: Length Class Mode table 16 -none- numeric kappa 1 -none- numeric Not what I was expecting, this was the case when my data was in the following

Re: [R] Cohen's kappa, unequal score ranges

2009-12-23 Thread Scot W. McNary
Julia, I was able to get your toy data to work, using the n * 2 data frame, which is what ckappa takes as input. > chk <- data.frame(matrix(c( +3, 4, +4, 3, +2, 1, +2, 1, +5, 1, +2, 4), byrow = TRUE, ncol = 2)) > chk X1 X2 1 3 4 2 4 3 3 2 1 4 2 1

[R] Rgraphviz on mac 10.6.2

2009-12-23 Thread kulwinder banipal
Rgraphviz Install works fine (http://www.bioconductor.org/packages/release/bioc/html/Rgraphviz.html) Latest version of graphviz is installed as well however I get following error when loading Rgraphviz (on Mac 10.6.2) Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared l

[R] Fw: R and minitab : Censored Kendall

2009-12-23 Thread Dennis J Low
Dennis J. Low Hydrologist U.S.Geological Survey - WRD 215 Limekiln Road New Cumberland, Pennsylvania 17070 (717) 730-6959 I have been running censored Kendall on Minitab, and more recently on R because of the large data size (about 4,000 samples). I have begun to compare the two results an

[R] animated R plots

2009-12-23 Thread Zd Gibbs
 Hi, I want to be able to save the following animated plot as a flash.  There are ultimately 6 plots, but when I run this and save it as a flash all I get is the last (6th) plot, not all six different plots in order by "year".  I’ve tested the flash commands and they work; so it has to be m

Re: [R] Cohen's kappa, unequal score ranges

2009-12-23 Thread Julia Myatt
Hi, Thank you so much for your help with this, I can get it to work now! Just one more question, my actual data will involve hundreds of rows in my two columns, is there anyway to import my data frame via read.table, and then convert this into the matrix ready for ckappa, rather than typing th

Re: [R] Cohen's kappa, unequal score ranges

2009-12-23 Thread Julia Myatt
Hi, Sorry, it's ok I've figured it out using the as.matrix function! Cheers, Julia. From: Scot W. McNary [smcn...@charm.net] Sent: 23 December 2009 15:41 To: Julia Myatt Cc: r-help@r-project.org Subject: Re: [R] Cohen's kappa, unequal score ranges Ju

Re: [R] trouble with model.tables SE means

2009-12-23 Thread Jon Prince
David and Peter, Thanks so much for all of your help, I think I understand R much better as a result. Omitting the Error() term in my aov does indeed allow me to get SE means, so I guess that was the issue. I suppose I can go back and calculate the SE values for each p*t entry (averaged across

Re: [R] COnfidence intervals for estimates of linear model

2009-12-23 Thread Greg Snow
Using 1.95996 is appropriate when you know the population standard deviation or the sample size is approximately infinity. Otherwise you should use the t-distribution, the qt function is useful for that. Or if you want intervals for the coefficients look at the confint function, if you want in

Re: [R] how to create normal qqplot with the 95% confidence interval

2009-12-23 Thread Greg Snow
What type of confidence limits do you want? And more importantly, why do you want them? -- 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...@

Re: [R] trouble with model.tables SE means

2009-12-23 Thread David Winsemius
On Dec 23, 2009, at 11:45 AM, Jon Prince wrote: David and Peter, Thanks so much for all of your help, I think I understand R much better as a result. Omitting the Error() term in my aov does indeed allow me to get SE means, so I guess that was the issue. I suppose I can go back and calcu

[R] iid.test

2009-12-23 Thread Gavriel & Esti Zoladz
I downloaded the iid.test, but I can't run it. I get the following message: Error: could not find function "iid.test" Where am I supposed to save this package in order that it works? Thanks, EZ [[alternative HTML version deleted]] __ R-help@r-p

Re: [R] iid.test

2009-12-23 Thread Brian G. Peterson
Gavriel & Esti Zoladz wrote: I downloaded the iid.test, but I can't run it. I get the following message: Error: could not find function "iid.test" Where am I supposed to save this package in order that it works? ?install.packages __ R-help@r-proje

Re: [R] Problem with "Cannot compute correct p-values with ties"

2009-12-23 Thread Greg Snow
Adding random noise to data in order to avoid a warning is like removing the batteries from a smoke detector to silence it rather than investigating the what is causing the alarm to go off. If the function is giving a warning it is best to investigate why, it is possible that you can ignore the

Re: [R] Rgraphviz on mac 10.6.2

2009-12-23 Thread Martin Morgan
Please report issues with Bioconductor packages on the Bioconductor mailing list. http://bioconductor.org/docs/mailList.html The top of the page you mention has Users on all platforms must install graphviz; see the README file, available in the source distribution of this file, for details

Re: [R] ?OT: Probabilistic Simulation

2009-12-23 Thread Greg Snow
A simple starting place for your simulation could be something like: peg <- rnorm(10, 1.8, 0.01) hole <- rnorm(10, 1.8, 0.02) fits <- peg < hole table(fits) mean(fits) Then you would need to add in the effects of temperature, drill type, etc. to account for the different pieces. If yo

Re: [R] iid.test

2009-12-23 Thread Brian G. Peterson
You're clearly in over your head. install.packages("iid.test") will install the package. you may replace this with any "packagename" ?functionname in R on any 'functionname' will bring up the Fine Manual which you are to then use to Read the Fine Manual (RTFM) before posting questions to

Re: [R] iid.test

2009-12-23 Thread William Dunlap
After the install.packages("iid.test") you will need to call library("iid.test") The former you do once and the latter you do each time you restart R. install.packages() puts the package onto your computer's disk and library() loads it into an R session. Bill Dunlap Spotfire, TIBCO Software

Re: [R] Problem with "Cannot compute correct p-values with ties"

2009-12-23 Thread Ted Harding
The issue of ties in the Mann-Whitney test needs some thought. The distribution function of the Mann-Whitney test is derived on the assumption that (in effect) the data are continuous variables so that (theoretically) there should be no ties. Whn ties occur, then this assumjtion has failed. 1. If

[R] colorkey: Placement of labels and ticks in center positions of "color key" of a levelplot ?

2009-12-23 Thread Zia Ahmed
Is it possible to place labels and ticks in center positions of each cut of "color key" of levelplot ? Here is my code that I used to produce a map (see attachment). The two labels and ticks were placed in left sides each segment of the color key. I want to show them center of the cut

[R] loading data into ZOO

2009-12-23 Thread Stephen J. Barr
Hello, I have a simple question. I am trying to load data into a zoo object. I have the data in CSV format as follows SYMBOL DATE TIME PRICE XXMMDD HH:MM:SS n.nn and there are multiple symbols in this one data frame. My question is, do I need to mer

Re: [R] loading data into ZOO

2009-12-23 Thread Brian G. Peterson
Stephen J. Barr wrote: I have a simple question. I am trying to load data into a zoo object. I have the data in CSV format as follows SYMBOL DATE TIME PRICE XXMMDD HH:MM:SS n.nn and there are multiple symbols in this one data frame. My question is,

[R] Unwanted association between a function and a namespace

2009-12-23 Thread p_connolly
I can't understand how the plyr package is turning up here: > sessionInfo() R version 2.10.1 (2009-12-14) i686-pc-linux-gnu locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 [7] LC_P

Re: [R] 3-D barplot

2009-12-23 Thread Terry Therneau
> I would like to create a 3 dimensional barplot of 16 odds ratios that > demonstrate an interaction between two variables (CD14 and CD23). odds<- matrix( c(1, 0.61, 2.1, 6.1, 1.5 ,1.3, 3.5, 11.9, 1.8 ,1.5, 4.4, 10.6,

[R] Resent: colorkey: Placement of labels and ticks in center positions of "color key" of a levelplot ?

2009-12-23 Thread Zia Ahmed
Is it possible to place labels and ticks in center positions of each cut of "color key" of levelplot ? Here is my code that I used to produce a map (see attachment). The two labels and ticks were placed in left sides each segment of the color key. I want to show them center of the cut

Re: [R] Unwanted association between a function and a namespace

2009-12-23 Thread Duncan Murdoch
On 23/12/2009 3:36 PM, p_conno...@slingshot.co.nz wrote: I can't understand how the plyr package is turning up here: sessionInfo() R version 2.10.1 (2009-12-14) i686-pc-linux-gnu locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8 [5] L

Re: [R] loading data into ZOO

2009-12-23 Thread Gabor Grothendieck
read.zoo does not allow multiple index fields but you can read it in using read.table, combine the index fields and use read.zoo to read the resulting data frame. Note that read.zoo can read data frames, not just files. Also note the aggregate= argument on read.zoo allows direct handling of situat

[R] prcomp : plotting only explanatory axis arrows

2009-12-23 Thread milton ruser
Dear all, I have a very large dataset (1712351 , 20) and would like to plot only the arrows that represent the contribution of each variables. On the sample below I woild like to plot only the explanatory variables (Murder, Assault..) and not the sites. prcomp(USArrests) # inappropriate prcomp(U

Re: [R] as.Date question

2009-12-23 Thread Marek Janad
Look at documentation ?as.Date as.Date first represents time in UTC, what gives: as.POSIXlt(zzz1, tz="UTC") HTH 2009/12/20 MAL : > All! > > This piece of code: > > zzz1 <- as.POSIXct("1999-03-18", tz="CET") > zzz2 <- as.POSIXlt("1999-03-18", tz="CET") > zzz1 == zzz2 > as.Date(zzz1) > as.Date(z

Re: [R] Unwanted association between a function and a namespace

2009-12-23 Thread p_connolly
Quoting Duncan Murdoch : On 23/12/2009 3:36 PM, p_conno...@slingshot.co.nz wrote: I can't understand how the plyr package is turning up here: [...] I can understand that lattice would be using grid without having that package loaded, but I can't understand how plyr got there. One way fo

[R] ensemble for machine learning

2009-12-23 Thread Jack Lewis
Could someone give a better example of how to extend [1] to create an ensemble from multiple machine learning classifiers to improve the classification performance? [1] http://heuristically.wordpress.com/2009/12/23/compare-performance-machine-learning-classifiers-r/ Berry and Linoff ("Maste

Re: [R] Unwanted association between a function and a namespace

2009-12-23 Thread Duncan Murdoch
On 23/12/2009 6:08 PM, p_conno...@slingshot.co.nz wrote: Quoting Duncan Murdoch : On 23/12/2009 3:36 PM, p_conno...@slingshot.co.nz wrote: I can't understand how the plyr package is turning up here: [...] I can understand that lattice would be using grid without having that package loaded

Re: [R] Unwanted association between a function and a namespace

2009-12-23 Thread Patrick Connolly
Tnanks for the help, Duncan. On Wed, 23-Dec-2009 at 06:57PM -0500, Duncan Murdoch wrote: > On 23/12/2009 6:08 PM, p_conno...@slingshot.co.nz wrote: [...] >> I thought I'd found a fix when I specifically loaded plyr and then >> detached it with the unload = TRUE argument. sessionInfo() seems to

[R] Column naming issues using read.table

2009-12-23 Thread arthurbeer01
Hi, this is my first post so please be gentle. I quite new to R and using it for my biology degree. My problem is. Im trying to import data from a .csv file using the read.table command. The .csv file header starts on row 2 but is contained in column 1, i have 600 data files and for future ease w

Re: [R] Column naming issues using read.table

2009-12-23 Thread jim holtman
This reads in your posted data: > x <- read.table(textConnection("Samplerate = 2 samps/sec + Nr Cnt1X Cnt1Y Cnt2X Cnt2Y sec100 hour + 153 84 43 2 22 12 + 290 155 74 0 72 12 + 390 155 74 0 121

[R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Francesco Napolitano
Hi all, I'm trying to learn R after years of Matlab's experience. Here is an issue I couldn't solve today. Consider the following piece of code (written by memory): for(i in 1:n){ submat <- data[1:i,] C <- colSums(submat) } The problem is that at the first iteration, data[1:1,] redu

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Benilton Carvalho
replace data[1:i,] by data[1:i,drop=FALSE]. b On Dec 24, 2009, at 12:46 AM, Francesco Napolitano wrote: > Hi all, > > I'm trying to learn R after years of Matlab's experience. Here is an > issue I couldn't solve today. > > Consider the following piece of code (written by memory): > > for(i in

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Henrik Bengtsson
Use submat <- data[1:i,, drop=FALSE] /H On Wed, Dec 23, 2009 at 6:46 PM, Francesco Napolitano wrote: > Hi all, > > I'm trying to learn R after years of Matlab's experience. Here is an > issue I couldn't solve today. > > Consider the following piece of code (written by memory): > > for(i in 1

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread David Winsemius
On Dec 23, 2009, at 9:46 PM, Francesco Napolitano wrote: Hi all, I'm trying to learn R after years of Matlab's experience. Here is an issue I couldn't solve today. Consider the following piece of code (written by memory): for(i in 1:n){ submat <- data[1:i,] C <- colSums(submat) }

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Francesco Napolitano
So I see from help([[) what's happening: thank you very much. Il giorno mer, 23/12/2009 alle 18.59 -0800, Henrik Bengtsson ha scritto: > Use > >submat <- data[1:i,, drop=FALSE] > > /H > > On Wed, Dec 23, 2009 at 6:46 PM, Francesco Napolitano > wrote: > > Hi all, > > > > I'm trying to learn

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Francesco Napolitano
Il giorno mer, 23/12/2009 alle 22.14 -0500, David Winsemius ha scritto: > > for(i in 1:n){ > >submat <- data[1:i,] > >C <- colSums(submat) > >} > > > In R the loop is not necessary, even confusing as you are demonstrating: > > > mat <- matrix(rnorm(100), nrow=10) # 10 x 10 > > colSu

Re: [R] Resent: colorkey: Placement of labels and ticks in center positions of "color key" of a levelplot ?

2009-12-23 Thread Deepayan Sarkar
On Wed, Dec 23, 2009 at 12:54 PM, Zia Ahmed wrote: > Is it possible to place  labels  and ticks  in center positions of each cut > of  "color key" of levelplot ? Here is my code that I used to produce a map > (see attachment). The two  labels and ticks  were  placed  in  left  sides > each segment

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Brian G. Peterson
Francesco Napolitano wrote: Il giorno mer, 23/12/2009 alle 22.14 -0500, David Winsemius ha scritto: for(i in 1:n){ submat <- data[1:i,] C <- colSums(submat) } In R the loop is not necessary, even confusing as you are demonstrating: > mat <- matrix(rnorm(100), nrow=10) # 10 x 10 >

Re: [R] Newbie: colSums() compared with Matlab's sum()

2009-12-23 Thread Francesco Napolitano
On Thu, Dec 24, 2009 at 4:38 AM, Brian G. Peterson wrote: > > data <- matrix(rnorm(100), nrow=10) # 10 x 10 > apply(data,2,cumsum) > > if you want to store the partial sums > (that was an important bit of information you left out) > > Thank you for your suggestion. However, I gave information onl

Re: [R] how to create normal qqplot with the 95% confidence interval

2009-12-23 Thread CJ Rubio
thank you very much! this is just what my professor was asking for. Walmes Zeviani wrote: > > Rubio, > > > Look at library(fBasics) the function qqnormPlot(). Below an example: > > qqnormPlot(rnorm(100)) > > Best's > > Walmes Zeviani, Brasil. > > > > > > CJ Rubio wrote: >> >> hi eve

Re: [R] animated R plots

2009-12-23 Thread Yihui Xie
Hi, saveSWF() (and other save*() functions) will only record the plots by *high-level* plotting commands, and this is because most graphics devices can only record high-level plots by default. In your animation code, there is actually one plot generated, and the rest are produced by low-level plot

[R] What is ".Machine$double.eps"?

2009-12-23 Thread Saji Ren
Hello,everyone: I met this notation when I read the original code of function "quantile".There is one sentence as below: >eps <- 100 * .Machine$double.eps when I input ".Machine$double.eps" in R, it returns "[1] 2.220446e-16". Can anyone show me the exact meaning about that? thanks Saji from S

[R] Innocentive challenge

2009-12-23 Thread Jim Lemon
Hi all, Apologies if this is inappropriate, but I found this challenge on Innocentive: http://click.email.innocentive.com/?ju=fe5c1c777563037c7717&ls=fe0417717764067b711d7075&m=fef91270706102&l=fe8815797c62037d72&s=fe1e16787c640d797c1371&jb=ffcf14&t= Algorithms Optimized for Efficiency and Spee

Re: [R] What is ".Machine$double.eps"?

2009-12-23 Thread Berend Hasselman
Saji Ren wrote: > > Hello,everyone: > > I met this notation when I read the original code of function > "quantile".There is one sentence as below: >>eps <- 100 * .Machine$double.eps > > when I input ".Machine$double.eps" in R, it returns "[1] 2.220446e-16". > Can anyone show me the exact meani

Re: [R] Help with makeClusters for Snow

2009-12-23 Thread SVM
Hi, Would it help if I was to install Ubuntu, and try loading it there? I've seen quite a few comments about running rmpi on Linux, through google search. Do you know if I have to install any specific packages/libraries before trying to install RMpi? I was recommended to install LAM 7.1 before