Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread BONACHE Adrien via R-help
Hi again, In fact the pwr package do not provide exact test power. So you should definitely use g*power. Sincerely, Adrien. De : Munjal Patel À : r-help@r-project.org Envoyé le : Mardi 7 juin 2016 18h26 Objet : [R] Power Calculation:2-sided exact equivalence test for Binomial Proporti

Re: [R] Power Calculation: Binomial Proportions (2 sided exact test for equivalence)

2016-06-07 Thread BONACHE Adrien via R-help
Dear Munjal, You should use the pwr package with the following code :library(pwr) # lower pwr.2p.test(h =ES.h(0.3, 0.2), n = 500, sig.level = 0.05, power = NULL, alternative = c("two.sided"))#upperpwr.2p.test(h =ES.h(0.3, 0.4), n = 500, sig.level = 0.05, power = NULL, alternative = c("two.sided")

Re: [R] Reading and converting time data via read.table

2016-06-07 Thread Ek Esawi
Thanks for the responses i received from David and Spencer. As for the package Ecfun, i looked at it briefly, but it's long and i am new to R; so this just adds to my confusion. David suggested the sub function which i will try and see what i get.i s As i said earlier i have no problems reading da

Re: [R] Faster Multivariate Normal

2016-06-07 Thread Jeff Newmiller
The help file ?dmvnorm is your friend. Read about "x". ghv <- matrix( gh[ as.vector( idx ) ], ncol = dm ) adjFactor2 <- dmvnorm( ghv, mean = mu, sigma = sigma ) -- Sent from my phone. Please excuse my brevity. On June 7, 2016 10:19:53 AM PDT, "Doran, Harold" wrote: >Thanks, Duncan. Not sure I f

Re: [R] Reading and converting time data via read.table

2016-06-07 Thread Spencer Graves
Have you considered asNumericDF in the Ecfun package? This can convert to a date-time column into POSIXct -- not POSIXlt -- using the POSIX argument to identify the POSIX columns and format to provide the "format" argument for as.POSIXct. I suggest you try the R-Forge version: in

Re: [R] Cut Dates into bins

2016-06-07 Thread TJUN KIAT TEO
It does not seem to work for me. I will show you exactly my data format SMA$TIME_DATE "28/9/2014" "17/6/2014" "19/9/2014" "17/1/2015" "13/1/2015" "21/10/2014" Shipment$DateRequire "2014-06-09" "2014-06-16" "2014-06-16" "2014-06-16" "2014-06-17" "2014-06-23" What I would like to do is the

Re: [R] Reading and converting time data via read.table

2016-06-07 Thread Ek Esawi
Thanks Petr! I am still unable to come up with a conversion formula/trick to convert my time data to POSIXlt which then can be used in read.table. Below are again a few lines from my file. Since as you see there are several columns of time data ONLY (no date), I am hoping that there is a way to

[R] Revolutions blog: May 2016 roundup

2016-06-07 Thread David Smith via R-help
Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. And in case you

Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread Leonardo Ferreira Fontenelle
Em Ter 7 jun. 2016, às 13:26, Munjal Patel escreveu: > Dear R-Users, > I am an intermediate level R user. > > I am performing the power calculations for the Binomial proportions (2 > sided). > I want to find the Power using the Exact test for the Equivalence of > Binomial proportions. > I do have

Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread Bert Gunter
Please search before posting, and if your search fails to get what you want, tell us why. I got what appeared to be many relevant hits on rseek.org using the search term "binomial exact power computations" . -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming

Re: [R] Power Calculation: Binomial Proportions (2 sided exact test for equivalence)

2016-06-07 Thread Leonardo Ferreira Fontenelle
Em Ter 7 jun. 2016, às 13:08, Munjal Patel escreveu: > Dear R-Sig-teaching users, > I am an intermediate level R user. You posted both emails to the same mailing list. Please remember that "cross-posting is considered to be impolite" and that "you should configure your e-mail software in such a w

Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread Cade, Brian
MJ: I think the EnvStats package has various power functions for binomial applications (also confidence interval half-widths). Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 22

[R] Power Calculation: Binomial Proportions (2 sided exact test for equivalence)

2016-06-07 Thread Munjal Patel
Dear R-Sig-teaching users, I am an intermediate level R user. I am performing the power calculations for the Binomial proportions (2 sided). I want to find the Power using the Exact test for the Equivalence of Binomial proportions. I do have the SAS code which is generating the Power for me but i

[R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread Munjal Patel
Dear R-Users, I am an intermediate level R user. I am performing the power calculations for the Binomial proportions (2 sided). I want to find the Power using the Exact test for the Equivalence of Binomial proportions. I do have the SAS code which is generating the Power for me but i am unable to

Re: [R] Faster Multivariate Normal

2016-06-07 Thread Doran, Harold
Thanks, Duncan. Not sure I follow, however. The call to dmvnorm(), in this instance, takes in a 4 x 1 vector of nodes (in addition to the mean and covariances matrices), such as in the example below which uses the original sample code. That vector of nodes changes for each iteration of the loop

Re: [R] Faster Multivariate Normal

2016-06-07 Thread Duncan Murdoch
On 07/06/2016 9:06 AM, Doran, Harold wrote: I am computing a complex multidimensional integral (4 dimensions) using Gauss-legendre and am looking to optimize one piece of code that is currently very expensive. The integral is evaluated for K individuals in my data and once this piece is comput

Re: [R] evaluate the daily mean with date in "POSIXct" "POSIXt" class

2016-06-07 Thread Stefano Sofia
Thank you Jim, you gave me the right hint. Stefano Da: Jim Lemon [drjimle...@gmail.com] Inviato: martedì 7 giugno 2016 9.39 A: Stefano Sofia Cc: r-help@r-project.org Oggetto: Re: [R] evaluate the daily mean with date in "POSIXct" "POSIXt" class Hi Stefano,

[R] Faster Multivariate Normal

2016-06-07 Thread Doran, Harold
I am computing a complex multidimensional integral (4 dimensions) using Gauss-legendre and am looking to optimize one piece of code that is currently very expensive. The integral is evaluated for K individuals in my data and once this piece is computed it can be stored and recycled over all K in

Re: [R] Antwort: Re: Merging variables

2016-06-07 Thread PIKAL Petr
Hi > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of > g.maub...@weinwolf.de > Sent: Tuesday, June 7, 2016 8:19 AM > To: Michael Dewey > Cc: r-help@r-project.org > Subject: [R] Antwort: Re: Merging variables > > Hi Michael, > > yes, I was astonished ab

Re: [R] Antwort: RE: Merging variables

2016-06-07 Thread PIKAL Petr
Hi see in line > -Original Message- > From: g.maub...@weinwolf.de [mailto:g.maub...@weinwolf.de] > Sent: Tuesday, June 7, 2016 8:35 AM > To: PIKAL Petr > Cc: Michael Dewey ; r-help@r-project.org > Subject: Antwort: RE: [R] Merging variables > > Hi Petr, > > I would like to describe the d

Re: [R] evaluate the daily mean with date in "POSIXct" "POSIXt" class

2016-06-07 Thread Jim Lemon
Hi Stefano, I might be missing something, but try this: MteBove<-read.table(text="posix_date posix_time snowtemp 2010-01-19 23:30:00 45 NA 2010-01-20 00:30:00 10 2.7 2010-01-20 03:00:00 45 NA 2010-01-20 03:30:00 44 NA 2010-01-20 04:00:00 44 NA 2010-01-20 04:30:00 44 NA 2010-01-20 05:0