[R] Correlation between two multi-dimensional matrices

Hi, Are there any methods for computing the correlation between two multi-dimensional matrices ? Will transforming the matrices into vectors and applying pearson be fine ? Any blind spots that I should be aware ? Thanks Stanley [[alternative HTML version deleted]] __

[R] preallocating matrices and rda read-back objects

I've read in Phil Spector's new book that it's a good idea to preallocate a big matrix, like u <- matrix(0,nrow,ncol) # (1) Now, I read contents of a huge matrix from a Fortran binary dump. u <- readBin(con,what="double",n=nrow*ncol) # (2) If I do (1) and then (2), u is a vector, obviously it

Re: [R] R in a cluster environment?

On Tue, 8 Apr 2008, Tornes, Ivan E wrote: I'm a cluster administrator and one of my users has requested that I install R on our system. Can R be used in a cluster computing environment that utilizes batch processing, i.e. no graphical output during processing? Yes, easily. Also, can R be p

Re: [R] question about nlminb

Have you considered optimizing over x1 = x[1:(length(x)-1]? You could feed a wrapper function 'f2(x1, ...)' that computes xFull = c(x1, 1-sum(x1)) and feeds that to your 'fn'. If this makes sense, great. Else, if my answer is not useful, be so kind as to PLEASE do read the posti

[R] R in a cluster environment?

I'm a cluster administrator and one of my users has requested that I install R on our system. Can R be used in a cluster computing environment that utilizes batch processing, i.e. no graphical output during processing? Also, can R be parallelized to run on many nodes/processors simultaneously?

[R] problem in ls.ranking.capa.ident mehod

Hi R-experts, I apologize for the multiple posting... My previous post seems to have been eaten by the server because of embedded code. First of all, I am working on Windows XP, using R version 2.6.2. I got this error message " Error in drop(.Call("La_dgesv", a, as.matrix(b)

Re: [R] post hoc statistical power

Many thanks to John Fox and Charles C. Berry. Both point to the paper by Hoenig and Heisey (email from John Fox below): > Another paper critical of post-hoc power calculations is Hoenig and Heisey, > 2001. "The Abuse of Power: The Pervasive Fallacy of Power Calculations for > Data Analysis." The

Re: [R] "preferred" version of Linux for R?

On 8 April 2008 at 19:06, Mitchell Maltenfort wrote: | If your criterion is use of resources, then you might be more | interested in which flavor of Linux makes it easy to turn the X-window | server off. No windows or mouse, more RAM and CPU cycles for R. As I | recall, the only Linux where that

Re: [R] Overlapping distributions (populations) - assigning an individual to a population?

> Rolf, > > > On Wed, 2008-04-09 at 10:57 +1200, Rolf Turner wrote: >> On 9/04/2008, at 10:30 AM, Phil Rhoades wrote: >> >> > People, >> > >> > Say a particular measure of an attribute for individuals in different >> > populations gives a set of overlapping normal distributions (one >> > distributi

Re: [R] post hoc statistical power

Dear Yuelin, Another paper critical of post-hoc power calculations is Hoenig and Heisey, 2001. "The Abuse of Power: The Pervasive Fallacy of Power Calculations for Data Analysis." The American Statistician, 55, 19-24. Hoenig and Heisey show that so-called "observed power" is simply a re-expression

Re: [R] Overlapping distributions (populations) - assigning an individual to a population?

Rolf, On Wed, 2008-04-09 at 10:57 +1200, Rolf Turner wrote: > On 9/04/2008, at 10:30 AM, Phil Rhoades wrote: > > > People, > > > > Say a particular measure of an attribute for individuals in different > > populations gives a set of overlapping normal distributions (one > > distribution per popul

[R] chi-square test

Hi R-users, I would like to find the goodness of fit using Chi-suare test for my data below: xobs=observed data, xtwe=predicted data using tweedie, xgam=predicted data using gamma > xobs <- c(223,46,12,5,7,17) > xtwe <- c(217.33,39,14,18.33,6.67,14.67) > xgam <- c(224.67,37.33,12.33,15.33,5.33,15)

Re: [R] "preferred" version of Linux for R?

There are binaries on CRAN for several distros. There is the R-sig-debian list for debian users. As far as I know, there are no mailing lists for other distros. This may indicate a preference for debian among R users. Or it may mean that debian is so difficult to use with R that it requires a maili

Re: [R] permutation test assumption?

On Wed, 9 Apr 2008, João Fadista wrote: > But in that paper they refer to another > (http://www.bepress.com/cgi/viewcontent.cgi?article=1053&context=ucbbiostat) > where they say that "permutation distribution produces an asymptotically > correct null distribution if (i) the sample sizes are equal.

Re: [R] Split

Assuming that this data is character, then this will work: (it removes the last period in the string) > x <- scan('clipboard', what='') Read 4 items > x [1] "1234." "2345." "3435." "4343." > sub("\\.$", "", x) [1] "1234" "2345" "3435" "4343" > On Tue, Apr 8, 2008 at 11:30 AM, Kris Ghosh <[EMAIL

Re: [R] Combining many csv files into one and adding a column with an id of each csv file read

I blew it. Forgot to use the parameter. Here is the revised copy: filenames <- Sys.glob("*.csv") # however you get the list of file allData <- lapply(filenames, function(.file){ dat<-read.csv(.file, header=F) dat<-dat[c(-1:-3),c(-1,-4,-5,-6,-7,-9,-10,-11,-12)] names(dat)<-c('time','T1

Re: [R] Combining many csv files into one and adding a column with an id of each csv file read

Here is one way of doing it. Read the data into a list and then use 'do.call(rbind...': filenames <- Sys.glob("*.csv") # however you get the list of file allData <- lapply(filenames, function(.file){ dat<-read.csv(filename, header=F) dat<-dat[c(-1:-3),c(-1,-4,-5,-6,-7,-9,-10,-11,-12)]

Re: [R] Separate y-limits in xYplot panels II

Hi John, I think even with scales=list(y=list(relation="free") you have to specify the y-limits for each panel, i.e. add something like ylim=list(c(0,1200),c(0,1)) to your code. Best regards, Christoph Wednesday, April 9, 2008, 12:37:43 AM, you wrote: > Hello, > I have been trying to use

Re: [R] "preferred" version of Linux for R?

If your criterion is use of resources, then you might be more interested in which flavor of Linux makes it easy to turn the X-window server off. No windows or mouse, more RAM and CPU cycles for R. As I recall, the only Linux where that might be a problem is Ubuntu, because there's no convenient w

Re: [R] findGlobals on apply

On Wed, 9 Apr 2008, Christophe Genolini wrote: Hi Luke Thanks for all these explanation, things are clearer. Let me go back on my initial problem, that was, as a programmer, I would like to have a tool to detect typo by detecting globals variables: I get that findGlobals is not design for th

Re: [R] Overlapping distributions (populations) - assigning an individual to a population?

On 9/04/2008, at 10:30 AM, Phil Rhoades wrote: > People, > > Say a particular measure of an attribute for individuals in different > populations gives a set of overlapping normal distributions (one > distribution per population). If I then measure this attribute in > a new > individual - how d

Re: [R] permutation test assumption?

> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of João Fadista > Sent: Tuesday, April 08, 2008 3:10 PM > To: ONKELINX, Thierry; r-help@r-project.org > Subject: Re: [R] permutation test assumption? > > Dear Thierry, > > Thanks for the reply. But as y

Re: [R] permutation test assumption?

But in that paper they refer to another (http://www.bepress.com/cgi/viewcontent.cgi?article=1053&context=ucbbiostat) where they say that "permutation distribution produces an asymptotically correct null distribution if (i) the sample sizes are equal...". Best regards, João ___

[R] Separate y-limits in xYplot panels II

Hello, I have been trying to use xYplot in Hmisc to graph plots, allowing each panel to have a y-axis dependent on the data in the panel. Following the advice from the R-Help list, message: [R] separate y-limits in xYplot panels on Wed, 30 May 2007 08:12:06 -0700 (PDT), I used scales=list(y=list(

[R] Overlapping distributions (populations) - assigning an individual to a population?

People, Say a particular measure of an attribute for individuals in different populations gives a set of overlapping normal distributions (one distribution per population). If I then measure this attribute in a new individual - how do I assess the likelihood of this new individual belonging to ea

Re: [R] permutation test assumption?

On Wed, 9 Apr 2008, João Fadista wrote: > Dear Thierry, > > Thanks for the reply. But as you may read in the paper > http://bioinformatics.oxfordjournals.org/cgi/content/abstract/22/18/2244 > when the sample sizes are not the same there may be an increase in the > Type I error rate. No, this is f

Re: [R] permutation test assumption?

Dear Thierry, Thanks for the reply. But as you may read in the paper http://bioinformatics.oxfordjournals.org/cgi/content/abstract/22/18/2244 when the sample sizes are not the same there may be an increase in the Type I error rate. Comments will be appreciated. Best regards, João Fadista

Re: [R] findGlobals on apply

Hi Luke Thanks for all these explanation, things are clearer. Let me go back on my initial problem, that was, as a programmer, I would like to have a tool to detect typo by detecting globals variables: I get that findGlobals is not design for that. I did not realy understand the use of checkUsa

Re: [R] Problem with NA data when computing standard error

On Tue, Apr 8, 2008 at 12:44 PM, LeCzar <[EMAIL PROTECTED]> wrote: > > Hey, > > I want to compute means and standard errors as two tables like this: > > se<-function(x)sqrt(var(x)/length(x)) > > The missings are not your main problem. The command var computes the variance-covariance matrix.

Re: [R] save multiple pages lattice-plot

On 4/8/08, Katell HAMON <[EMAIL PROTECTED]> wrote: > Dear R-helper, > > I am trying to save some plots I made with lattice, but I have several > pages printed for each call. > > My layout is something like c(col=2,row=4,pages=12) and I didn't find a > way to save those plots, because I am usual

[R] "preferred" version of Linux for R?

Is there a recommended/preferred version of Linux for using with R? Is there one version of Linux that R-users prefer, and/or that works "better" with R? I am working with "large" datasets, and hope to take advantage of as much RAM as reasonable (8-32gb?). Thanks in advance!

Re: [R] reading in a Fortran binary array

See ?readBin The format of a Fortran unformatted file is compiler- and OS-dependent, but if you know what it is, readBin() can read it. On Tue, 8 Apr 2008, Alexy Khrabrov wrote: > Greetings -- I'd like to avoid converting a Fortran array of floats > into ASCII and back reading it in R. Further

Re: [R] Problem with NA data when computing standard error

Hi there, Perhaps se<-function(x)sqrt(var(x,na.rm=T)/sum(!is.na(x))) object1<-as.data.frame.table(tapply(Data[Year=="1999"],list(Group[Year=="1999"],Season[Year=="1999"]),mean)) object2<-as.data.frame.table(tapply(Data[Year=="1999"],list(Group[Year=="1999"],Season[Year=="1999"]),se)) Hope this

Re: [R] how to check if a variable is preferentially present in a sample

Hi Tania, I think it could be. I tried a solution based on your data set using a chi-squared approach. Here is what I got: # # Data set set.seed(123) d <- data.frame(cbind(val=rnorm(1:10)^2, group=sample(LETTERS[1:5],100,repl=TRUE))) d[,"val"]<-as.numeric(as.character(d$val)) #

Re: [R] fit points to plane

Bert Gunter wrote: But probably not a good idea, as there are likely to be some scientfically "interesting" local inhomogeneities that an appropriate nonparametric smoother (e.g. splines,...) could reveal. ok... im working on laser scanner cloud of points... the model im working with is som

Re: [R] fit points to plane

Rolf Turner wrote: On 9/04/2008, at 8:28 AM, Luca Penasa wrote: Is there a function for obtaining the best-fitting plane from a large number of points (something like 25.000 points)? ?lm what does it mean? ## Attention

Re: [R] diagonally fill a rectangle with color gradient

Hi tom soyer wrote: > Hi, > > Is it possible to diagonally fill a rectangle with a color gradient? I > noticed that the gradient.rect of plotrix could fill a rect either up and > down or from side to side. I am looking for something similar but fills > diagonally instead, e.g., from the upper le

Re: [R] distance matrix as text file - how to import?

On Tue, Apr 8, 2008 at 1:50 PM, Hans-Jörg Bibiko <[EMAIL PROTECTED]> wrote: > Dear all, > > I have -hopefully- a tiny problem. > > I was sent a text file containing a distance matrix à la: > > 1 > 2 3 > 4 5 6 > Try this! I put your test data in "text.txt" and voila: mat <- matrix(0, 3,3)

[R] Smoothing 3 D data

Hey. I have a set of data points (x1,y1,z1; x2,y2,z2;...xn,yn,zn). I need to smooth these in 3D. For example if these were in 2 D then one would use inverse distance weighting or moving averages. Does anyone know of any funtion in R that can be used to do this (Using 3D dat

Re: [R] fit points to plane

But probably not a good idea, as there are likely to be some scientfically "interesting" local inhomogeneities that an appropriate nonparametric smoother (e.g. splines,...) could reveal. -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rolf

[R] Problem with NA data when computing standard error

Hey, I want to compute means and standard errors as two tables like this: se<-function(x)sqrt(var(x)/length(x)) object1<-as.data.frame.table(tapply(Data[Year=="1999"],na.rm=T,list(Group[Year=="1999"],Season[Year=="1999"]),mean)) object2<-as.data.frame.table(tapply(Data[Year=="1999"],na.rm=

[R] Separate y-limits in xYplot panels II

Hello, I have been trying to use xYplot in Hmisc to graph plots, allowing each panel to have a y-axis dependent on the data in the panel. Following the advice from the R-Help list, message: [R] separate y-limits in xYplot panels on Wed, 30 May 2007 08:12:06 -0700 (PDT), I used scales=list(y=li

[R] save multiple pages lattice-plot

Dear R-helper, I am trying to save some plots I made with lattice, but I have several pages printed for each call. My layout is something like c(col=2,row=4,pages=12) and I didn't find a way to save those plots, because I am usually using savePlot() which from what I saw can only save one page

[R] Split

I have data in the form of a column such as 1234. 2345. 3435. 4343. I want to have the data in this form ..i.e to remove the "dot" at the end of each number above. 1234 2345 3435 4343 I am trying to use split but it is not working. Any suggestions? [[alternative HTML version deleted]] __

Re: [R] fit points to plane

On 9/04/2008, at 8:28 AM, Luca Penasa wrote: > Is there a function for obtaining the best-fitting plane from a > large number of points (something like 25.000 points)? ?lm ## Attention:\ This e-mail message is privileged an

Re: [R] simple graphing question

On Tue, Apr 8, 2008 at 2:18 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > #copy and paste this into R > f <- (structure(list(TKN = c(0.103011025, 0.018633208, 0.104235702, > 0.074537363, 0.138286096), RM = c(215, 198, 148, 119, 61)), .Names = c("TKN", > "RM"), class = "data.frame", row.names =

Re: [R] simple graphing question

This is maybe not the most elegant way, but it does de job. You first put f$RM values in negative form. Then you plot your graph without the x axis labels. After, you create the labels you want. Try this : f <- (structure(list(TKN = c(0.103011025, 0.018633208, 0.104235702, 0.074537363, 0.1382

Re: [R] simple graphing question

On Tue, Apr 8, 2008 at 3:52 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > Sorry, to have the x-axis to go from 200 to 0 or to reverse the x > points- the line starts on the left hand side of the graph at x=215, > y=0.10301103 ... and end with x=61, y=0.13828610. does this make > sense?

Re: [R] simple graphing question

Try: plot(f$RM~f$TKN, type="b") On Tue, Apr 8, 2008 at 4:18 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > #copy and paste this into R > f <- (structure(list(TKN = c(0.103011025, 0.018633208, 0.104235702, > 0.074537363, 0.138286096), RM = c(215, 198, 148, 119, 61)), .Names = > c("TKN", > "RM")

Re: [R] plot function / par() settings

On Tue, Apr 8, 2008 at 12:47 PM, Pedro Mardones <[EMAIL PROTECTED]> wrote: > Dear all; > > I'm trying to create a 2 x 3 plot (something I know like lattice can > do better) using the plot function. However; I'm not sure how to make > the width of the plots to be the same on each column. I guess

Re: [R] Metropolis acceptance rates

Hi Michael, Actually I was wrong in the initial mail. rbprobitGibbs uses the Gibbs Sampler. I think the rhierbinLogit does the Logit model which uses a Metropolis hastings approach uses the acceptance rate and the function outputs the rejection rate in the output list. So you are not required t

[R] simple graphing question

#copy and paste this into R f <- (structure(list(TKN = c(0.103011025, 0.018633208, 0.104235702, 0.074537363, 0.138286096), RM = c(215, 198, 148, 119, 61)), .Names = c("TKN", "RM"), class = "data.frame", row.names = 25:29)) plot(f$TKN~f$RM, type="b") I would like to reverse the X-Axis. How do I do

[R] reading in a Fortran binary array

Greetings -- I'd like to avoid converting a Fortran array of floats into ASCII and back reading it in R. Furthermore it's much faster to dump large arrays in binary, as they take up much less space with full precision -- many decimal points take up many bytes in ASCII versus four or eight

[R] Combining many csv files into one and adding a column with an id of each csv file read

Dear R experts, I have been looking into the help-pages and old questions from the R-Help site, but the options offered there don't seem to work in my case. First of all, I am working on Windows XP, using R version 2.6.2. I am attaching two csv files as an example of how the data I am tra

[R] distance matrix as text file - how to import?

Dear all, I have -hopefully- a tiny problem. I was sent a text file containing a distance matrix à la: 1 2 3 4 5 6 Now I wanted to import these data into a dist object to, let's say, do 'plot(hclust(v))'. My first naïve approach was to scan the text file in order to get a vector v. Then I

[R] diagonally fill a rectangle with color gradient

Hi, Is it possible to diagonally fill a rectangle with a color gradient? I noticed that the gradient.rect of plotrix could fill a rect either up and down or from side to side. I am looking for something similar but fills diagonally instead, e.g., from the upper left corner to the bottom right. Doe

[R] fit points to plane

Is there a function for obtaining the best-fitting plane from a large number of points (something like 25.000 points)? Thank you... Luca -- Email.it, the professional e-mail, gratis per te: http://www.email.it/f Sponsor: Gioca con i Supereroi Marvel sul cellulare! Clicca qui: http://adv.email

Re: [R] plotmath "overstrikes" in output on a Linux system

On Tue, Apr 8, 2008 at 4:29 AM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > Looks likes the the laptop is using different fonts with incorrect font > metrics. This could happen because it has a different screen resolution, or > one of the systems is set to use scalable fonts or it is giving met

[R] plot function / par() settings

Dear all; I'm trying to create a 2 x 3 plot (something I know like lattice can do better) using the plot function. However; I'm not sure how to make the width of the plots to be the same on each column. I guess the answer maybe obvious but I haven't been able to figure it out. I'll appreciate any

Re: [R] lme and confidence intervals

Cristian Carranza hotmail.com> writes: > After fitting a mixed effects model to repeated measurements data set, and after several unsuccessful > atempts to make a simple plot of the confidence interval for the fitted model, I gave up and now I am asking > for help in this useful list. > > Could a

Re: [R] Metropolis acceptance rates

Hi Michael, If you are using the standard logit or probit model, it is fairly easy to save the acceptance rate after each draw. I would recommend using the "bayesm" package as the source code is easy to manipulate. For instance in the probit function (rbprobitGibbs), you need to include a co

Re: [R] findGlobals on apply

Yes, based on that it seems understandable although initially suprising. On Tue, Apr 8, 2008 at 12:56 PM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > On 4/8/2008 7:42 AM, Gabor Grothendieck wrote: > > I think what he is referring to is that findGlobals lists mean under > > variables rather than fu

Re: [R] Help me about singular error nls

jarod_v6 libero.it> writes: > > I have that problem: > this is file : PBS.txt > Time RFU > 0.27.3021 > 0.080026.1565 > rm(list=ls()) > print("ls()") > #carica Dati con file di testo > > b <-read.table("PBS.txt", sep="\t", header= TRUE) > print("b") > b > nlmod1 <- nls(R

Re: [R] post hoc statistical power

On Tue, 8 Apr 2008, Yuelin Li wrote: > I remember reading about post hoc statistical power on R-help. But I > can't seem to find them with RSiteSearch("post hoc statistical power") > and variations of it. > > I would like to learn more about post hoc statistical power, its > meaningfulness, advan

[R] nonlinear regression model to data with power exponential distributed errors

Hi, are you aware of any R packages that can fit a nonlinear regression model to a dataset with power exponential distributed random errors? Daniel __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posti

Re: [R] findGlobals on apply

On 4/8/2008 7:42 AM, Gabor Grothendieck wrote: > I think what he is referring to is that findGlobals lists mean under > variables rather than functions when you do this with his f: > >> findGlobals(f, FALSE) > $functions > [1] "{" "apply" > > $variables > [1] "mean" Yes, I understood that, a

Re: [R] Practical Data Limitations with R

Millions of rows can be a problem if all is loaded into memory, depending on type of data. Numeric should be fine but if you have strings and you would want to process based on that column (string comparisons etc) then it would be slow. You may want to combine sources outside - stored procedures

[R] post hoc statistical power

I remember reading about post hoc statistical power on R-help. But I can't seem to find them with RSiteSearch("post hoc statistical power") and variations of it. I would like to learn more about post hoc statistical power, its meaningfulness, advantages and disadvantages. I thought the issue was

[R] RSPerl & threads II

My previous post seems to have been eaten by the server because of embedded code. You can now find that code at: http://tinyurl.com/3o88j2 I aoplogize for the double-posting... -Jay Hi all, I'm trying to use the RSPerl module in a script that uses threads. I am able to call R functions withou

[R] Goodness of fit tests

Dear all, Is it possible to use chisq.test(x,p) (where x = c(N1,...,Nk) are the counts in the k "bins" and p is a vector of the probabilities that a random observation will fall in the various bins) to test the goodness of fit of a certain distribution for which we had to estimate the paramete

[R] RSPerl & threads

Hi all, I'm trying to use the RSPerl module in a script that uses threads. I am able to call R functions without a problem when I don't use threads. However, using threads, I get varying errors depending on how I call the R functions. If I call the R::initR function in the "boss" thread and then

Re: [R] re garding Garch prediction mechanism

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. In particular, I suggest you use a published data set from a package like tseries, fEcofin, or FinTS. If you want to compare results with Matlab,

Re: [R] Change the position of panel strips in a lattice plot.

On 4/8/08, Gustaf Rydevik <[EMAIL PROTECTED]> wrote: > Hi all, > > > In lattice plots, is there any option to position the panel strips > with text below each subgraph, instead of above? No. -Deepayan > i.e. in: > Depth <- equal.count(quakes$depth, number=8, overlap=.1) > xyplot(lat ~ long

[R] Help me about singular error nls

Dear List I have that problem: this is file : PBS.txt TimeRFU 0. 27.3021 0.0800 26.1565 0.1600 25.4836 0.2400 25.2104 0.3200 24.7967 0.4000 24.6310 0.4800 24.3472 0.5600 24.2905 0.6400 24.1394 0.7200 24.0244 0.8000 23.8532 0.8800 23.7908 0.9600 23.6668 1.0400 23.7458 1.1200

[R] Help me about singular error nls

Dear List I have that problem: this is file : PBS.txt TimeRFU 0. 27.3021 0.0800 26.1565 0.1600 25.4836 0.2400 25.2104 0.3200 24.7967 0.4000 24.6310 0.4800 24.3472 0.5600 24.2905 0.6400 24.1394 0.7200 24.0244 0.8000 23.8532 0.8800 23.7908 0.9600 23.6668 1.0400 23.7458 1.1200

Re: [R] Practical Data Limitations with R

> We are new to R and evaluating if we can use it for a project we need to > do. We have read that R is not well suited to handle very large data > sets. Assuming we have the data prepped and stored in an RDBMS (Oracle, > Teradata, SQL Server), what can R reasonably handle from a volume > per

[R] how to check if a variable is preferentially present in a sample

Dear All, I do apologise if this question is out of place for this list but I've tried searching mailing lists and read "Introductory Statistics with R" by Peter Dalgaard, but couldn't find any hints on solving my question below: I have a data frame (d) of values which I will rank in decrea

Re: [R] Practical Data Limitations with R

On Tue, Apr 08, 2008 at 09:26:22AM -0500, Jeff Royce wrote: > We are new to R and evaluating if we can use it for a project we need to > do. We have read that R is not well suited to handle very large data > sets. Assuming we have the data prepped and stored in an RDBMS (Oracle, > Teradata, SQL S

Re: [R] Practical Data Limitations with R

Dear Jeff, R works fine for 22 rows that i tested on a home PC with XP . Memory is limited to hardware that you have. I suggest beefing up RAM to 2 GB and hard disk space and then working it out. I evaluated R too on my site www.decisionstats.com and I found it comparable if not better to SPS

Re: [R] Question Regarding 'pipe'

On Tue, Apr 08, 2008 at 06:52:33AM -0500, [EMAIL PROTECTED] wrote: > Reason, I need to do this in awk and not R: > > Let's say 'x' is the tabular representation of a sparse contingency table > > x >x10 x11 x12 x13 x14 x15 x16 x17 x18 x19 Freq > 10 0 0 2 2 2 0 2 0 2 54 >

[R] Practical Data Limitations with R

We are new to R and evaluating if we can use it for a project we need to do. We have read that R is not well suited to handle very large data sets. Assuming we have the data prepped and stored in an RDBMS (Oracle, Teradata, SQL Server), what can R reasonably handle from a volume perspective? Ar

Re: [R] How to access the attributes of a ggplot?

> Thanks for replying. I know that I can directly modify some of the > properties of the plot object, but I was more interested in querying > the current plot properties, something like: > > xlimits <- getggopts(pobj, "x_scale_limits") > > Is there anything like this implemented? Ah, that's m

[R] Porting from 1 stats platform to another

Dear All, I was recommended this book " http://oit.utk.edu/scc/RforSAS&SPSSusers.pdfand found it very good as a 3-in-1 reference guide , and also as migration guide from 1 platform to another. Many Thanks to the people who helped me in this. It took me some time to find out the R GUI though, tha

Re: [R] How to pack my stuff into a package (library, collection)?

> "TL" == Tribo Laboy <[EMAIL PROTECTED]> > on Tue, 8 Apr 2008 22:15:13 +0900 writes: TL> okey-dokey, one more problem resolved. TL> Keeping one documentation .Rd file for each R source file. not at all ... you still think "Matlab" One R source file typically contains severa

[R] Change the position of panel strips in a lattice plot.

Hi all, In lattice plots, is there any option to position the panel strips with text below each subgraph, instead of above? i.e. in: Depth <- equal.count(quakes$depth, number=8, overlap=.1) xyplot(lat ~ long | Depth, data = quakes) ,is there any way to make "Depth" appear below the subgraphs, ins

Re: [R] Mode Vs Class

'mode' is a mutually exclusive classification of objects according to their basic structure. The 'atomic' modes are numeric, complex, charcter and logical. Recursive objects have modes such as 'list' or 'function' or a few others. An object has one and only one mode. 'class' is a property assig

[R] basehaz and newdata

The 'basehaz' function is just a wrapper for survfit, and includes only some of the arguments of the former. It's main reason for existence is that another more well known (but inferior :-) package uses that term. I don't understand the final comment in the exchange, however: > I had been hopi

Re: [R] permutation test assumption?

Dear João, You can do permutation tests on an unbalanced design. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwa

[R] permutation test assumption?

Dear all, Can I do a permutation test if the number of individuals in one group is much bigger than in the other group? I searched the literature but I didin´t find any assumption that refers to this subject for permutation tests. Best regards João Fadista Ph.d. student

Re: [R] How to pack my stuff into a package (library, collection)?

okey-dokey, one more problem resolved. Keeping one documentation .Rd file for each R source file. Thanks! TL __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/post

Re: [R] How to pack my stuff into a package (library, collection)?

On 08/04/2008 7:08 AM, Tribo Laboy wrote: > Thanks all for the help and suggestions. I am little by little finding > my way. I have another question to the people who use the R packaging > system. Say I have a function called "myfun.R". I guess you mean you have a source file myfun.R, containing

[R] Metropolis acceptance rates

Is there a way to recover Metropolis-step acceptance rates AFTER completing posterior draws? The immediate application is in the probit.bayes and logit.bayes models used by Zelig... which I believe is merely calling MCMCpack. So one strategy, to which I am fixing to resort, is to call, say,

Re: [R] Install spatstat package on gentoo

Luca Penasa wrote: > Anybody know how to install R-spatstat package on a gentoo based system? Same as on any other system, see the manual R Installation and Administration and please do read the posting guide. Best, Uwe Ligges > Thanks you > Luca > > > > -- > Email.it, the professional

[R] Mode Vs Class

Hi R, Just came across the 'mode' of an object. What is the basic difference between ?class and ?mode...For example: d=data.frame(a=c(1,2),b=c(5,6)) class(d) [1] "data.frame" mode(d) [1] "list" But, c=c(2,3,5,6,7) class(c) [1] "numeric" mode(c) [1] "numeric" Could anyone

Re: [R] sample size for log-rank test with more than 2 groups

Essioux, Laurent wrote: > Hi everyone, > > I was wondering whether extension of the current spower function for > Hmisc were existing? > > My current focus is to calculate sample size based on the log-rank test > with more than 2 groups (with/without trend) > > Taking into account the loss of fo

Re: [R] Question Regarding 'pipe'

Reason, I need to do this in awk and not R: Let's say 'x' is the tabular representation of a sparse contingency table > x x10 x11 x12 x13 x14 x15 x16 x17 x18 x19 Freq 10 0 0 2 2 2 0 2 0 2 54 51 0 0 2 2 2 0 2 1 21 61 0 1 0 2 2 0

Re: [R] findGlobals on apply

On Tue, 8 Apr 2008, Christophe Genolini wrote: f <- function(x){apply(x,2,mean)} findGlobals(f) mean is a global variable, so findGlobals gets it right. That sound strange to me: a "variable" is something that vary... mean does not vary. maen will ge an argument that is a line of x and will

Re: [R] findGlobals on apply

I think what he is referring to is that findGlobals lists mean under variables rather than functions when you do this with his f: > findGlobals(f, FALSE) $functions [1] "{" "apply" $variables [1] "mean" On Tue, Apr 8, 2008 at 7:14 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > Christophe Ge

Re: [R] How to pack my stuff into a package (library, collection)?

Yesterday I wrote: > I took a look at this today. You get an error message but the package > is still installed without the CHM compiler, so that's less urgent. > > I did add a menu entry to install a source package from a directory: > > Packages | Install source package from local folder... > >

[R] Install spatstat package on gentoo

Anybody know how to install R-spatstat package on a gentoo based system? Thanks you Luca -- Email.it, the professional e-mail, gratis per te: http://www.email.it/f Sponsor: Una BMW Z4 Roadster a 10€? Prova gratis su Bidplaza.it! * Clicca qui: http://adv.email.it/cgi-bin/foclick.cgi?mid=7673&d

  1   2   >