similar to: function for the average or expected range?; CORECTION

Hi, All: ** My previous email on this subject seemed to contain an error; check the correction below; please excuse my errors. Is there a function in R to compute the expected range of a sample of size n from some distribution? I ask, because I was recently asked about the control chart constant 'd2', which is the expected range for a sample of size n from a standard normal.

Hi, All: Is there a function in R to compute the expected range of a sample of size n from some distribution? I ask, because I was recently asked about the control chart constant 'd2', which is the expected range for a sample of size n from a standard normal. There is a fairly simple formula for the expected value of the range, given, e.g., in Kendall and Stuart (1969)

I am embarrassed to have to ask this but can anyone tell me of a Tukey multiple comparisons procedure available for R? I have looked through the search page, through the FAQ, and in the index of V&R (1999), and I still can't find such a thing. I see there is a ptukey function and a qtukey function but that is as far as I got. Do I need to roll my own? -- Douglas Bates

Dear all, I'm using R to perform multiple comparison testing on agriculture genotype trials. To perform the Duncan's MRT, I use the qtukey() function with the following syntax: qtukey(p = ((1 - 0.05) ^ (pos - 1)), nmeans = pos, df = ni) I experience a strange behaviour when the number of means in the trial and the number of residual degrees of freedom (ni) becomes high

Well, it is clear enough that the problem is in interpreting the documentation. However, when you claim you tested something, and found it inconsistent with tables, it would be advisable to back it up with examples! The description in the help files and in the sources is admittedly confusing. The original paper has this, rather more clear, description in the abstract: "We consider the

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after

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,

Dear all, I wanted to compare Bonferroni vs TukeyHSD correction over a range of groups and group sizes, and wanted to use the function qtukey. In the help documentation it says qtukey(p, nmeans, df, nranges = 1, lower.tail = TRUE, log.p = FALSE) Arguments q vector of quantiles. p vector of probabilities. nmeans sample size for range (same for each group). df degrees of freedom for s (see

We cannot help you understand what you are doing if you do not show us what you are doing. Here are some discussions about how to communicate questions about R [1][2][3]. [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example [2] http://adv-r.had.co.nz/Reproducibility.html [3] https://cran.r-project.org/web/packages/reprex/index.html -- Sent from my phone.

On a recent FreeBSD 8.0-CURRENT (i386) building R (any version) breaks with the following messages: ---------------------------------------------------------------------- [...snip...] gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c wilcox.c -o wilcox.o gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include

Hi, I'm building Tukey's table using qtukey function. It happens that I can't get the values of Tukey's one degree of freedom and also wanted to eliminate the first column. The program is: Trat <- c(1:30) # number of treatments gl <- c(1:30, 40, 60, 120) # degree freedom tukval <- matrix(0, nr=length(gl), nc=length(Trat)) for(i in

-- View this message in context: http://r.789695.n4.nabble.com/Generalized-DCC-GARCH-ML-estimation-tp2245125p2245125.html Sent from the R help mailing list archive at Nabble.com.

I was wondering whether there exists a naming convention for row and column names in R data frames and matrices. E.g: PriceSwap or PRICESWAP or PRICE.SWAP Many thanks. Reinhold [[alternative HTML version deleted]]

Please note that building R on Windows is documented in "R Installation and Administration", including links to external software. Particularly there is a link to texinfo which is part of Rtools. The documentation is maintained and it is a sufficient source of information for building R on Windows. https://cran.r-project.org/doc/manuals/r-release/R-admin.html

Hi, I am having trouble using the ave function with a POSIXct object. For example: x<-Sys.time()+0:9*3600 dat<-data.frame(id=rep(c('a',' b','c'),each=10),dt=rep(x,3),i=rep(1:10,3)) dat # This is what I want to do: dat$time.elapsed<-unsplit(lapply(split(dat$dt,dat$id),function(x) x-x[1]),f=dat$id) dat # The above code does the trick, but from the standpoint of

Dear all, I'm writing my manuscript to publish after analysis my final data with ANOVA, ANCOVA, MANCOVA. In a section of my result, I did correlation of my data (2 categirical factors with 2 levels: Quantity & Quality; 2 dependent var: Irid.area & Casa.PC1, and 1 co-var: SL). But as some traits (here Irid.area) are significantly influenced by the covariate (standard length, SL), I

Dear Eliza, Try this: Lines1<-readLines(textConnection("1911.01.01?????? 7.87 1911.01.02?????? 9.26 1911.01.03?????? 8.06 1911.01.04?????? 8.13 1911.01.05????? 12.90 1911.02.06?????? 5.45 1911.02.07?????? 3.26 1911.03.08?????? 5.70 1911.03.09?????? 9.24 1911.04.10?????? 7.60 1911.05.11????? 14.82 1911.05.12????? 14.10 1911.06.13?????? 7.87 1911.06.14?????? 9.26

# R for Windows will not send your bug report automatically. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs at r-project.org # ###################################################### This report is joint from Richard Heiberger <rmh at temple.edu> and Burt Holland <bholland at temple.edu>. Burt Holland is the coauthor

I am not sure if my message made it through, so here it is again! Hi Rer's, R-1.1.0 I have two questions for you: 1) I am trying to complete a multiple comparison test after completing a one-way ANOVA on some data. I think this is pretty reasonable. aov(MetricSubset ~ GeneNameFactor) works pairwise.t.test(MetricSubset,GeneNameSubset,p.adjuxt.method=bonferroni,p ool.sd=FALSE)

I am interested in a statistical comparison of multiple (5) time series' generated from modeling software (Hydrologic Simulation Program Fortran). The model output simulates daily bacteria concentration in a stream. The multiple time series' are a result of varying our representation of the stream within the model. Our main question is: Do the different methods used to represent a