Re: [R] Fwd: about plantbreeding library

2013-08-17 Thread Jeff Newmiller
Please keep in mind that this is the R-help mailing list, not the plantbreeding-help mailing list, so most (perhaps all) of us don't know what the theory behind your work is. The fact that you are attempting to load data sets and use functions that are not accessible yet you are forging on

Re: [R] latin1 encoding in WriteXLS

2013-08-17 Thread Hugo Varet
Yes, it also occurs with WriteXLS version 3.2.1. This test on several computers always leads to the same error. Hugo Varet 2013/8/17 Rainer Hurling rhur...@gwdg.de Am 13.08.2013 19:40, schrieb Hugo Varet: Dear R users, I've just updated the WriteXLS package (on R 3.0.1) and I now have

[R] Amelia, Zelig and latex in R

2013-08-17 Thread Francesco Sarracino
Dear listers, I am running some OLS on multiply imputed data using Amelia. I first imputed the data with Amelia. than I run a OLS using Zelig to obtain a table of results accounting for the multiply imputed data-sets. And I'd like to do this for various models. Finally, I want to output all the

Re: [R] is there anything wrong with doing a multinomial logistic like it is a binary logistic regression?

2013-08-17 Thread Michael Dewey
At 04:13 17/08/2013, Eunice Chou wrote: My outcome variable (y) is 3 categories. Is there anything bad about using the following code to get a parameter estimate for my bivariate model? publicfit = glm(y ~ public, data=dataSPSS.vmj, family=binomial) Have you tried it? What did R tell you?

Re: [R] Is it possible to avoid copying arrays when calling list()?

2013-08-17 Thread MRipley
Thanks for the input, but it looks like I found a simple solution. Turns out that if you assign to lists by name, then R doesn't make extra copies: x-double(10^9) mylist-list() system.time(mylist[[1]]-x) user system elapsed 2.992 3.352 6.364 x-double(10^9) mylist-list()

Re: [R] Amelia, Zelig and latex in R

2013-08-17 Thread Francesco Sarracino
Hi Christopher, thanks for your reply. Unfortunately, that's not what I am looking for. I would like to have a table with the results of the two models (lm.imputed1 and lm.imputed2) in two separate columns. According to stargazer syntax I should type something like: stargazer(lm.imputed1,

Re: [R] Amelia, Zelig and latex in R

2013-08-17 Thread Christopher Desjardins
Does this do what you want? library(Amelia) library(Zelig) library(stargazer) library(xtable) data(africa) m = 10 imp1 - amelia(x = africa,cs=country,m=m) imp2 - amelia(x = africa,cs=country,m=m) lm.imputed1 - zelig(infl ~ trade + civlib, model=ls,data = imp1) lm.imputed2 - zelig(infl ~ trade +

Re: [R] Error in Corpus() in tm package

2013-08-17 Thread Milan Bouchet-Valat
Le vendredi 16 août 2013 à 19:35 -0700, Ajinkya Kale a écrit : I am trying to use the text mining package ... I keep getting this error : rm(list=ls()) library(tm) sourceDir - Z:\\projectk_viz\\docs_to_index ovid - Corpus(DirSource(sourceDir),readerControl = list(language = lat)) Error

Re: [R] Amelia, Zelig and latex in R

2013-08-17 Thread Christopher Desjardins
What do you mean by results? Do you want just the estimated parameters? And are you looking for one big table with all the estimated parameters from all imputation runs? Chris On Sat, Aug 17, 2013 at 11:18 AM, Francesco Sarracino f.sarrac...@gmail.com wrote: Hi Christopher, thanks for your

Re: [R] Amelia, Zelig and latex in R

2013-08-17 Thread Christopher Desjardins
Oh and are you looking for just the summarized results over all the imputed runs? i thought you wanted them from each iteration. On Sat, Aug 17, 2013 at 11:38 AM, Christopher Desjardins cddesjard...@gmail.com wrote: What do you mean by results? Do you want just the estimated parameters? And

Re: [R] regex challenge

2013-08-17 Thread Frank Harrell
Bill I found a workaround: f - ff(formula, lab) f - as.formula(gsub(`, , as.character(deparse(f Thanks for your elegant solution. Frank -- Thanks Bill. The problem is one of the results of convertName might be 'Heading(Age in Years)*age' (this is for the

Re: [R] Error in Corpus() in tm package

2013-08-17 Thread Ajinkya Kale
It contains all text files which were converted from doc, docx, ppt etc. using libreoffice. Some of them are non-english text documents. Sorry I cannot share the corpus.. but if someone can shed light on what might cause this error then I can try to eliminate those documents if some specific docs

[R] Runtime error using ncdf

2013-08-17 Thread Camilo Mora
HI everyone, I have encountered a problem while using ncdf to open nc files in R. I found in the internet several comments in the past but no solution. I could not find a direct solution but I found the source of the problem, if anyone may know where the solution could be, and an

Re: [R] Runtime error using ncdf

2013-08-17 Thread David W. Pierce
Hi Camilo, you don't say what platform you are running on, but the version of the underlying netcdf library installed on your machine must have been properly compiled to enable access to files greater than 2GB. Otherwise, the R interface to the netcdf library cannot work with such files. You

Re: [R] Error in Corpus() in tm package

2013-08-17 Thread Ajinkya Kale
Funny, it works fine if I use VectorSource ovid - Corpus(VectorSource(list.files(sourceDir)[1:1253]), readerControl = list(language = lat)) So I tried only executing DirDource(sourceDir) and that fails with the error i mentioned earlier. So its not a problem with Corpus() which I thought

Re: [R] Error in Corpus() in tm package

2013-08-17 Thread Ajinkya Kale
I think I know why it works faster, cause VectorSource in above code only takes the files names as a corpus and not the contents of the files :D duh! Any suggestions to create a vector source out of contents of the txt files ? On Sat, Aug 17, 2013 at 1:59 PM, Ajinkya Kale kaleajin...@gmail.com

[R] remove from list

2013-08-17 Thread Tony Paredes
[[alternative HTML version deleted]] __ 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/posting-guide.html and provide commented, minimal, self-contained,

[R] Appropriateness of R functions for multicore

2013-08-17 Thread Hopkins, Bill
Has there been any systematic evaluation of which core R functions are safe for use with multicore? Of current interest, I have tried calling read.table() via mclapply() to more quickly read in hundreds of raw data files (I have a 24 core system with 72 GB running Ubuntu, a perfect platform for

[R] A function taken out of its original environment performs more slowly.

2013-08-17 Thread Xiao He
Hi dear R-users, I encountered an interesting pattern. Take for example the function combn(), I copied and pasted the function definition and saved it as a new function named combn2() (see the end of this email). As it turned out, combn2() seems to be substantially slower than the original

Re: [R] A function taken out of its original environment performs more slowly.

2013-08-17 Thread Uwe Ligges
On 18.08.2013 01:05, Xiao He wrote: Hi dear R-users, I encountered an interesting pattern. Take for example the function combn(), I copied and pasted the function definition and saved it as a new function named combn2() (see the end of this email). As it turned out, combn2() seems to be

Re: [R] Appropriateness of R functions for multicore

2013-08-17 Thread Jeff Newmiller
In most threaded multitasking environments it is not safe to perform IO in multiple threads. In general you will have difficulty performing IO in parallel processing so it is best to let the master hand out data to worker tasks and gather results from them for storage. Keep in mind that just

Re: [R] Fwd: about plantbreeding library

2013-08-17 Thread Kevin Wright
If you have properly installed the plantbreeding package, then each time you start R, you need to type: library(plantbreeding) before you can access the data or the functions in the package. Kevin On Sat, Aug 17, 2013 at 12:18 AM, Waqas Shafqat waqas1...@gmail.com wrote: -- Forwarded

Re: [R] remove from list

2013-08-17 Thread David Winsemius
Please read the details offered in the resources below. On Aug 17, 2013, at 2:27 PM, Tony Paredes wrote: [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] remove from list

2013-08-17 Thread bioprogrammer
Sent from my T-Mobile 4G LTE Device Original message From: David Winsemius dwinsem...@comcast.net Date: 08/17/2013 9:39 PM (GMT-07:00) To: Tony Paredes tgal...@gmail.com Cc: r-help@r-project.org Subject: Re: [R] remove from list Please read the details offered in the

[R] First time r user

2013-08-17 Thread Dylan Doyle
Hello R users, I have recently begun a project to analyze a large data set of approximately 1.5 million rows it also has 9 columns. My objective consists of locating particular subsets within this data ie. take all rows with the same column 9 and perform a function on that subset. It was