similar to: Fwd: Bad \usage lines question

Dear colleagues in R, Has anybody implemented the 1) (Goodman & Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the asymmetric and symmetric forms), or is anyone aware that they might reside in some package? A search in the R archives does indicate that they are (somehow) part of the CoCo package, but I would rather not start learning how to transform my

Dear colleages in R, I have earlier been working with R in Linux, where reading in a table containing Scandinavian letters ("?", "?", and "?") in the header as part of variable names has not caused any problem whatsoever. However, when trying to do the same in R running on new MacOS-X (with an Intel processor) with the same original text table does not seem to

On Wed, 1 Aug 2007, Upasna Sharma <upasna at iitb.ac.in> wrote: > From: "Upasna Sharma" <upasna at iitb.ac.in> > Subject: [R] Goodman Kruskal's tau > > I need to know which package in R calculates the Goodman Kruskal's > tau statistic for nominal data. Also is there any implementation for > multiple classification analysis (Andrews at al 1973) in R?

Dear all, I've been delighted to just notice that Cohen's formulas for Effect Size 'w' and the associated power have been implemented in the 'pwr' package (thanks to St?phane Champely and others).. There is one aspect, though, that perplexes me. I'm doing some last minute post hoc analyses, meaning that my sample size (N=3404) has been long fixed, and I'm

I'm wondering how odds ratio is computed. I thought that it is (n11/n12)/(n21/n22), but it is not what fisher.test() computes. Could somebody let me know? > n11=3 > n12=1 > n21=1 > n22=3 > > n1_=n11+n12 > n2_=n21+n22 > > n_1=n11+n21 > n_2=n12+n22 > > x=rbind(c(n11,n12),c(n21,n22)) > > threshold=dhyper(n11,n1_,n2_,n_1) >

Hello everybody, I use the sweetpotato database included in R package: data(sweetpotato) This dataset contains two variables: yield(continous variable) and virus(factor variable). Due to Levene test is significant I cannot assume homogeneity of variances and I apply Welch test in R instead of one-way ANOVA followed by Tukey posthoc. Nevertheless, the problems come from when I apply posthoc

dear all, I want to perform a posthoc test for my ANCOVA: a1<-aov(seeds~treatment*length) With summary(glht(a1, linfct = mcp(treatment = "Tukey"))) R tells me: "covariate interactions found -- please choose appropriate contrast" How do I build these contrasts? Ideally, I would like to have the posthoc test for the ANCOVA including a block-effect

Dear R experts, I'm beginner. My question about ANOVA for unequal sample sized data should be obsolete but I can not clarify it. I have a dataset from 23 males and 18 females. I measured one condition('cond') with 4 levels. So I'd like to see main effect of gender, cond and gender by cond interaction and also postHoc test. (In fact, I have to do anova 90 times) * 1. Question

Hi, I have a R output that looks as follow: Rad:0 Rad1:2 Rad3:3 I want to make a new matrix that looks like : sample size is 2400 Variable    n11  n12 Rad            0     2400-0=2400 Rad1          2       2400-2 Rad3  3      2400-3   Thanks a lot for your time and help:) Best,Farnoosh Sheikhi [[alternative HTML version deleted]]

Hi everyone, I'm using the function 'HWE.exact' of 'genetics' package to compute p-values of the HWE test. My data set consists of ~600 subjects (cases and controls) typed at ~ 10K SNP markers; the test is applied separately to cases and controls. The genotypes are stored in a list of 'genotype' objects, all.geno, and p-values are calculated inside the loop over all

I did 2 zip files containing each: - a working tftp-root with syslinux 3.75 - a not working tftp-root wit syslinux 3.81 They have ben build and tested using the instructions given in my first message. The bug appears starting syslinux 3.81, NOT 3.80 (I tried too many versions) The zip files are available here: http://t.bugier.free.fr/tftp-root--fail.zip

Hello R-Users, I have been using (nnet) by Ripley to train a neural net on a test dataset, I have obtained predictions for a validtion dataset using: PP<-predict(nnetobject,validationdata) Using PP I can find the -2 log likelihood for the validation datset. However what I really want to know is how well my nueral net is doing at classifying my binary output variable. I am new to R and I

Hi, I have a genotype data for both case and controls and would like to calculate the HW p-value. However, since the number of one genotype is 0, I got wired result. Would someone help me to figure it out? Or confirm it's right? Thanks a lot. ============ > library( "genetics" ) NOTE: THIS PACKAGE IS NOW OBSOLETE. The R-Genetics project has developed an set of enhanced

Hello I'm working on my RIS server; I got a problem when I use pxelinux.0 to load setupldr.exe from windows XP 32bits (both home and corporate) Here is how I configure a working server, then how to modify it to reproduce the bug: ========================= I suppose you got -a dhcp server configured to load pxelinux.0 to client computer. -a working tftp server (mine is tftpd-hpa, with a rule

Dear R-users, Is there a function in R that classifies data in more than 2 groups using logistic regression/classification? I want to compare the c-indices of earlier research (lrm, binary response variables) with new c-indices obtained from 'multiple' (more response variables) logistic regression. Best regards, Stephanie Delalieux Department Biosystems M?-BIORES Group of Geomatics

*Hi all, * *Assume that I have the following data set with tow variables and I want count the number of observation with identical values * ** *x1 x2* * 1 1 * * 1 0 * * 0 1* * 0 1* * 0 0* * 1 1* * 0 1 * I want the following output ** * * *n1=3 # number of identical observation between x1 and x2 variables* *n2=4 # number of different observation* How do I do it in

Hello all, Quick question: I want to do posthoc contrasts for a linear mixed effects model. However, when trying to use the csimtest function in multcomp package I receive an error message saying it cannot find the function, even after installing and loading package multcomp. Any pointers would be greatly appreciated Daniel

gpxelinux.0 (3.73 and 3.74) can't work with microsoft RIS (startrom.n12),but pxelinux.0 (3.73 and 3.74) can work with microsoft RIS well. keeppxe can't work in 3.74.

Hi, I am doing a retrospective analysis on a cohort from a designed trial, and I am fitting the model fit<-glmD(survived ~ Covariate*Therapy + confounder,myDat,X=TRUE, Y=TRUE, family=binomial()) My covariate has three levels ("A","B" and "C") and therapy has two (treated and control), confounder is a continuous variable. Also patients were randomized to

Dear R experts, I'm new in R and a beginner in terms of statistics. It should be simple question, but definitely difficult to solve it by myself. I'd like to see main effect of group(gender: sample size is different(M:F=23:18) and one of condition(cond) and the interaction at each subset from 90 datasets So I perform anova 90 times using a command like below; for(i in 1:90)