similar to: power of a multiway ANOVA

Hi, does anyone know of a function to calculate odds ratios in multiway tables (stratified) (+ the other usual statistics involved) i mean: say we have a table r*c*d, For every d (depth) we have a r*c table, and in this table the odds ratio's are calculated for every 2*2 subtable in it. logically this function would look like): ORs(multiwaytable) or ORs(data$var1r,data$var2c,data$var3d)

R has amazing capabilities, but percentage tables are a weak spot IMHO. There's prop.table but that's rather unwieldly, especially for multiway tables. CrossTable by Marc Schwartz in the gregmisc library makes percentage tables a breeze but is limited to two-way tables. So I decided to try my own hand at writing an R-function that would make it easy to produce nicely formatted percentage

Full_Name: Michael J. Lew Version: R version 2.9.1 (2009-06-26) OS: OS X Submission from: (NULL) (210.49.195.149) The default function power.t.test gives inaccurate values, particularly with large sig.level. There are two reasons for me to propose that power.t.test is inaccurate: 1. When sig.level approaches 1, so should the power. However,

Dear R experts, I have conducted a power calculation in order to estimate the number of subjects needed to detect an effect size of d=0.28 (cohen's d) for a difference between two independent groups (alpha level should be 0.05 and the effect should be detected with 80% probability). The results from the code below indicates that I would need n=400 subjects (200 in each group). This is seems

Craig, Thanks for your follow-up note on using the asypow package. My problem was not only constructing the "constraints" vector but, for my particular situation (Poisson regression, two groups, sample sizes of (1081,3180), I get very different results using asypow package compared to my other (home grown) approaches. library(asypow) pois.mean<-c(0.0065,0.0003) info.pois <-

PTAk_1.1-1 ( Principal Tensor Analysis on k modes) has been released on CRAN A multiway method to decompose a tensor (array) of any order, as a generalisation of SVD also supporting non-identity metrics and penalisations. 2-way SVD with these extensions is also available. The package includes also some other multiway methods: PCAn (Tucker-n) and

PTAk_1.1-1 ( Principal Tensor Analysis on k modes) has been released on CRAN A multiway method to decompose a tensor (array) of any order, as a generalisation of SVD also supporting non-identity metrics and penalisations. 2-way SVD with these extensions is also available. The package includes also some other multiway methods: PCAn (Tucker-n) and

Hello R-users, I'm currently facing a pretty hard problem which I'm hopping you'll be able to help me with. I'm using R to create images. That alone is not the problem, the problem is that I'm using R to create 168 000 images... My code (which is given below) use different package (raster and rgdal) to import a image (size 20gig) and divide it into 168 000 pictures that are

hello, When I compil the samba new version samba 3.6.5, on Debian jessie 8.8, most of the time I've got an error message. I can compile without problems samba 3.6.4 on the running servers. It doesn't happen on a fresh install Thanks for your help, - serge. 'configure' finished successfully (1m56.906s) WAF_MAKE=1 python ./buildtools/bin/waf build Waf: Entering directory

Hello, I use a type III anova ("car" package) to analyse an unbalanced data design. I have two factors and I would have the effect of the interaction. I read that the result could be strongly influenced by the contrasts. I am really not an expert and I am not sure to understand indeed about what it is... Consequently, I failed to properly used the fit.contrast function (gregmisc

I have an unbalanced design I would like to run a power analysis on. What I have been able to find has pointed me to using the pwr.f2.test function as described below. My problem is that I don't know how to appropriately define the numerator and demoninator df. If someone can help here is some more info about my design. It is an unbalanced 2^3 x 3 design where the factor with 3 levels is a

Bonjour, Je veux cherché l’estimateur de vraisemblance maximal (MLE)d’une fonction à 3 paramètre inconue étant donné une échantillon de taille 50 (les observations des valeurs de x) alors comment je peux procédé La fonction de densité est définie par : f(x)= 1/3(g(a1)+g(a2)+g(a3)) avec g(ai)=(exp(ai)*ai^x)/x! pour i=1,2,3. Je vous remercie beaucoup. A. Elhabti

Dear all, my colleague (a statistician - not yet using R) aked me, if there is a R package comparable to the software "PASS" (Power Analysis and Sample Size Procedures), comprising about 150 procedures. I found the R package "pwr". Are there other (more comprehensive) power analysis packages? Many thanks! Udo

Hi, I am writing a decompiler. I was wondering if some of LLVM could be used for a decompiler. There are several stages in the decompiler process. 1) Take binary and create a higher level representation of it. Like RTL. 2) The output is then broken into blocks or nodes, each block ends in a CALL, JMP, RET, or 2-way or multiway conditional JMP. 3) The blocks or nodes are then analyzed for

Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Strangely, I never get access to the pwr object in R. I tried installing it to /usr/local/lib/R/site-library and ~/R/x86_64-pc-linux-gnu-library. I also tried installing it

Hi all, I would like to estimate power and necessary sample size for a GLM with a response variable that has a poisson distribution. Do you have any suggestions for how I can do this in R? Thank you for your help. Sincerely, Craig -- Craig A. Faulhaber Department of Forest, Range, and Wildlife Sciences Utah State University 5230 Old Main Hill Logan, UT 84322 (435)797-3892

Dear Any, Is there a fonction in R to change a string to uppercase ? Thanks for all your help

This fixes: btrfs-list.c: Dans la fonction «ino_resolve» : btrfs-list.c:511: attention : déclaration implicite de la fonction « «strndup» » btrfs-list.c:511: attention : incompatible implicit declaration of built-in function «strndup» make: *** [btrfs-list.o] Erreur 1 and: btrfs.c: Dans la fonction «split_command» : btrfs.c:168: attention : déclaration implicite de la fonction « «strndup» »

Hello, Is there any R function which does power calculation for unbalanced groups (n1 neq n2)? Since power.t.test has n Number of observations (per group). Many thanks, Samuel [[alternative HTML version deleted]]

My apologies for the statistical naivete of my question but... Is there an established method or calulating the statistical power of a correlation test? And if so is there a method in R for it? Thank you, Collin Lynch.