similar to: simple Q: required sample size & usage of power.t.test()

Dear Reader, I need a horizontal barplot but with the lowest value on the y-axis on top. How to get this? Thank you, osiander at 24on.cc -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the

Hi everyone, I'm working with the Epical library, specicatly using the power test in proportions. I think this test is not working like in the book: Biostatistical Analysis (4th Edition): Jerrold H. Zar In the example 23.25. (I attach this Pic) It's not the same answer. Using the follow command don't give the same answer. library(epicalc) power.for.2p(0.75, 0.50, 50, 45, alpha =

Hello R-help, I have a data.frame where I want to scatterplot variable1 versus variable2 using a third variable as a factor. For each factor I want to have a regression line in the same plot with different color and symbols. Var1, and var2 contain a lot of NA,sometimes the dependent var sometimes the predictor, sometimes both together. Three problems occure: Makeing subsets to split the data

Dear R users: Suppose I have two different response variables y1, y2 that I regress separately on the different explanatory variables, x1 and x2 respectively. I need to compare points on two regression lines. These are the x and y values for each lines. x1<-c(0.5,1.0,2.5,5.0,10.0) y1<-c(204,407,1195,27404313) x2<-c(2.5,5.0,10.0,25.0) y2<-c(440,713,1520,2634) Suppose we need to

Full_Name: Osiander Meixner Version: 2.9.0 (2009-04-17) OS: Windows Vista Enterprise, 32bit, SP1 Submission from: (NULL) (15.195.185.82) package XML_2.5-1 >library(XML) throws an error window saying: " RGui: RGui.exe - Unable to Locate Component This application has failed to start becuase iconv.dll was not found. Re-installing the application may fix this problem. " RGui then

Please install from the "CRAN extras" repository (where you will find a working copy) and do NOT report a bug in R that is not a bug in R but in a contributed package / the package repository. Best, Uwe Ligges osiander.meixner at hp.com wrote: > Full_Name: Osiander Meixner > Version: 2.9.0 (2009-04-17) > OS: Windows Vista Enterprise, 32bit, SP1 > Submission from: (NULL)

I am using nginx with unicorn as a reverse proxy. Has been fantastic, but needed to look through some code for first time. Using ruby 1.9.3p194 with rails 3.2.9, development mode with ssl enabled, I start unicorn with -d ...and it does not get served up by nginx. Any ideas, or more information needed? Thanks! Jet Below is my unicorn.rb for developent: # config/unicorn.rb env =

Hi, I'm having difficulty to replicate in R a nested (hierarchical) anova example found in p. 308 of Zar, J.H. 1996. Bostatistical Analysis. Prentice Hall. 3rd ed. The example (15.1) is as follows: The variable is blood cholesterol concentration in women (in mg/100 ml of plasma). This variable was measured after the administration of one of three different drugs, each drug having been

Hi everyone, I have a data frame (d), wich has the results of mosquitoes trapping in three different places. I suspect that one of these places (Local=='Palm') is biased by low numbers and will yield slower slopes in the variance-mean regression over the areas. I wonder if these slopes are diferents. I've looked trought the support list for methods for comparing slopes and found the

[This was filtered (looking like spam) and manually approved by your list maintainer, MM ] EvolDir wrote: > Colleagues, > > A few days ago I posted a message inquiring what folks were using for > teaching biostats. I received a stunning 38 replies offering various > opinions, some quite surprising, along with links to web pages, pointers for > coverage, commentaries on

I am trying to do a Nonparametric Tukey-type multiple comparison post-hoc test to determine which groups are significantly different. I have read the dialogue on this topic from the R-help, and am still not clear why no statistical packages include this test as an option? Is it not an appropriate test to conduct on non-normally distributed data? Is the only option to calculate it by hand

Hello all, Here goes one of my first functions. I want to make a nonparametric multiple sample comparison with unequal sample sizes (see Zar?s Biostatistical Analysis, 3rd. Ed., pg. 201 Example 10.11, pg. 288 Example 11.10). In the real world, I want to compare samples of fish length captured with different fishing gears. After using the Kruskal-Wallis test I want to check the differences

I've written a g.test() aka log-likelihood ratio test function for my opwn use. It's something I've seen requested (and looked to find myself) on this list a few times. It has the same basic syntax as chisq.test(). It does both goodness of fit tests and tests of independence. Yates' and Williams' corrections are implemented. I've put some examples from Sokal & Rohlf

Error file begins: Unhandled exception: unimplemented function msvcp80.dll.??0?$basic_ifstream at DU?$char_traits at D@std@@@std@@QAE at PAU_iobuf@@@Z called in 32-bit code (0x7b83aa52). I downloaded a new copy of this and placed it in the /fakedlls dirctory, where it replaced what was there, but it still doesn't work. The site I d/l'd the file from gave the impression that there is

Hi, I'm newbie in R package, but i'm extremely gratefull for these list! is there possible to make kruskal.test summarize the sum of ranks?? e.g ##### z<-c(1.2,3.4,0,0,0,0,0) ranks(z) # like in octave package ans= 6 7 3 3 3 3 3 sum(ranks(z)) ans= 28 ### then apply multiple comparisions for many other data SE<-sqrt(n*(n*k)*(n*k+1)/12) # according Zar,JH. Bioestatistical

I need to specify error terms properly in a mixed-effects anova model. I know you can add error terms in aov using Error and can specify random factors in lme but I am not sure how these get treated. When making the calculations for fixed and random factors, are the correct error terms used and how can you get aov or lme to use different error terms for fixed and random effects? I'm

How I make the below matrix | a^0 0 0 ... 0 | | a^1 a^0 0 ... 0 | | a^2 a^1 a^0 ... 0 | | . | | . | | . | | a^n a^(n-1) a^(n-2) ... a^0 | with no loops, where "a" is a constant? ======================================================================== C?zar de Freitas Depto. de Estat?stica - UFPE Recife - Pernambuco - Brasil

I'm testing out quantstrat using a simple one-indicator strategy. The error I get after applyStrategy(...) is Error in .xts(e, .index(e1), .indexCLASS = indexClass(e1), .indexFORMAT = indexFormat(e1), : index length must match number of observations In addition: Warning messages: 1: In match.names(column, colnames(data)) : all columns not located in roc_15 for STOXX.Open STOXX.High

Hi. I searched the list and didn't found nothing similar to this. I simplified my example like below: #I need calculate correlation (for example) between 2 columns classified by a third one at a data.frame, like below: #number of rows nr = 10 #the third column is to enforce that I need correlation on two variables only dataf =

is there a way to use the actual index value for plotting zoo objects this is the way that the index is set up and a sample range of what I would like to plot 01/01/06 00:00:00 - 01/01/06 23:45:00 { library(zoo) # chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) }} x <- structure(c(15.57, 15.5,