similar to: coin package (conditional inference / permutation): parameter teststat

A general question- and excuse me if you find it irrelevant Are you tracking how many people finally download a specific R version by counting the number of downloads through log/ analytics software like GA etc? If so- can we see some numbers Ajay Websites- http://decisionstats.com http://dudeofdata.com Linkedin- www.linkedin.com/in/ajayohri

Dear R-list members, I am trying to put a legend just outside the plotting area of a graph. 1. Is there some way to plot symbols on the figure margin, in the same way as function "mtext" can write on the margin? 2. Is there some way to use the "legend" function (or some equivalent function) to plot a legend on the figure margin? 3. Any other way to plot either symbols or a

To the R-help list, In the documentation of the prop.trend.test function in the stats package, no bibliography has been provided which would allow one to find out the theoretical basis of that function and/or details of its implementation. May I suggest that some bibliography be included, as it generally happens with other statistical functions. I currently use R 2.8.0 running on Windows XP.

Dear R list members, How to produce barplots anchored to the x-axis (not floating above the x-axis) with a box around? With both following codes, the lower horizontal line of the box is below the y = 0 line: # first code x <- c(1,2,3,4) barplot(x,yaxs='i') box() # second code x <- c(1,2,3,4) op <- par(yaxs='i') barplot(x) box() par(op) The parameter yaxs='i'

Dear R users, The parameters xaxs and yaxs (function par, package graphics) seem not to work with the function hist (package graphics), even when the parameters xlim and ylim are defined. Is there any way to make yaxs="i" and xaxs="i" work properly with the function hist, mainly to produce histograms that "touch" the horizontal axis? The R documentation and the R

Dear R-list members, Once a gam (package gam) model has been fitted with family=poisson, is there some R function that could extract the diagonal elements of the smoother matrix S, to be used in a cross-validation for the selection of the best smoothing parameter, following equation 3.19 on page 48 of the book T.J. Hastie and R.J. Tibshirani, Generalized Additive Models, Chapman and Hall/CRC,

Dear R-list members, I have a data file with thousands of lines (cases), where each line contains the values of several variables. I would like to separate these lines in small groups, with each group followed by a blank line, to ease the visual inspection of the data in some situations. I am writing the output files with function write.fwf in library gdata, for correct column alignment. Below

Dear R-list members, I have just downloaded R version 2.11.0. Then I installed package gdata. At the end of the installation, a warning message said that package gtools could not be found. Package gdata was, after the installation, included in the "Package Index" in the help files. Messages during installation (which was made through the menu Packages / Instal package(s), choosing the

Dear R users, I need to save a graph in jpeg format. After plotting the graph, when the graphics window is active, in the File menu, the Save as / Jpeg / 100% quality correctly saves the graph in jpeg format. But I would like to know, how could I control the resolution (in dpi) of the saved jpeg file? I need to produce a jpeg at 1200 dpi. I have tried also the jpeg function in the package

Dear R helpers, I am trying to understand how to use the independence_test function in the coin package. I think I suffer from a misunderstanding about what the package does. Either that or I do not understand how to use it properly. Specifically, I cannot understand if I can test independence of arbitrary statistics. Take the following example: set.seed(10) d <- data.frame(y = c(rnorm(10,

Anyone know if coin can run a permutation test based on a (user-defined) statistic other than the mean difference? The function independence_test does the permutation t-test via difference in means. I'm wondering if it's possible to use independence_test to run a permutation test for some other statistic than the difference in means. For example, I'd like to run a permutation test

On 5/2/07, r-help-request@stat.math.ethz.ch < r-help-request@stat.math.ethz.ch> wrote: > > Send R-help mailing list submissions to > r-help@stat.math.ethz.ch > > To subscribe or unsubscribe via the World Wide Web, visit > https://stat.ethz.ch/mailman/listinfo/r-help > or, via email, send a message with subject or body 'help' to >

Dear R-list members, I would like to pose a question about the use and results of the glm() function for logistic regression calculations. The question is based on an example provided on p. 229 in P. Dalgaard, Introductory Statistics with R, 2nd. edition, Springer, 2008. By means of this example, I was trying to practice the different ways of entering data in glm(). In his book, Dalgaard

I am using "apply" to run exact permutation t-tests by columns using the coin package. For example: library(coin) dat <- matrix(rnorm(7*35),7,35) fun <- function(x) { pvalue( independence_test(x~f, data=data.frame(x, f=factor(c(rep("a",4),rep("b",3)))), distribution = "exact")) } p.vals <- apply(dat, 2, fun) Some small-scale

Dear List, I'm interested in performing a permutation test on the interaction between a binary treatment indicator and a covariate (either continuous or categorical). I'm interested in the p-value of the interaction effect from a permutation test, and I'm using the coin package for that purpose. As I haven't seen any examples like this in the package documentation (or anywhere

Dear sir, I am Shibu John from Thrombosis Research Institute India. It is a multidisciplinary organisation concerned with the interrelated problems of thrombosis and atherosclerosis. I was searching for Cochran armitage trend test program in R. Then I had seen your R coding for C-A trend test. I tried that in the R software. But I can?t run the program due the [Error: could not find function

Dear fellow R users, Keywords: Kruskal-Wallis, Post-Hoc, pair-wise comparisons, Nemenyi-Damico-Wolfe-Dunn test, coin package, oneway_test I am using the "oneway_test" function in the R package "coin" and I am obtaining results which I cannot believe are accurate. I do not wish to waste anyone's time and so if the following problem is rather trivial, I apologize, however I

Hello, Would you tell my how to extract a result from a test - it's justified because I need to run this test many times. Here is an example from authors' test: > library("coin") > 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)), rep(c(0, 1), c(33, 15)))) > ca.test<-independence_test(tumor ~

I partly answered my question since independence_test() function in coin package apparently do Cochran-Armitage trend test just like Eric Lecoutre's function tabletrend() - slightly modified here: > independence_test(pheno ~ geno, data = dat2, teststat = "quad", scores = list(geno = c(0, 1, 2))) Asymptotic General Independence Test data: pheno by groups 1 < 2

Using the package ?exactRankTests? one can execute a one-sample permutation test for a hypothesized location parameter of 0 like: perm.test(rnorm(30,0)) The package ?exactRankTests? seems now to be deprecated in favor of the ?coin? package which as I understand is a superset of ?exactRankTests? in terms of functionality. The ?coin? package allows one to run a two-sample permutation test using