Re: [R] YSI data in R

2013-10-02 Thread Jim Lemon
On 10/03/2013 12:13 PM, Orlen wrote: My apologies. I want to test for normality to determine whether I need to use parametric or nonparametric tests to analyze the data. A prof suggested this was a good place to start with my data. I have attached a copy of part of my data. There are five tra

Re: [R] Command line r

2013-10-02 Thread Tomek R
Hi, I wanted to create even something simple-r than little-r or Rscript. I even use Rscript in simple-r (but with clunky github interface it is difficult to decipher it). The idea is that the first command that I always end up executing in R is d<-read.table('in.txt', ...). Then, I either plot t

Re: [R] YSI data in R

2013-10-02 Thread Orlen
My apologies. I want to test for normality to determine whether I need to use parametric or nonparametric tests to analyze the data. A prof suggested this was a good place to start with my data. I have attached a copy of part of my data. There are five transects, and five weeks, and each transe

Re: [R] YSI data in R

2013-10-02 Thread Bert Gunter
... and I would suggest that you read and follow The Posting Guide. -- Bert On Wed, Oct 2, 2013 at 6:44 PM, Rolf Turner wrote: > On 10/03/13 11:51, Orlen wrote: >> >> I am dealing with a huge data set collected from an YSI water monitoring >> probe - the data includes water depth, transect locat

Re: [R] YSI data in R

2013-10-02 Thread Rolf Turner
On 10/03/13 11:51, Orlen wrote: I am dealing with a huge data set collected from an YSI water monitoring probe - the data includes water depth, transect location, temperature, conductivity, salinity, pressure etc etc. I am trying to do basic summary statistics on my data and analyze the data. M

[R] YSI data in R

2013-10-02 Thread Orlen
I am dealing with a huge data set collected from an YSI water monitoring probe - the data includes water depth, transect location, temperature, conductivity, salinity, pressure etc etc. I am trying to do basic summary statistics on my data and analyze the data. My data is organized by water dep

Re: [R] foreign package problem with stata 13

2013-10-02 Thread Pablo Menese Camargo
The thing is, the oldest version in wich stata 13 save, is stata12. I can use this, but would be a good momento to upgrade foreign package. thanks a lot anyways!!! On Wed, Oct 2, 2013 at 8:04 PM, Richard Kwock wrote: > Use "saveold" instead of "save" when you save your dataset in stata. > > I do

Re: [R] foreign package problem with stata 13

2013-10-02 Thread Richard Kwock
Use "saveold" instead of "save" when you save your dataset in stata. I don't think foreign library in R supports stata 13 yet. Richard On Wed, Oct 2, 2013 at 3:11 PM, Pablo Menese Camargo wrote: > foreign package does not support dataset saved at stata 13. > anyone knows any wayt to make it wor

Re: [R] rJava does not see the htmlunit java package

2013-10-02 Thread Joe Byers
Joe Byers aim.com> writes: Sorry I forgot a subject. This thread is for rJava not seeing the htmlunit java package __ 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.o

[R] foreign package problem with stata 13

2013-10-02 Thread Pablo Menese Camargo
foreign package does not support dataset saved at stata 13. anyone knows any wayt to make it works? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Drawing garbled

2013-10-02 Thread Jim Lemon
On 10/03/2013 01:51 AM, yy wrote: Hi: I am Chinese, I am developing a java application, and deploy it to tomcat 7.0.42, I use rJava to use the R, when I use the command line to start tomcat, the R drawing well done, see attachment histview.png,but when i use windows service to start

Re: [R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread David Winsemius
On Oct 2, 2013, at 1:33 PM, Filipe Correia wrote: > David, > > On Wed, Oct 2, 2013 at 6:56 PM, David Winsemius > wrote: > This is not an appropriate question for R-help. You should ask your > statistics instructor or pose the question in a venue where explanations of > purely statistical pro

[R] fminsearch usage

2013-10-02 Thread Edward Patzelt
R Help - I'm attempting to use fminsearch and it continues to give me errors. I have a normal Rescorla Wagner function that has 2 parameters and returns total error. for some reason fminsearch keeps saying the type of errors below, I've put the function and data after it using dput > sol <- fmins

Re: [R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread Filipe Correia
David, On Wed, Oct 2, 2013 at 6:56 PM, David Winsemius wrote: > > This is not an appropriate question for R-help. You should ask your > statistics instructor or pose the question in a venue where explanations of > purely statistical problems are on-topic. (Just because you used R to get > results

Re: [R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread Filipe Correia
Hi, Thank you for your answer Charles. On Wed, Oct 2, 2013 at 6:50 PM, Charles Determan Jr wrote: > When you chose a different alternative argument you are asking a different > null hypothesis. I was considering the null hypothesis was expressed through the "mu" argument. I think I got this ide

Re: [R] R Excel - Microsoft Excel is waiting for OLE action

2013-10-02 Thread Jeff Newmiller
There is a mailing list devoted to supporting RExcel that you should be posting on at http://rcom.univie.ac.at/. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. #

Re: [R] Drawing garbled

2013-10-02 Thread Jeff Newmiller
Off topic here, this list is about R and your question is about system configuration. FWIW a common problem with execution of programs as services is permissions issues, but I have no desire to debug such OS-specific issues for you. BTW, both images look the same to me. -

Re: [R] ggplot2 multiple legends modification (resent witout attachments)

2013-10-02 Thread David Winsemius
On Oct 2, 2013, at 11:12 AM, rolf.kem...@renesas.com wrote: > > I resent this as I recognized that attachments are not allowed. Attachments are allowed, but only certain types . To attach data it should be MIME-text and the usual way to get attachments recognized as such jointly by your mail-

[R] Prevalence Assessment in R (Free Webinar)

2013-10-02 Thread Ray DiGiacomo, Jr.
Hello R Users, Are you interested in using R to perform prevalence assessment studies? This free one-hour webinar will provide a general introduction to "prevalence", an R package that provides a collection of functions used to perform prevalence assessment studies. Webinar Format: - Introductio

[R] ggplot2 multiple legends modification (resent witout attachments)

2013-10-02 Thread rolf . kemper
I resent this as I recognized that attachments are not allowed. Dear ALL, please find attached a small test file. (code is copied in this mail too) Basically (as shown) I would like to combine several results in one graph. I could succeed to plot the data correctly, but the legend generation

Re: [R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread David Winsemius
On Oct 2, 2013, at 8:33 AM, Filipe Correia wrote: > Hello everyone, > > I'm having some trouble interpreting the results of a Wilcoxon > (Mann-Whitney U) test. Hope you can help. > > This is the R script that I am running: > > a <- c(1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 5, 1, 1, 1, 3, 1,

Re: [R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread Charles Determan Jr
Filipe, When you chose a different alternative argument you are asking a different null hypothesis. You are looking at a two-tailed, lesser than, and greater than hypotheses. Which one you chose is dependent upon your initial question. Are you asking generically if your two populations (a and b

Re: [R] Bioconductor / AnnotationDbi: Why does a GOAllFrame contain more rows than its argument GoFrame?

2013-10-02 Thread Martin Morgan
On 10/02/2013 09:28 AM, Asis Hallab wrote: Dear Bioconductor Experts, This will be responded to on the Bioconductor mailing list; please address any follow-ups there. http://bioconductor.org/help/mailing-list/ Martin thank you for providing such a useful tool-set. I have a question reg

Re: [R] Bioconductor / AnnotationDbi: Why does a GOAllFrame contain more rows than its argument GoFrame?

2013-10-02 Thread Jeff Newmiller
You have addressed the wrong mailing list. See http://www.bioconductor.org/help/mailing-list/ --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go...

[R] ggplot2 multiple legends modification

2013-10-02 Thread rolf . kemper
Dear ALL, please find attached a small test file. (code is copied in this mail too) Basically (as shown) I would like to combine several results in one graph. I could succeed to plot the data correctly, but the legend generation does not fit. For the dfStartPoints and the dfEndpoints I just wan

[R] Drawing garbled

2013-10-02 Thread yy
Hi: I am Chinese, I am developing a java application, and deploy it to tomcat 7.0.42, I use rJava to use the R, when I use the command line to start tomcat, the R drawing well done, see attachment histview.png,but when i use windows service to start tomcat, the R drawing bad, see attach

Re: [R] Group all the consecutive days

2013-10-02 Thread arun
Hi, Try: dat1<- read.table(text='"year" "month" "day" "rainfall" "landslide" "3" 2007 6 6 1.6 0 "4" 2007 6 7 1.8 0 "6" 2007 6 12 4.6 0 "8" 2007 7 5 6.6 0 "9" 2007 7 10 3 0 "10" 2007 7 11 1.2 0 "11" 2007 8 3 6.4 0 "12" 2007 8 10 2.8 0 "14" 2007 9 4 5.4 0 "15" 2007 9 5 1 0 "16" 2007 9 10 2.8 0 "17"

[R] Interpreting the result of a Wilcoxon (Mann-Whitney U) test

2013-10-02 Thread Filipe Correia
Hello everyone, I'm having some trouble interpreting the results of a Wilcoxon (Mann-Whitney U) test. Hope you can help. This is the R script that I am running: a <- c(1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 5, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 1, 1) b <- c(1, 2, 1, 1, 2, 3, 2, 2, 1, 2, 1, 1, 1

[R] R Excel - Microsoft Excel is waiting for OLE action

2013-10-02 Thread Johnson, Alex
Hi, I am hoping you may be able to help please? I am trying to use R Excel to run an R script using "RInterface.RunRFile" I have tested this command on basic scripts and they work fine. However when I try to run quite a long script which runs a number of forecasts using the forecast package,

[R] (no subject)

2013-10-02 Thread Joe Byers
I am trying to get rJava to load the htmlunit java package (http://htmlunit.sourceforge.net/). rJava does not seem to be able to load and find the class in this package. This issue is the same on Windows 7 and Fedora 17 computers. rJava will load and init because I can see java.util.* classe

[R] Help with vapply() loop

2013-10-02 Thread Axel Urbiz
Hello, Although the example below doesn't necessary make any sense from a statistical perspective, it is just a close enough example to hopefully get your help upon. For my purpose, I'm particularly interested to know if there is a way to replace the results from the vapply() function below by a

[R] Bioconductor / AnnotationDbi: Why does a GOAllFrame contain more rows than its argument GoFrame?

2013-10-02 Thread Asis Hallab
Dear Bioconductor Experts, thank you for providing such a useful tool-set. I have a question regarding the package AnnotationDbi, specifically the classes GOFrame and GOALLFrame. During a GO Enrichment Analysis I create a data frame with Arabidopsis thaliana GO annotations and from that first a

Re: [R] missclassification in estimating proportions

2013-10-02 Thread David Winsemius
On Oct 2, 2013, at 12:15 AM, pattebegga wrote: > Hi > > Assume that I want to compare two methods, say skinfold measurement and BMI, > in how they classify subjects into four categories. In a 2x2 table I would > simply calculate sensitivity and specificity and NPV an PPV, using the > standard fo

Re: [R] Group all the consecutive days

2013-10-02 Thread jim holtman
try this: > rain <- read.table(text = '"year" "month" "day" "rainfall" "landslide" + "3" 2007 6 6 1.6 0 + "4" 2007 6 7 1.8 0 + "6" 2007 6 12 4.6 0 + "8" 2007 7 5 6.6 0 + "9" 2007 7 10 3 0 + "10" 2007 7 11 1.2 0 + "11" 2007 8 3 6.4 0 + "12" 2007 8 10 2.8 0 + "14" 2007 9 4 5.4 0 + "15" 2007 9 5 1 0

[R] update.packages fails on pdf write - CSAgent

2013-10-02 Thread Tim Howard
All, I work in a building where our computers are running Windows XP with Cisco Security Agent (CSA) activated, locking many automated actions, especially file writes. We (lowly staff) are allowed to install software, so on upgrading to R 3.0.2, I tried the standard approach for updating my package

Re: [R] leveneTest - output

2013-10-02 Thread S Ellison
> I have a question concerning the output of leveneTest. I don't > understand > the "7" in my output It's the residual degrees of freedom for the corresponding ANOVA table. You'd get the same value if you did a one-way ANOVA by group.ID. Having said that, your test isn't going to tell you much

Re: [R] Basic help on DF creation row by row

2013-10-02 Thread Rui Barradas
Hello, Inline. Em 02-10-2013 15:16, rolf.kem...@renesas.com escreveu: Hello Rui, thanks a lot for the quick and exhaustive reply ! It shows that my understanding of DF is totally wrong. I'm familiar with relational DBs where the columns of table have specific data types. Hence , I thought a

[R] How to create an ROC with three possible classes?

2013-10-02 Thread Mary Kindall
I have a toy example reproduced below in which the response variable has three possible classes. I am trying to create an ROC but not sure how to deal with it when there are three classes. library(ipred) control = rpart.control(maxdepth = 20, minsplit = 20, cp = 0.01, maxsurrogate=2, surro

Re: [R] How to pass an object between two independant environment?

2013-10-02 Thread Duncan Murdoch
On 02/10/2013 10:12 AM, Frederic Fournier wrote: Hello everyone, I am trying to pass an object cleanly between two independant environments, any ideas on how to do that? I have something like this: calling.function <- function(my.object) { ... ... local(source("function_to_which_I_want_to_pass

[R] How to pass an object between two independant environment?

2013-10-02 Thread Frederic Fournier
Hello everyone, I am trying to pass an object cleanly between two independant environments, any ideas on how to do that? I have something like this: calling.function <- function(my.object) { ... ... local(source("function_to_which_I_want_to_pass_my_object.Rda", local = new.env(parent = .Glob

[R] Group all the consecutive days

2013-10-02 Thread Stefano Sofia
Dear R-users, I have a data frame where in each row there is the daily rainfall cumulative; missing days mean that in that days rainfall has been zero. I need to group all the consecutive days in a single row and store in the field "rainfall" the sum of these consecutive days. Is there a reasonabl

Re: [R] optim evils

2013-10-02 Thread Duncan Murdoch
On 02/10/2013 9:18 AM, Michael Meyer wrote: Greetings, In obedient deference to the demands of the collective I emailed a BUG report containing code and data to r-b...@r-project.org but found subsequently that I am unable to load the page http://bugs.r-project.org/ to check on the status of th

Re: [R] transferring a graph from R into word

2013-10-02 Thread Flavio Barros
If you will use the graph on latex i recommend save as pdf. Ex: pdf(file='graph.pdf') plot( ... ) dev.off() and your graph will be saved on the working directory. For Microsoft Word i prefer png. You can set the resolution to the highest possible in order to get a resizable figure. Flavio Barr

[R] In which order should the "p" argument within chisq.test() be specified when working in 3 dimensions?

2013-10-02 Thread Peter Davidsen
Dear List, I'm in the process of setting up a chi square test with 3 groups (A, B and C). Hence, I have generated an array object like this: my.obj <- array(c(overlap_ABC, overlap_AB, overlap_AC, unique_A, overlap_BC, unique_B, unique_C, not_in_any), dim=c(2,2,2)) Now I'm in the process of gener

[R] r - Forecast model

2013-10-02 Thread Stephen Donnelly
Hi All, I'm currently completing my MSc dissertation, using R to build a realised volatility model using HAR. There is a great guide on using high frequency in R which has been invaluable to me. I was wondering though if you could all help with one issue? I've successfully achieved page 16 exam

Re: [R] package:nlme unexpected varIdent behaviour

2013-10-02 Thread Kevin Wright
See this post: http://r.789695.n4.nabble.com/varIdent-error-using-gam-function-in-mgcv-td2260556.html Kevin On Tue, Oct 1, 2013 at 4:37 PM, Sonya Odsen wrote: > I'm running R version 3.0.2 (2013-09-25) on a Windows computer. > Here is a simplified version of my data: > > a <- expand.grid(SITE

Re: [R] Vector from csv file.

2013-10-02 Thread arun
Hi, It is not clear whether you read the file correctly or not.  Please check str(tbl_all) The code you showed should work if the file was read correctly. For example: set.seed(468) tbl_all<- as.data.frame(matrix(sample(1:10,10*22,replace=TRUE),ncol=22,dimnames=list(NULL,c("A","V",LETTERS[3

[R] optim evils

2013-10-02 Thread Michael Meyer
Greetings, In obedient deference to the demands of the collective I emailed a BUG report containing code and data to r-b...@r-project.org but found subsequently that I am unable to load the page http://bugs.r-project.org/ to check on the status of this report. Can anyone else load this page?

Re: [R] Removing time and min from Date

2013-10-02 Thread arun
Hi Farnoosh, Use ?as.Date() dat1<- read.table(text="V1 2012-01-12 08:23:00 2012-01-19 15:29:00",sep=",",header=TRUE,stringsAsFactors=FALSE) dat1$V1<- as.Date(dat1$V1)  dat1 #  V1 #1 2012-01-12 #2 2012-01-19 A.K. From: farnoosh sheikhi To: "smartpink

Re: [R] missclassification in estimating proportions

2013-10-02 Thread Frank Harrell
This is dichotomania and is an inappropriate use of continuous variables. Use an information-preserving approach such as rank correlation. Also, this is not an R question. Frank --- Assume that I want to compare two methods, say skinfold measurement and BMI, in how they classify subj

Re: [R] Command line r

2013-10-02 Thread peter dalgaard
On Oct 2, 2013, at 11:07 , C.H. wrote: > Have a look at the littler. > > http://dirk.eddelbuettel.com/code/littler.html > Or Rscript (comes standard with R) for that matter: $ echo '2+2' | Rscript - [1] 4 -pd > On Sun, Sep 29, 2013 at 3:36 AM, Tomek R wrote: >> Hi, >> I have found myself

Re: [R] over dispersion plot

2013-10-02 Thread Michael Dewey
At 19:15 01/10/2013, Sagan Friant wrote: I want to show over dispersion in my data by plotting the variance against the mean for 8 groups. I used the following code to plot by type (parasite). scatterplot(variance~mean | parasite, reg.line=lm, smooth=FALSE, spread=FALSE, boxplots=FALSE, span=0.

Re: [R] Command line r

2013-10-02 Thread C.H.
Have a look at the littler. http://dirk.eddelbuettel.com/code/littler.html On Sun, Sep 29, 2013 at 3:36 AM, Tomek R wrote: > Hi, > I have found myself often doing simple statistical analysis using Linux > command line on a single dataset. Therefore, I put a perl script together, > which makes

Re: [R] transferring a graph from R into word

2013-10-02 Thread Rui Barradas
Hello, Left click on the graph and then on the file menu, choose "save as". You will have a choice of graphics file formats. Word accepts several of them, so it's up for you to choose. PNG or BMP, for instance. Hope this helps, Rui Barradas Em 02-10-2013 07:32, Charles Thuo escreveu: I hav

Re: [R] transferring a graph from R into word

2013-10-02 Thread PIKAL Petr
Hi I am not sure if right clicking was meant as standard way to transferring graph to Word. Anyway there shall be menu with several options for saving graph or copying them to clipboard. You probably shall also check ?pdf, ?png, ?postsript or other devices. Regards Petr > -Original Mes

[R] missclassification in estimating proportions

2013-10-02 Thread pattebegga
Hi Assume that I want to compare two methods, say skinfold measurement and BMI, in how they classify subjects into four categories. In a 2x2 table I would simply calculate sensitivity and specificity and NPV an PPV, using the standard formulas, but can the same be directly applied in the 4x4 table