search for: mypow

Hi, I am trying to create a power curve to show how the power of a t-test varies depending on the mean. Any ideas how I should go about this? Louisa _________________________________________________________________ [[alternative HTML version deleted]]

I am trying to store regression objects in a data.table df <- data.frame(x = rnorm(20)) df[, "y"] <- with(df, x + 0.1 * x^2 + 0.2 * rnorm(20)) mydt <- data.table(mypower = c(1, 2), myreg = list(lm(y ~ x, data = df), lm(y ~ x + I(x^2), data = df))) mydt #?? mypower??? myreg #???? <num>?? <list> #1:?????? 1 <lm[12]> #2:?????? 2 <lm[12]> But mydt[1, 2] has only the coeffients of the first regression. mydt[2, 2] has residuals of the first re...

...21, 2024 at 2:25?PM Naresh Gurbuxani <naresh_gurbuxani at hotmail.com> wrote: > > I am trying to store regression objects in a data.table > > df <- data.frame(x = rnorm(20)) > df[, "y"] <- with(df, x + 0.1 * x^2 + 0.2 * rnorm(20)) > > mydt <- data.table(mypower = c(1, 2), myreg = list(lm(y ~ x, data = df), > lm(y ~ x + I(x^2), data = df))) > > mydt > # mypower myreg > # <num> <list> > #1: 1 <lm[12]> > #2: 2 <lm[12]> > > But mydt[1, 2] has only the coeffients of the first regressio...

...parameter for negative binomial, infinity correspond to poisson # #example power analysis #return three power, power1 for zero part, power2 for negative binomial part #power3 for joint test,significance level can be set, default is 0.05 #M is simulation time #require pscl package #library(pscl) mypower <- function(n, beta0, beta1, alpha0, alpha1, theta, siglevel=0.05, M=1000){ myfun <- function(n,beta0,beta1,alpha0,alpha1,theta,siglevel){ data <- mysim(n,beta0,beta1,alpha0,alpha1,theta) require(pscl, quietly = TRUE) res <- hurdle(y ~ trt, data = data, dist = "negbin&quo...

...sh Gurbuxani > <naresh_gurbuxani at hotmail.com> wrote: >> I am trying to store regression objects in a data.table >> >> df <- data.frame(x = rnorm(20)) >> df[, "y"] <- with(df, x + 0.1 * x^2 + 0.2 * rnorm(20)) >> >> mydt <- data.table(mypower = c(1, 2), myreg = list(lm(y ~ x, data = df), >> lm(y ~ x + I(x^2), data = df))) >> >> mydt >> # mypower myreg >> # <num> <list> >> #1: 1 <lm[12]> >> #2: 2 <lm[12]> >> >> But mydt[1, 2] has only t...

Dear all, I have problems to compile a DLL for Windows XP with R version 2.5.1 (2007-06-27). See output below: C:\AZ_DATEN\C, C++\BE_speedup>R CMD SHLIB dmypow.c Goto undefined subroutine &DynaLoader::bootstrap_inherit at C:\Programme\R\R-2.5 .1\share\perl/XSLoader.pm line 80. Compilation failed in require at c:/Programme/Perl/lib/IO/Handle.pm line 262. BEGIN failed--compilation aborted at c:/Programme/Perl/lib/IO/Handle.pm line 262 . Compilation fail...

> R code works but the shared objects built from fortran code do not find > the appropriate libs with functions like __pow_ii or __epx at runtime. I used to have endless trouble with __pow_ii which as I recall could be be fixed by the correct incantation of the mix of compiler libraries and LD_LIBRARY_PATH. It seemed to be extremely sensitive and never easy to get right. The problem is