[R] biglm.big.matrix: Problem with weighting

Hello colleagues, I have tried to use the package bigmemory, biganalytics and biglm. I want to specify a multivariate regression with a weight. I have imported a large dataset with the library(bigmemory). I load the library (biglm) and specified a regression with a weight. But I get

Re: [R] Package testing

Duncan, thanks, if the header is skipped, then it is much better. It seems that I will be able to run the tests without actually running 'R CMD check' and extract the results, so it all looks good now. Best Regards, Gabor On Tue, Jun 15, 2010 at 7:25 PM, Duncan Murdoch murdoch.dun...@gmail.com

[R] fgev + negative location parameter for positive data?

Hi R, I have a vector with all positive values and have fit the GEV distribution (shape parameter=0) to the data and I get the negative value of the location parameter. Is this possible? library(evd) r=read.table(clipboard,header=F)[,1] length(r) [1] 2087 sum(r=0) #All data points

Re: [R] Error when callin g C-Code

This is not the appropriate list for C-level questions. Please study the posting guide and repost on the right list. Also, you might want to also show how you actually call this from R. Are you by any chance using .External instead of .Call ? Romain Le 15/06/10 23:20, Fabian Zäpernick a

[R] disabling scientific number format

Hello R users, I am a bit puzzled by the fact that options(digits=22) x0 - 69880 y0 - 26185 x1 - 69600 y1 - 22937 number - paste(x0,y0,x1,y1,sep=) number - as.numeric(travel) print(number) gives number in scientific notation (some precision is wasted by having a e+19 at the tail):

[R] Exporting multiple plots

Sir, I want to export 10, say graphs using 'for' loop. What will be the command and format for exporting multiple plots. Regards, Suman Dhara [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Fitting Gamma distributionm

I'm looking for goodness of fit tests for gamma distributions with large data sizes. I have a matrix with around 10,000 data values in it and i have fitted a gamma distribution over a histogram of the data. The problem is testing how well that distribution fits. Chi-squared seems to be used

Re: [R] Graphics question: How to create a changing smudge factor for overlapping lines?

Hi Tal, I you use ggplot you can use the alpha command to make lines transparent. The nice thing is that when they overlap, the transparency adds up. I use this a lot to visualize outcomes from ensemble modelling (e.g. time series of RMSE). A small example: library(ggplot2) dat =

Re: [R] an alternative to R for nonlinear stat models

On 06/16/2010 07:35 AM, benedikt.g...@ieu.uzh.ch wrote: Hi I implemented the age-structure model in Gove et al (2002) in R, which is a nonlinear statistical model. However running the model in R was very slow. So Dave Fournier suggested to use the AD Model Builder Software

[R] Odp: Matrix to Vector

Hi r-help-boun...@r-project.org napsal dne 15.06.2010 19:51:39: Hi All, I am trying to turn a Matrix into a vector for analysis purposes. I need to select only certain columns from the entire matrix for the vector (intraday time intervals). Also I need to transpose the Matrix (so times

[R] data frame

Dear list, I have the following problem. I have a data frame like this CLUSTERYEAR variableDelta R_pivot M1 2005 EC01 NA NA M1 2006 EC012

[R] Backslash in paste() function

Hi, I'm trying to build a vector of latex commands. However, I need the command strings to begin with a backslash \. I have: test - c('foo','bar') and I need to rebuild the array, encapsulating the text items with latex stuff, like this: paste(\parbox[b]{3cm}{, test, }, fill=TRUE)

[R] Fitting Gamma distribution

I'm looking for goodness of fit tests for gamma distributions with large data sizes and for different data. I have a matrix with around 4.000 data values in it and i have fitted a gamma distribution with fitdistr. You can see the example: fitdistr(corpo,gamma,lower=0.001) Errore in

Re: [R] multiple plots without for loops

Hi I use an option with PDF, which can have multi page output. You can pdf(some_graphs.pdf) # do all your plotting let say in cycle something like for( i in 1:n) plot(something[,1], something[,i+1], ...) dev.off() You will get one multipage pdf document in working directory which you can

Re: [R] meta analysis with repeated measure-designs?

Dear Gerrit, If the correlations of the dependent effect sizes are unknown, one approach is to conduct the meta-analysis by assuming that the effect sizes are independent. A robust standard error is then calculated to adjust for the dependence. You may refer to Hedges et. al., (2010) for more

[R] prompt question

Hi,   R has a number of functions to aid documentation: prompt (for functions) promptMethods promptClass package.skeleton   My question: does package.skeleton generate the basis for a documentation for functions, methods and (S4) classes? Or only for (generic) functions? I cannot this this here

[R] Question

Dear Mr. for writing program about Gibbs sampling, i have a question. if i want to generate data from Exponential distribution but range of X is restricted, how can i do? regards, A.Rashidi   [[alternative HTML version deleted]] __

[R] Antwort: Re: an alternative to R for nonlinear stat models

Hi Paul Oh ok sorry, I send you below a copy of the R code I was using. It's very possible that my programming is slow as well. I'm very happy to learn. Sorry the code is a bit long, but its structured like this and cotains three main functions and the optimization: 1.

Re: [R] Graphics question: How to create a changing smudge factor for overlapping lines?

Hi Paul, Thank you for the ggplot2 suggestions. Please note that during the last few hours, Hadley has put up his own ggplot2 implementation of the clustergram. I updated about it in the original post:

[R] Decoding raw vectors every other byte

I have to decode raw vectors read from a connection like this: +++ fl-file.choose() fb-file.info(fl)$size fc-file(fl) open(fc,open='rb') seek(fc,where=offset) dat-readBin(fc,'raw',n=(fb-offset),size=1,signed=FALSE) close(fc) ++ The decoding performs a combination

Re: [R] Decoding raw vectors every other byte

x - 1:10; dim(x) - c(length(x)/2, 2); print(x); [,1] [,2] [1,]16 [2,]27 [3,]38 [4,]49 [5,]5 10 str(x[,1]); str(x[,2]); If length(x) is odd, some care is needed, e.g. x - matrix(x, ncol=2, byrow=FALSE); /H On Wed, Jun 16, 2010 at 1:15 PM, Alex van der

Re: [R] cross validation of SVM

From ?svm: cross if a integer value k0 is specified, a k-fold cross validation on the training data is performed to assess the quality of the model: the accuracy rate for classification and the Mean Squared Error for regression Uwe Ligges On 15.06.2010 23:14, Amy Hessen wrote: hi,

Re: [R] regression analysis p value

On 15.06.2010 23:17, Liang Cheng wrote: Hello:) I have 4 points in a 2-dimension coordinate. how can I do the regression analysis and calculate the p-value£¿ Well, you can do by uising lm(), but you should not. Reporting the values of your 4 points makes more sense than calculating

Re: [R] Backslash in paste() function

Just double all the backslashes and you are fine. In order to see the outcome, use cat() (not print). Uwe Ligges On 16.06.2010 09:49, Stefan Petersson wrote: Hi, I'm trying to build a vector of latex commands. However, I need the command strings to begin with a backslash \. I have:

Re: [R] disabling scientific number format

On 16.06.2010 08:59, David Croll wrote: Hello R users, I am a bit puzzled by the fact that options(digits=22) x0- 69880 y0- 26185 x1- 69600 y1- 22937 number- paste(x0,y0,x1,y1,sep=) number- as.numeric(travel) print(number) gives number in scientific notation (some precision is wasted

Re: [R] Exporting multiple plots

On 16.06.2010 08:03, suman dhara wrote: Sir, I want to export 10, say graphs using 'for' loop. What will be the command and format for exporting multiple plots. Depends on the graphics device you are going to use, e.g. if you want all plots in one pdf or as separate bitmaps or ... See the

Re: [R] Backslash in paste() function

2010/6/16 Uwe Ligges lig...@statistik.tu-dortmund.de: Just double all the backslashes and you are fine. In order to see the outcome, use cat() (not print). Is R ever going to get an alternate way to specify literal character strings? This has been a real annoyance for years.

[R] prcomp() and the lenght of PC:s

Hi, I would like to know whether there is some deeper rationale behind or is it just an established practice that the lenghts of principal components, giving for example by prcomp-function, are normalised to 1? Best regards, Atte Tenkanen University of Turku, Finland Department of Musicology

[R] Binning function in R

Is there some function in R that does what Interactive Grouping node in SAS Enterprise Miner does? It makes variable binning using WOE(weight of evidence) for scoring modeling. I’ve found 3 binning function in R : rattle::binning

[R] memory management in R

I have volunteered to give a short talk on memory management in R to my local R user group, mainly to motivate myself to learn about it. The focus will be on what a typical R coder might want to know ( e.g. how objects are created, call by value, basics of garbage collection ) but I want

Re: [R] working with zoo time index ??

Hi thanks Let say data are written like this: 1990-01-01 10:01:00 , 0.910 1990-01-01 10:03:00 , 0.905 Would it be ok to read it with theses lines or is better to use your way? tmp - read.table(demo2.txt, sep = ,) z - zoo(tmp[, 2], as.Date(as.chron(tmp[, 1]), format = %Y-%m-%d %H:%M:%S))

[R] generating samples from multivariate distributions

Sir, I want to draw random from any multivariate disrtibution. Is there any function in R to do this? Regards, Suman Dhara [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] working with zoo time index ??

I can't see where you check the date. -- View this message in context: http://r.789695.n4.nabble.com/working-with-zoo-time-index-tp2255804p2257226.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] Problem with the recode function

Dear John, Thanks a lot for the time you spent on my problem. I don't believe you can do something to avoid this kind of problem. I don't know if it is technically possible but I wonder if when we load Rcmdr plug-ins from the Rcmdr menu, it wouldn't be possible during the restart of Rmcdr

Re: [R] working with zoo time index ??

On Wed, Jun 16, 2010 at 7:58 AM, skan juanp...@gmail.com wrote: Hi thanks Let say data are written like this: 1990-01-01 10:01:00 ,  0.910 1990-01-01 10:03:00 ,  0.905 Would it be ok to read it with theses lines or is better to use your way? tmp - read.table(demo2.txt, sep = ,) z -

Re: [R] memory management in R

You might want to mention/talk about packages that enhance R's ability to work with less RAM / more data, such as package SOAR (transparently moving objects between RAM and disk) and ff (which allows vectors and dataframes larger than RAM and which supports dense datatypes like true boolean,

Re: [R] Parallel computing on Windows (foreach)

On 15-Jun-10 17:07, Sergey Goriatchev wrote: Hello, I am reading Using The foreach Package document and I have tried the following: - sessionInfo() R version 2.10.1 (2009-12-14) i386-pc-mingw32 locale: [1]

Re: [R] Parallel computing on Windows (foreach)

Here's what I do to get processors on my Windows machine (it should work for dual cores as well) fully engaged in a calculation (wow I like that) require(snow) require(doSNOW) require(foreach) #change the 8 to however many cores/phys processors you have on your machine cl.tmp =

Re: [R] Problem with the recode function

Dear Alain, -Original Message- From: Alain Guillet [mailto:alain.guil...@uclouvain.be] Sent: June-16-10 8:30 AM To: John Fox Cc: r-help@r-project.org Subject: Re: [R] Problem with the recode function Dear John, Thanks a lot for the time you spent on my problem. I don't believe

[R] replicate a grob n times

Dear all, I use the following to create a list of identical grobs, require(grid) rep.grob - function(g, n){ replicate(n, g, simplify=FALSE) } This approach suffers two problems: 1- R CMD check is not happy about the S3-like name. How can / Should I make this a real S3 method? 2- I don't know

[R] R software

Hello, I have downloaded several times your program R-2.11.1 32.exe, after the download when I click on the exe file I always receive the following error: The setup files are corrupted. Please obtain a new copy of the program. Could you help me to download English or Russian version for Windows

Re: [R] Backslash in paste() function

Just double all the backslashes and you are fine. In order to see the outcome, use cat() (not print). Uwe Ligges On 16.06.2010 09:49, Stefan Petersson wrote: Hi, I'm trying to build a vector of latex commands. However, I need the command strings to begin with a backslash \. I

Re: [R] working with zoo time index ??

On Wed, Jun 16, 2010 at 8:36 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Wed, Jun 16, 2010 at 7:58 AM, skan juanp...@gmail.com wrote: Hi thanks Let say data are written like this: 1990-01-01 10:01:00 ,  0.910 1990-01-01 10:03:00 ,  0.905 Would it be ok to read it with

Re: [R] an alternative to R for nonlinear stat models

I'd echo Paul's sentiments that we need to know where and when and how R goes slowly. I've been working on optimization in many environments for many years (try ?optim and you'll find me mentioned!). So has Dave Fournier. AD Model Builder has some real strengths and we need Automatic

[R] Problems loading RWeka and rJava under R 2.10.1

Dear R-users, I am running R 2.10.1 under Windows XP Professional. Moreover, I run Java 5.0 and the RWeka package 0.4-1 and rJava 0.8-4. are successfully installed. I get two error messages while loading both packages. library(RWeka) Error in get(Info[i, 1], envir = env) :

Re: [R] prcomp() and the lenght of PC:s

What do you mean by 'lenght'? It is part of the definition that the coefficient vector has Euclidean length one: a principal component is a projection. See for example MASS p.302. I don't see anything that has length 1 in the R sense. On Wed, 16 Jun 2010, Atte Tenkanen wrote: Hi, I would

[R] data frame

Dear list, I have the following problem. I have a data frame like this CLUSTERYEAR variableDelta R_pivot M1 2005 EC01 NA NA M1 2006 EC012

Re: [R] MANOVA proportion of variance explained

I think you are looking for a multivariate measure of association, analogous to R^2 for a univariate linear model. If so, there are extensions of eta^2 from univariate ANOVAs for each of the multivariate test statistics, e.g., for Pillai (-Bartlett) trace and Hotelling-Lawley trace and a given

Re: [R] R software

This is not a problem of R, but of your own computer. If it is a corrupted download, it might be stored in the temporary internet files. Delete all temporary internet files and try to download again. If that doesn't work, you'll have to dive into the Windows itself. That error can have numerous

Re: [R] generating samples from multivariate distributions

You can try the function Mvdc from the package copula, the function rmnorm from the package mnormt, or the function rmvnorm from the package mvtnorm. Or you can just type in Google R random multivariate normal distribution and get the same answer. Cheers Joris On Wed, Jun 16, 2010 at 1:41 PM,

[R] Latex and r

Dear R´ers I´m trying to get a summary table using latex and summary in the rms package to no avail. I´m running R 2.10.1, Mac OS X snow leopard and I have the mactex 2009 distribution installed. Any obvious things I´m missing? //M options(digits=3) set.seed(173) sex -

[R] Sampling with replacement

Dear all, good morning, I have a population, let say members are tagged with some simple number like 1,2,3,...20. I want to draw a sample with replacement of size 4 (say, can be more than 20 also). Is there any R function which will show me all such possible samples? Thanks

Re: [R] Sampling with replacement

sample(1:20,4,replace=TRUE) should do it. Jun On Wed, Jun 16, 2010 at 9:20 AM, Somnath Somnath somnath700...@gmail.comwrote: Dear all, good morning, I have a population, let say members are tagged with some simple number like 1,2,3,...20. I want to draw a sample with replacement of size 4

Re: [R] Sampling with replacement

Try sample(20, 4, replace = TRUE) HTH, Jorge On Wed, Jun 16, 2010 at 10:20 AM, Somnath Somnath wrote: Dear all, good morning, I have a population, let say members are tagged with some simple number like 1,2,3,...20. I want to draw a sample with replacement of size 4 (say, can be more

Re: [R] Question

Two possibilities : rescale your random vector, or resample to get numbers within the range. But neither of these solutions will give you a true exponential distribution. I am not aware of truncated exponential distributions that are available in R, but somebody else might know more about that. #

Re: [R] Latex and r

moleps wrote: Dear R´ers I´m trying to get a summary table using latex and summary in the rms package to no avail. I´m running R 2.10.1, Mac OS X snow leopard and I have the mactex 2009 distribution installed. Any obvious things I´m missing? file587f83cb.log. sh: xdvi: command not found

Re: [R] R software

On 06/16/2010 03:08 PM, Gulbanu Rys wrote: Hello, I have downloaded several times your program R-2.11.1 32.exe, after the download when I click on the exe file I always receive the following error: The setup files are corrupted. Please obtain a new copy of the program. Could you help me to

Re: [R] Multiple regressions

Hi, On Mon, Jun 14, 2010 at 4:24 PM, Amy Hessen amy_4_5...@hotmail.com wrote: Hi Steve, Thank you very much for your reply. Could you please tell me how I can know that a given dataset is non linear regression? This is the type of question you should really consult a (local) stat expert

[R] Mean variance plot of a data frame

Hello, I have a data frame that looks like Position, Relevance 0, 0.151531117072265 0, 0.245532564696541 0, 0.285207130323724 0, 0.302718099889669 0, 0.308213700400889 0, 0.341562706518953 0.0151515151515152, 0.368991215730364 0.0263157894736842, 0.256406702156839 0.0263157894736842,

Re: [R] Sampling with replacement

If you for some reason want to be shown all the possible combinations, try expand.grid(1:20,1:20,1:20,1:20) (ugly code). Don't use this for sampling. hth Rafael 2010/6/16 Jorge Ivan Velez jorgeivanve...@gmail.com Try sample(20, 4, replace = TRUE) HTH, Jorge On Wed, Jun 16, 2010 at

Re: [R] Sampling with replacement

On Jun 16, 2010, at 10:20 AM, Somnath Somnath wrote: Dear all, good morning, I have a population, let say members are tagged with some simple number like 1,2,3,...20. I want to draw a sample with replacement of size 4 (say, can be more than 20 also). Already answered on the list. Is

[R] error message stops for command

Hello all, I am using a for statement to run a loop in R but I get an error message and the loop stops. I need that the loop continues til the end. I know I can use the NEXT command but I don´t know how the sintax woul be for this particular situation. Any suggestions? Thanks, Diana Buitrago

Re: [R] Question about user define function

On 15/06/10 21:39, GL wrote: Have the following function that is called by the statement below. Trying to return the two dataframes, but instead get one large list including both tables. ReadInputDataFrames- function() { dbs.this= read.delim(this.txt, header = TRUE, sep = \t, quote=\,

Re: [R] error message stops for command

for (i in 1:100){ if (i%%5==0) {print(And provide reproducible code.);next} print(Before posting to the mailing list, I will read the posting guide.) } Cheers Joris On Wed, Jun 16, 2010 at 5:37 PM, diana buitrago dianacam...@gmail.com wrote: Hello all, I am using a for statement to run a

Re: [R] Sampling with replacement

Hi Rafael, You might try: r - expand.grid(rep(list(1:20), 4)) dim(r) [1] 16 4 HTH, Jorge 2010/6/16 Rafael Björk If you for some reason want to be shown all the possible combinations, try expand.grid(1:20,1:20,1:20,1:20) (ugly code). Don't use this for sampling. hth Rafael

[R] RegExp question

Dear all, I'm trying to filter out the number of leaves (it should be 1 in the example below) from the following string: string [1] Java-Object{J48 pruned tree\n--\n: 0 (15.0/3.0)\n \nNumber of Leaves : \t1\n\nSize of the tree : \t1\n} Any idea how to do that as simple as

[R] Is there any function to read genbank flat files.

Hi I am trying to read genbank est flat files. Is there any function in R to do that. I was looking for something similar to SeqIO.parse(input_handle, genbank) in bioPython that convert the genebank sequence into fasta sequence. (One option is to call biopython from R but I am not familiar with

Re: [R] Is there any function to read genbank flat files.

Hi, I think it's better to ask on Bioconductor mailing list. Regards On Jun 16, 6:19 pm, Fahim Md fahim...@gmail.com wrote: Hi I am trying to read genbank est flat files. Is there any function in R to do that. I was looking for something similar to SeqIO.parse(input_handle, genbank)  in

Re: [R] Mean variance plot of a data frame

Hi: Look at the function plotmeans() in package gplots and run example(plotmeans). I tried something like plotmeans(Relevance ~ round(Position, 3), data = df, xlab = 'Position') but you may want to play with some of the options to get what you want. HTH, Dennis On Wed, Jun 16, 2010 at 8:25 AM,

Re: [R] RegExp question

On Jun 16, 2010, at 12:04 PM, Andrej wrote: Dear all, I'm trying to filter out the number of leaves (it should be 1 in the example below) from the following string: string [1] Java-Object{J48 pruned tree\n--\n: 0 (15.0/3.0)\n \nNumber of Leaves : \t1\n\nSize of the tree :

Re: [R] RegExp question

Thanks David for your fast reply, but now I realized tat string is of type: class(string) [1] jobjRef attr(,package) [1] rJava so I get an error when i try with gsub or sub: sub(^.+\\t(\\d+)\\n.+$, \\1, string) Error in as.character.default(x) : no method for coercing this S4 class to a

Re: [R] RegExp question

Publicly produce something we can work with. I have no idea how to create an example that will match such an object. ?dput ?dump Read Posting Guide. -- David. On Jun 16, 2010, at 12:54 PM, Andrej wrote: Thanks David for your fast reply, but now I realized tat string is of type:

Re: [R] error message stops for command

Hi, On Wed, Jun 16, 2010 at 11:37 AM, diana buitrago dianacam...@gmail.com wrote: Hello all, I am using a for statement to run a loop in R but I get an error message and the loop stops. I need that the loop continues til the end. I know I can use the NEXT command but I don´t know how the

Re: [R] RegExp question

Sorry, I apologize. Below is the minimal example. library(RWeka) model - J48(as.factor(Species)~., data = iris) model J48 pruned tree -- Petal.Width = 0.6: setosa (50.0) Petal.Width 0.6 | Petal.Width = 1.7 | | Petal.Length = 4.9: versicolor (48.0/1.0) | | Petal.Length

[R] Is there a non-parametric repeated-measures Anova in R ?

Hello Prof. Harrell and dear R-help mailing list, I wish to perform a non-parametric repeated measures anova. If what I read online is true, this could be achieved using a mixed Ordinal Regression model (a.k.a: Proportional Odds Model). I found two packages that seems relevant, but couldn't find

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

It's possible to use the ordinal regression model if your data are ordered categories. The standard non-parametric test is the Friedman test. ?friedman.test Jeremy On 16 June 2010 10:22, Tal Galili tal.gal...@gmail.com wrote: Hello Prof. Harrell and dear R-help mailing list, I wish to

Re: [R] Parallel computing on Windows (foreach) (Sergey Goriatchev)

foreach (or virtually anything you might use for concurrent programming) only really makes sense if the work the clients are doing is substantial enough to overwhelm the communication overhead. And there are many ways to accomplish the same task more or less efficiently (for example, doing blocks

Re: [R] Latex and r

moleps wrote: xdvi is installed in the same location as yours. I even did a reinstallment of mactex. Still doesnt work. But since I´m now convinced its related to my latex distribution I´ll take the problem elsewhere.. And the directory that xdvi is located in is in your path?

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

Hello Jeremy, Thank you for replying. I came across friedman test (I even wrote and published R code to easily perform a post-hoc analysis of friedman testhttp://www.r-statistics.com/2010/02/post-hoc-analysis-for-friedmans-test-r-code/ ). But what I am after is *multi-way* repeated-measures

Re: [R] RegExp question

On Jun 16, 2010, at 1:05 PM, Andrej wrote: Sorry, I apologize. Below is the minimal example. library(RWeka) model - J48(as.factor(Species)~., data = iris) model J48 pruned tree -- Petal.Width = 0.6: setosa (50.0) Petal.Width 0.6 | Petal.Width = 1.7 | | Petal.Length =

Re: [R] Latex and r

moleps wrote: Dear R´ers I´m trying to get a summary table using latex and summary in the rms package to no avail. I´m running R 2.10.1, Mac OS X snow leopard and I have the mactex 2009 distribution installed. Any obvious things I´m missing? //M options(digits=3) set.seed(173) sex -

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Tal Galili Sent: Wednesday, June 16, 2010 10:44 AM To: Jeremy Miles Cc: r-help@r-project.org Subject: Re: [R] Is there a non-parametric repeated-measures Anova in R ? Hello

[R] Reading data from xls..........please help

Can anyone help me how to read xls file into R. I have tried following library(gdata) xlsfile - file.path(.path.package('gdata'),'xls','iris.xls') read.xls(xlsfile) I got following error: Converting xls file to csv file... Error in system(cmd, intern = !verbose) : perl not found Error in

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

The Friedman test lacks power. When there are only 2 blocks it reduces to the inefficient sign test. Frank On 06/16/2010 12:43 PM, Tal Galili wrote: Hello Jeremy, Thank you for replying. I came across friedman test (I even wrote and published R code to easily perform a post-hoc analysis of

[R] Antwort: Re: an alternative to R for nonlinear stat models

Hi I think my first answer doesn't seem to have gone through due to the attachement. So I copy pasted my last answer including the code below. Sorry about that: Hi Paul Oh ok sorry, I send you below a copy of the R code I was using. It's very possible that my

Re: [R] Parallel computing on Windows (foreach) (Sergey Goriatchev)

Hi, I think I lost the reference email somewhere, but: Not only is the sequential foreach much slower than the simple for-loop (as least in this particular instance), but I am not quite sure how to make foreach run parallel. Where would I get this parallel backend? I looked at doMC and

Re: [R] an alternative to R for nonlinear stat models

As far as I can tell, Gove et al. (2002) might be a good example for benchmarking the optimization performance of R vs. ADMB. It would be great if expert R users/developers could tweak Beni's model so that the performance comparison is valid. The main purpose is not to see which is faster or

[R] Help asked for automated generation of ncdf variables

Hi, I'm using the ncdf library for creating ncdf files. But I need to create about 100 variables per file (e.g. single rivers), So I do not like to create each variable separately. Unfortunately I found no way to make this work, as I'm unable to create a correct list of class var.ncdf. ###code

Re: [R] [ADMB Users] an alternative to R for nonlinear stat models

Hi Arni (and others), My dissertation work involves use (and extension) of models of the same ilk (sometimes exactly the same) as those described by Nancy Gove and John Skalski in their 2002 article. I began with R, and moved to my own home-brewed C/C++ programs for the sake of of speed when

Re: [R] Latex and r

moleps mole...@gmail.com writes: Apparently you don't have xdvi installed on your system. HTH Georg Dear R´ers I´m trying to get a summary table using latex and summary in the rms package to no avail. I´m running R 2.10.1, Mac OS X snow leopard and I have the mactex 2009 distribution

Re: [R] Sampling with replacement

How about library(TeachingSampling) SupportWR(20,4) Tom -- View this message in context: http://r.789695.n4.nabble.com/Sampling-with-replacement-tp2257450p2257644.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] working with zoo time index ??

Hi I'll ask in a different way... I have all this in a file.txt 1990-01-01 10:00:00 , 0.900 # element 1 1990-01-01 10:01:00 , 0.910 # element 2 1990-01-01 10:03:00 , 0.905 # element 3 1990-01-01 10:04:00 , 0.905 # element 4 1990-01-01 10:05:00 ,

[R] Function argument as string

Hi, Suppose a write a function a_fn-function(arg1) { return(table(arg1)); } I have a column called AGE. Now I call the function c = a_fn(AGE); When a_fn is called, AGE is received in arg1. My question is, how do I access the actual name of the argument arg1? i.e, inside the function, i

[R] Accessing the elements of summary(prcomp(USArrests))

Hello again, I was hoping one of you could help me with this problem. Consider the sample data from R: summary(prcomp(USArrests)) Importance of components: PC1 PC2PC3 PC4 Standard deviation 83.732 14.2124 6.4894 2.48279 Proportion of Variance 0.966

Re: [R] working with zoo time index ??

I said taking the first element everyday, but that was just an example, I could need one every 2 hours or something more complicated such as one every hour if the former one was non null. -- View this message in context:

Re: [R] Sampling with replacement

Thanks for all those reply. Is there any general rule to determine how many samples I would get from a population of size n, I draw a sample of size m (m may be greater than n) if sample is drawn with replacement? Thanks, 2010/6/16 Jorge Ivan Velez jorgeivanve...@gmail.com Hi Rafael, You

Re: [R] Function argument as string

Vishwanath Sindagi wrote: Hi, Suppose a write a function a_fn-function(arg1) { return(table(arg1)); } I have a column called AGE. Now I call the function c = a_fn(AGE); When a_fn is called, AGE is received in arg1. My question is, how do I access the actual name of the argument arg1?

Re: [R] Latex and r

xdvi is installed in the same location as yours. I even did a reinstallment of mactex. Still doesnt work. But since I´m now convinced its related to my latex distribution I´ll take the problem elsewhere.. Regards, //M On 16. juni 2010, at 17.19, Prof Brian Ripley wrote: On Wed, 16 Jun

Re: [R] Reading data from xls..........please help

On Wed, Jun 16, 2010 at 2:29 PM, Christofer Bogaso bogaso.christo...@gmail.com wrote: Can anyone help me how to read xls file into R. I have tried following library(gdata) xlsfile - file.path(.path.package('gdata'),'xls','iris.xls') read.xls(xlsfile) I got following error: Converting xls

[R] mgcv, testing gamm vs lme, which degrees of freedom?

Dear all, I am using the mgcv package by Simon Wood to estimate an additive mixed model in which I assume normal distribution for the residuals. I would like to test this model vs a standard parametric mixed model, such as the ones which are possible to estimate with lme. Since the smoothing

  1   2   >