Re: [R] Is it possible to define another kind of NA

2014-11-13 Thread Marc Girondot
Le 13/11/2014 01:26, MacQueen, Don a écrit : Along the lines of what Bert Gunter said, the ideal way to represent I agree that LDL is a special case of what could be named ODL (Out of detection limit). To answer to Bert Gunter, indeed if LDL (or ODL) values are changed into NA, the results will

[R] metafor - code for analysing geometric means

2014-11-13 Thread Purssell, Ed
?Dear All I have some data expressed in geometric means and 95% confidence intervals. Can I code them in metafor as: rma(m1i=geometric mean 1, m2i=geometric mean 2, sd1i=geometric mean 1 CI /3.92, sd2i=geometric mean 2 CI/3.92...etc, measure="MD") All of the studies use geometric means.

Re: [R] metafor - code for analysing geometric means

2014-11-13 Thread Michael Dewey
On 13/11/2014 11:00, Purssell, Ed wrote: ?Dear All I have some data expressed in geometric means and 95% confidence intervals. Can I code them in metafor as: rma(m1i=geometric mean 1, m2i=geometric mean 2, sd1i=geometric mean 1 CI /3.92, sd2i=geometric mean 2 CI/3.92...etc, measure="

Re: [R] Is it possible to define another kind of NA

2014-11-13 Thread Keith Jewell
On 13/11/2014 11:08, Marc Girondot wrote: Le 13/11/2014 01:26, MacQueen, Don a écrit : Along the lines of what Bert Gunter said, the ideal way to represent I agree that LDL is a special case of what could be named ODL (Out of detection limit). To answer to Bert Gunter, indeed if LDL (or ODL) va

Re: [R] Is it possible to define another kind of NA

2014-11-13 Thread peter dalgaard
As you may know already, that design can be considerably improved by making the wheel triangular, reducing the bump count per revolution by a full 25% -pd On 13 Nov 2014, at 16:17 , Keith Jewell wrote: > You are re-inventing the wheel and yours will probably end up square! -- Peter Dalgaard

Re: [R] Is it possible to define another kind of NA

2014-11-13 Thread Rich Shepard
On Thu, 13 Nov 2014, peter dalgaard wrote: As you may know already, that design can be considerably improved by making the wheel triangular, reducing the bump count per revolution by a full 25% It has been written that those who go around in circles think of themselves as big wheels. Rich

Re: [R] Is it possible to define another kind of NA

2014-11-13 Thread Rich Shepard
On Thu, 13 Nov 2014, Keith Jewell wrote: You are re-inventing the wheel and yours will probably end up square! R already has facilities for handling censored data, e.g. Surv in the survival package (which despite its name is applicable to applications other than survival analysis). There i

[R] How to put inlined C code on a worker node?

2014-11-13 Thread Mikhail Umorin
Hello — I am using inline C functions within foreach %dopar% loop. On SMP (doParallel, doMP) it works but on MPI-based clusters (doMPI) it does’t. The reason, I think, is because the object code produced using the inline package, essentially an .so file, does not get copied onto the worker node

[R] number of weights in multinom

2014-11-13 Thread Franck Vermet
Dear colleagues, In the function multinom (package nnet), I get the following message after training for a model with 9 inputs and 6 classes (output) : # weights: 66 (50 variable) I understand that there are 50 variables in the model, but I don't understand the number 66. How can we interpret

[R] hi

2014-11-13 Thread Cox Lwaka
I have a bit of trouble here to program in r. I am anew user but i really enjoy working with it.I have a large number of variables in a matrix that are arranged sequentially on a line (chromosome). This order has to be maintained whatsoever. I am to develop an r algorithm that will develop group

[R] merge 2 data.frames in dependence of 2 values

2014-11-13 Thread Matthias Weber
Hello togehter, i have a little problem. Maybe anyone can help me. I have 2 data.frames, which look like as follows: First: NAMEMONTH BONUS 1 Andy 2014-10 100 2 Pete 2014-10200 3 Marc2014-10300 4 Andy2014-11400 Se

Re: [R] merge 2 data.frames in dependence of 2 values

2014-11-13 Thread Rui Barradas
Hello, See ?merge, in particular the argument 'all'. dat1 <- read.table(text = " NAMEMONTH BONUS 1 Andy 2014-10 100 2 Pete 2014-10200 3 Marc2014-10300 4 Andy2014-11400 ", header = TRUE, stringsAsFactors = FALSE)

Re: [R] merge 2 data.frames in dependence of 2 values

2014-11-13 Thread William Dunlap
merge(df1, df2, all=TRUE) Bill Dunlap TIBCO Software wdunlap tibco.com On Thu, Nov 13, 2014 at 6:02 AM, Matthias Weber < matthias.we...@fntsoftware.com> wrote: > Hello togehter, > > i have a little problem. Maybe anyone can help me. > > I have 2 data.frames, which look like as follows: > First:

Re: [R] merge 2 data.frames in dependence of 2 values

2014-11-13 Thread Michael Dewey
On 13/11/2014 14:02, Matthias Weber wrote: Hello togehter, i have a little problem. Maybe anyone can help me. I think you might find ?merge enlightening Indeed given that the word merge occurs in your subject line and your text it is surprising you have not already found it. I have 2 d

[R] rpart package: How can I save print(rpart)

2014-11-13 Thread Karim Mezhoud
Hi, the print of rpart fitting gives the summary of tree I would like to save the console text of: fit <- rpart(formula, data) summary <- print(fit) when I look in "summary" I did not find the same thing as in "print(rpart)" [1] "Clinical Data exists" [1] "merging samples from Clinical an

[R] rpart package: prp does not work with RStudio and R console

2014-11-13 Thread Karim Mezhoud
Hi, All right for rpart package but it seems there is a confusion in "text" function. prp is not found Error in fun() : could not find function "prp" Thanks Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ [[alternative HTML version deleted]] ___

Re: [R] rpart package: How can I save print(rpart)

2014-11-13 Thread William Dunlap
Use capture.output(), as in > junk <- capture.output(summary(1:10)) > junk [1] " Min. 1st Qu. MedianMean 3rd Qu.Max. " [2] " 1.003.255.505.507.75 10.00 " > cat(junk, sep="\n") Min. 1st Qu. MedianMean 3rd Qu.Max. 1.003.255.505.5

Re: [R] rpart package: How can I save print(rpart)

2014-11-13 Thread Karim Mezhoud
Yes Thanks! that works, but I loose the \n when I would like to save or edit it. getTextInWindows is a function that edits any text in editor. getTextInWindows(summary): without "\n" save (file= "junk.txt", junk):without "\n" getTextInWindow(capture.output(cat(junk, sep = "\n"))) :No works Thank

Re: [R] rpart package: How can I save print(rpart)

2014-11-13 Thread William Dunlap
Use paste(collapse="\n", junk) if you want it as a single string with \n's in it. Bill Dunlap TIBCO Software wdunlap tibco.com On Thu, Nov 13, 2014 at 1:13 PM, Karim Mezhoud wrote: > Yes Thanks! that works, > but I loose the \n when I would like to save or edit it. > > getTextInWindows is a fun

Re: [R] rpart package: How can I save print(rpart)

2014-11-13 Thread Karim Mezhoud
OK thanks That works. :) Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/ On Thu, Nov 13, 2014 at 10:41 PM, William Dunlap wrote: > Use paste(collapse="\n", junk) if you want it as a single string with \n's > in it. > > Bill Dunlap > TIBCO Software > wdunlap tibco

[R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread Zheng Da
Hello, I'm trying to generate plain text from the .Rd files. I run "R CMD Rd2txt xx.Rd -o xx.txt" to get the plain text file. When I open the text file, I see "_^H" in front of each character in all section titles. I can't figure out where the problem is. If I don't specify -o argument, "R CMD Rd2

[R] Help with ddply/summarize

2014-11-13 Thread John Posner
I have a straightforward application of ddply() and summarize(): ddply(MyFrame, .(Treatment, Week), summarize, MeanValue=mean(MyVar)) This works just fine: Treatment Week MeanValue 1MyDrug BASELINE 5.91 2MyDrugWEEK 1 4.68 3MyDrugWEEK 2 4.08 4MyDr

[R] Data Import to R

2014-11-13 Thread Ramesh Gautam
While importing .csv files into R, all data are converted to factor-by default. But, how can I preserve the original format of the data like numeric to numeric, integer to integer, character to character etc while importing from csv to R environment. I tried several ways, no thing helps. I used 's

Re: [R] Help with ddply/summarize

2014-11-13 Thread David Winsemius
On Nov 13, 2014, at 3:31 PM, John Posner wrote: > I have a straightforward application of ddply() and summarize(): > > ddply(MyFrame, .(Treatment, Week), summarize, MeanValue=mean(MyVar)) > > This works just fine: > > Treatment Week MeanValue > 1MyDrug BASELINE 5.91 > 2My

Re: [R] Data Import to R

2014-11-13 Thread David Winsemius
On Nov 13, 2014, at 4:28 PM, Ramesh Gautam wrote: > While importing .csv files into R, all data are converted to factor-by > default. But, how can I preserve the original format of the data like > numeric to numeric, integer to integer, character to character etc while > importing from csv to R e

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread Duncan Murdoch
On 13/11/2014, 2:51 PM, Zheng Da wrote: > Hello, > > I'm trying to generate plain text from the .Rd files. I run "R CMD > Rd2txt xx.Rd -o xx.txt" to get the plain text file. When I open the > text file, I see "_^H" in front of each character in all section > titles. I can't figure out where the pr

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread Rolf Turner
On 14/11/14 13:59, Duncan Murdoch wrote: On 13/11/2014, 2:51 PM, Zheng Da wrote: Hello, I'm trying to generate plain text from the .Rd files. I run "R CMD Rd2txt xx.Rd -o xx.txt" to get the plain text file. When I open the text file, I see "_^H" in front of each character in all section titles.

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread David Winsemius
On Nov 13, 2014, at 11:51 AM, Zheng Da wrote: > Hello, > > I'm trying to generate plain text from the .Rd files. I run "R CMD > Rd2txt xx.Rd -o xx.txt" to get the plain text file. When I open the > text file, I see "_^H" in front of each character in all section > titles. I can't figure out wher

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread Zheng Da
Thank you. Great help! Given your information, we can generate text files without _^H with the following command: echo "tools::Rd2txt(\"$input_file\", out=\"$output_file\", options=list(underline_titles=FALSE))" | R --no-save It's kind tedious though. I'm just curious. Which text editor can actual

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread Jeff Newmiller
> Which text editor can actually display _^H correct? Not aware of any. It is from very old line printer behavior, which most (all?) printers can still support even though few computers are set up to utilize it. --- Jeff Newm

[R] HELP ON NON-LINEAR MIXED MODEL

2014-11-13 Thread Popoola Daniel
Good Morning Sir/Ma, I am POPOOLA DANIEL a Forest Biometrician in making from the University of Ibadan, Ibadan Nigeria. Please Sir/Ma I am having issues on performing Non-linear mixed model on R (using maximum likelihood approach). I am trying to input four different measured variables which are

Re: [R] Data Import to R

2014-11-13 Thread PIKAL Petr
Hi Your original numeric data probably contain something which prevents read.* to accept them as numeric (decimal point, white space) what is result of str(imported.data) Petr Pikal > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On B

Re: [R] hi

2014-11-13 Thread PIKAL Petr
Hi 1. Do not post in html 2. Post example of data, preferably by dput function 3. Based on posted data explain what do you want to achieve 4. As it seems to be biological issue did you look at Bioconductor? Cheers Petr > -Original Message- > From: r-help-boun...@r-