Re: [R] funnel plot asymmetry

2016-08-10 Thread BONACHE Adrien via R-help
Dear Christos, Maybe you should read it before using Egger's test : http://www.cienciasinseso.com/en/tag/eggers-test/If you still want to perform Egger's test, use metafor : R: Test for funnel plot asymmetry |   | |   |   |   |   |   | | R: Test for funnel plot asymmetrymetabias.meta {meta} R Do

Re: [R] Forking and adapting an R package

2016-07-14 Thread BONACHE Adrien via R-help
Hi Timo, To perform the first point, you just have to write the function name in R without using parentheses and arguments after it. You will see the code of the function. Copy and Paste it on notepad, change the name of the function in the notepad, then change what you want to change in your fu

Re: [R] t-test for regression estimate

2016-06-28 Thread BONACHE Adrien via R-help
Dear Steven, I understand your request now. Maybe you could try to change the function lm and create another function enabling you to have the information you want really quickly. Use lm() to see the code of this command and change this one marginally to get what you want to. Once this is done,

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")

[R] Is my simulation to compute power of a multiple ordinal logistic regression right?

2016-05-03 Thread BONACHE Adrien via R-help
Good day, I was performing a power analysis of articles published in a journal of management using the pwr package in R. However, it seemed to be impossible to compute power for small, medium and large Effect Size for multiple ordinal logistic regression output. I have tried using G*power, but i