similar to: Help with cochran.test

Hi there fellow R-users, Can anyone tell me how to build a prediction interval for a gaussian log-link model for the reponse variable?? I can find the standard error of the predictions but I cant seem to find the prediction interval. Is there a way I can calculate the prediction interval from the standard errors?? Here's the example: logX<-rnorm(100)

Hello I want to calculate natural logarithm of sum of combinations as follow: (R code) { com_sum=choose(2000000,482)*choose(1000000,118)+choose(2000000,483)*choose(1000000,117)+...+choose(2000000,i)*choose(1000000,600-i)+...+choose(2000000,600)*choose(1000000,0) #calculate the sum result=log(com_sum) #calculate the log of the sum } But

Hello I have a time-comsuming program which need to simplify, I have tested the annotated program as follow: > #define function which will be call > calsta <- function(c, n=100000) + { + i <- seq(from=0, length=c) + logx <- lchoose(NT-n, CT-i) + lchoose(n, i) + logmax <- max(logx) + logmax + log(sum(exp(logx - logmax))) + } > CT=6000 #assignment to CT >

Hi all I was playing with termplot(), and came across what appears to be an inconsistency. It would appreciate if someone could enlighten me: > # First, generate some data: > y <- rnorm(100) > x <- runif(length(y),1,2) > # Now find the log of x: > logx <- log(x) > > # Now fit two models that are exactly the same, but specified differently: > m1 <-

Dear R-users, I am using R version 2.10.1 under windows. In a barplot, I want to mark one of the bars with a special border color. For example: barplot(c(3, 7, 11), border = c(NA, "red", NA)) But how to do this when the bars are stacked? for example: barplot(matrix(1:6, ncol=3)) # border of second bar (i.e. the one with total height = 7) should be red again, I try: barplot(matrix(1:6,

Hello T I want to order some calculation "result", there will be lots of "result" that need to calculate and order PS: the "result" is just a intermediate varible and ordering them is the very aim # problem: # For fixed NT and CT, and some pair (c,n). order the pair by corresponding result # c and n are both random variable CT<-6000 #assignment to CT

>>>>> William Dunlap via R-devel >>>>> on Sun, 23 Jun 2019 10:34:47 -0700 writes: >>>>> William Dunlap via R-devel >>>>> on Sun, 23 Jun 2019 10:34:47 -0700 writes: > include/Rmath.h declares a set of 'logspace' functions for use at the C > level. I don't think there are core R functions that call

I agree with many the sentiments about the wisdom of computing very small p-values (although the example below may win some kind of a prize: I've seen people talking about p-values of the order of 10^(-2000), but never 10^(-(10^8)) !). That said, there are a several tricks for getting more reasonable sums of very small probabilities. The first is to scale the p-values by dividing the

>>>>> jing hua zhao >>>>> on Mon, 24 Jun 2019 08:51:43 +0000 writes: > Hi All, > Thanks for all your comments which allows me to appreciate more of these in Python and R. > I just came across the matrixStats package, > ## EXAMPLE #1 > lx <- c(1000.01, 1000.02) > y0 <- log(sum(exp(lx))) > print(y0) ## Inf

R Version 2.2.0 Platform: Windows When I use barplot but select a ylim value greater than zero, the graph is distorted. The bars extend below the bottom of the graph. For instance the command produces a problematic graph. barplot(c(200,300,250,350),ylim=c(150,400)) Any help would be appreciated. Paul [[alternative HTML version deleted]]

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080611/2da35c97/attachment.pl>

Hi, I am trying to write the Realized GARCH model with order (1,1) The model can be describe bellow: r_t = sqrt( h_t) * z_t logh_t = w + b*logh_(t-1) + r*logx_(t-1) logx_t = c + q*logh_t + t1*z_t +t2*(z_t ^2 -1) + u_t and z follow N(0,1) , u follow N(0, sigma.u^2) But I'm troubled with the simulation check for my code. After I simulate data from the model and estimate the data, I

Hi List. It may be a very simple question, but I couldn’t find an answer on the internet. Which function (maybe in a specific package) would I use to perform a Cochran C Test for Homogeneity of Variances? Stats package have the mantelhaen.test, but I’m “almost sure” that It’s not what I want, it’s probably the Cochran Q Test. So, any ideas? Thank you for the pacience.

Hi! I am searching for the Cochran-Armitage-trend-test. Is it included in an R-package? Thank you! --

Dear R users, I am happy to announce the new package `samplesizeCMH` is now available on CRAN <https://cran.r-project.org/package=samplesizeCMH>. The `samplesizeCMH` package will compute power and sample size for the Cochran-Mantel-Haenszel test for stratified 2x2 tables. There are also several helper functions for working with 2x2 tables, such as converting a proportion into an odds

Dear R-enthusiasts, I am trying to do a Cochran-Armitage test for trend in R. After consulting google I found Torsten Hothorn's remark that the 'coin' library could be used. lungtumor <- data.frame(dose = rep(c(0, 1, 2), c(40, 50, 48)), tumor = c(rep(c(0, 1), c(38, 2)), rep(c(0, 1), c(43, 7)),

Dear List, I am looking for what B.S.Everitt refers to as Cochrane Method for testing independence in combined 2x2 contingency tables. Is it the same method as the Cochran-Mantel-Haenszel Chi-Squared Test for Count Data in R? Thanks, Serguei [[alternative HTML version deleted]]

Hello all I wish to perform a mcnemar.test() for a 5X5 matrix. Wikipedia tells me (http://en.wikipedia.org/wiki/Cochran_test) I should turn to cochran.test. The only place I found it was in the "outliers" package, but the command cochran.test() acts differently then mcnemar.test() , and doesn't take a table as input. Any ideas on how to use it ? #Example code: aa =

Dear R users, I am happy to announce the new package `samplesizeCMH` is now available on CRAN <https://cran.r-project.org/package=samplesizeCMH>. The `samplesizeCMH` package will compute power and sample size for the Cochran-Mantel-Haenszel test for stratified 2x2 tables. There are also several helper functions for working with 2x2 tables, such as converting a proportion into an odds

Hello, I've tried to analyze some data with a CMH test. My 3 dimensional contingency tables are 2x2xN where N is usually between 10 and 100. The problem is that there may be 2 strata with opposite counts (the 2x2 contigency table for these are reversed), producing opposite odds ratios that cancle out in the overall statistics. These opposite counts are very important for my analysis, since