similar to: Coin Toss Simulation

Is there a way to estimate a nominal response model? To be more specific let's say I want to calibrate: \pi_{v}(\theta_j)=\frac{e^{\xi_{v}+\lambda_{v}\theta_j}}{\sum_{h=1}^m e^{\xi_{h}+\lambda_{h}\theta_j}} Where $\theta_j$ is a the dependent variable and I need to estimate $\xi_{h}$ and $\lambda_{h}$ for $h \in {1...,m}$. Thank you, Mauricio Romero Quantil S.A.S. Cel: 3112231150

All, I've devised a way to toss a coin three times and record the number of heads, then use the table command to summarize the number of time 0, 1, 2, or 3 heads were obtained. n <- 100 x <- sample(c("H","T"),3*n,replace=TRUE) dim(x) <- c(3,n) num_heads <- apply(x,2,function(x) sum(x=="H")) table(num_heads)/n What I'd like to do next is run this

Dear R-help, Could someone please try to explain this paradox to me? What is more likely to show up first in a string of coin tosses, "Heads then Tails", or "Heads then Heads"? ##generate 2500 strings of random coin flips ht <- replicate(2500, paste(sample(c("H", "T"), 100, replace = TRUE), collapse =

Well, I've now tried just about every permutation I can think of - including using the ATLAS lapack (instead of generic), in various combinations with ACML, or generic BLAS, or ATLAS blas. No matter what I do, so long as I have lapack of any flavour, I get the 'stats error' with make check. And, it always seems to throw the error at the same point in the stats test file (an svd in

Hello, sorry if this is a duplicate message -- the R mail-server told me that my original post was being held for moderator approval. Whatever. I just want to know if it is possible to connect R to a running process and have realtime data updates to some plots? Thanks, Charles [[alternative HTML version deleted]]

Hello, i'm trying to find a more modern system to reproduce the functionality that was available through the Histoscope program (from Fermilab). Namely, the capability of connecting to a running process and having plots update in realtime in response to new data. Is this possible with R? Thank you, Charles Cosse [[alternative HTML version deleted]]

Given the following data, and hypothesized median M.0 I've found a method to implement the Wilcoxon signed-rank test. Data: (with one zero difference and tied ranks) x <- c(136, 103, 91, 122, 96, 145, 140, 138, 126, 120, 99, 125, 91,142, 119, 137) M.0 <- 119 > library(exactRankTests) Package ?exactRankTests? is no longer under development. Please consider using package ?coin?

Hi, R allows me to run a one sample Wilcoxon test like this: wilcox.test(c(1,3.5,2.1,4,1.5,5), mu=2, exact=TRUE) The function 'wilcoxsign_test' from the package 'coin' should (I suppose) be able to calculate exact p values even if there are ties in the ranks. However, I couldn't find information on how to run a one sample test using 'wilcoxsign_test' like in the

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

Hi, I am trying to make use of "approximate" option to obtain null distribution through Monte-Carlo resampling, as described in coin library documentation. Unfortunately, this does not play well with my data -- permutation process swallows astonishingly large amounts of RAM (4-5Gb) and runs far too long (30 min for B=10). Apparently, this is caused by the size of my dataset (see

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

Hi all, is it possible to perform a linear-by-linear association test (lbl_test, package coin) for a SINGLY ordered table ? And, if yes, how to manage the "scores" ? Thanks in advance Gilles - France [[alternative HTML version deleted]]

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

ProductPublic Mobile Phone Multi-charger (With/Without Coin Module) Functions1. Accept coins. 2. Allows conversation while charging. 2. 10 travel charger plugs for most models. 3. Light box advertising. 4. UV phone sanitizing InputAC 110 - 250V 50~60Hz Output< 12 V 1C 5 A x 10 Measurement29 cm x 20 cm x 51.6 cm Net Weight5 kg Contacthttp://www.YukiCharger.com sales@YukiCharger.com

Conditional Inference Procedures in a Permutation Test Framework The `coin' package implements a general framework for conditional inference procedures, commonly known as permutation tests, theoretically derived by Strasser & Weber (1999). The conditional expectation and covariance for a broad class of multivariate linear statistics as well as the corresponding multivariate limiting

Conditional Inference Procedures in a Permutation Test Framework The `coin' package implements a general framework for conditional inference procedures, commonly known as permutation tests, theoretically derived by Strasser & Weber (1999). The conditional expectation and covariance for a broad class of multivariate linear statistics as well as the corresponding multivariate limiting

Hi, I am running the following example which can be found on page 12 of the pdf file of COIN package wt<-wilcox_test(pd~age,data=water_transfer,distribution="exact", conf.int=TRUE) "wt" actually contains the estimate of difference in location and the confidence interval of it. I am just wondering how can I extract these values? From the examples, I understand that the

Dear R-list members, This is in fact a question about statistics, not directly about the R software. The coin package, for conditional inference through permutation methods, has as it main function the function independence_test. One of its parameters is teststat, about which the package documentation says: teststat: "a character, the type of test statistic to be applied: either a

Dear expeRts, I would like to use a oneway_test (from package coin) to test whether two groups differ on various variables. The variables are encoded within a data frame. Unfortunately, I obtained an error, that I don't understand. Could you please help me ? Example: library(coin) y <- as.data.frame(matrix(rnorm(200), ncol=2)) group <- as.factor(unif(100)) lapply(y,

Good Evening R-Help Community: I have attached a file that contains the output from sessionInfo() and a summary of my Win XP system. I am?running R 2.12.0 and using Tinn-R 2.3.6.2 as my interface. When I attempt to call either the doBy or coin packages R generates an error that I do not understand and have so far not been able to resolve by searching R resources. I exchanged a couple of