similar to: Multiple Comparisons-Kruskal-Wallis-Test: kruskal{agricolae} and kruskalmc{pgirmess} don't yield the same results although they should do (?)

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

Dear R helpers, I have been using the loadings function from the multiv library and I get the typical output (see below). When I try to export these results to a file using a write.table() I get the following error message "Error in as.data.frame.default(x[[i]], optional = TRUE) : can't coerce loadings into a data.frame" Any idea why write.table is doing that and any

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

Hi all, After running kruskal.test I have got results (p<0,005) pointing to reject the hypothesis that the samples were draw from the same population. Howerver when I run the kruskalmc there are no significant differences in any of the multiple comparisons. Is that possible? Some clarification? Thanks, Humber <https://sites.google.com/site/humberandrade> [[alternative HTML version

Hi all I am trying to generate a normal unbalanced data to estimate the coefficients of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to estimate the variance components and their standard errors. Further, I am trying to use the likelihood ratio test to test H0: sigma^2_b = 0 (random effects variance component), and the t-test to test H0:mu=0 (intercept of the model Yij = mu

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Dear list! I work with multiple Kruskal-Wallis test (kruskalmc, package pgirmess), which evaluates differences in medians among groups (5 groups). A result of a test is significant differences among some groups, while median values are the same for 4 groups (using tapply). Why? p.s.: number of samples in groups vary from 50 to 4900. Thanks to all, OV .

Dear all, I run a kruskal wallis test and found significant results. Then, I conducted all pairwise comparisons and found no significant results. Could anyone please give me a hint as to why this happens or redirect me towards a specific web page where I can find more info? (I used alpha=5% and made no bonferroni or other correction for the pairwise comparisons) Thank you   Dr. Iasonas

Hello R users, Can any one please help me to find a way to build phylogenetic tree by R program from any distance matrix. Suppose I have a data like : MATRIX [1] '1' 0.0 [2] '2' 0.071 0.0 [3] '3' 0.036 0.286 0.0 [4] '4' 0.429 0.75 0.714 0.0 [5] '5' 0.679 0.179 0.214 0.536 0.0 [6] '6' 0.893 0.929 0.964 0.464 0.357 0.0 [7]

I have a table like the following. I want to fit Cm to Vm like this: Cm ~ Cl+Q1*b1*38.67*exp(-b1*(Vm-Vp1)*0.03867)/(1+exp(-b1*(Vm-Vp1)*0.03867))^2+Q2*b2*38.67*exp(-b2*(Vm-Vp2)*0.03867)/(1+exp(-b2*(Vm-Vp2)*0.03867))^2 I use nls, with start=list(Q1=2e-3, b1=1, Vp1=-25, Q2=3e-3, b2=1, Vp2=200). But I always get 'singlular gradient' error like this. But in SigmaPlot I can get the result. How

Hi there, I've got a small question : is there any post-hoc test for Kruskal rank sum test integrated in R ? I know that the Nemenyi test is one of the post-hoc that can be used, but there's no (to by knowledge) R function for it. What should I do ? Thanks, Etienne

Hi, I have searched and found a response to a question similar to mine but when I tried the code, R says it's not an actual function so I thought I'd ask here. http://r.789695.n4.nabble.com/file/n4633035/Cookies.csv Cookies.csv I have attached the data I am using. I am trying to look at two things: how moisture content changes over time, and how it changes along the length of a log

Hi, I would like to find out if the kruskalmc command in the pgirmess package controls for the family-wise type-I error rate. Thank you Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-713178 Fax: +357-22-590539 Honorary Research Fellow Department of

Hi I have a dataframe of genomic positions something like below: chr start end strand 1a 2a 3a 1b 2b 3b 1c 2c 3c 1d 2d 3d chr1 201199522 201199523 - 51 62 63 46 53 58 53 78 37 27 65 62 chr1 201199584 201199585 - 42 49 35 54 62 55 43 56 26 21 64 63 chr1 201199603

Hola. Llevo poco tiempo con R y estoy atascado en un problemilla, a ver si podéis ayudarme. Estoy trabajando con una base de datos sobre personalidad en diferentes especies de aves. Esta base tiene datos de más de 7000 individuos y me gustaría crear un bucle mediante el cual realizar un test de Kruskal Wallis por cada especie (quiero ver si los comportamientos son consistentes dentro de cada

I am interested in using "multcompLetters" after running "kruskalmc" but I'm a newbie and I'm not having luck figuring it out. I can run "kruskalmc" just fine, but after studying the documentation for "multcompletters" for a long time, I cannot figure out how to make it work. Any ideas? R input is below > >> #here is my data set >>

Hi, I have got two groups of samples; and for every row, I wish to calculate Kruskal-Wallis' p-value. In the example below, and the stars (****) show where I am struggling to design and put things together. Any help would be appreciated. myFile <- data.frame(Sample_1a = 1:10, Sample_1b = 2:11, Sample_1c = 3:12, Sample_2a=4:13, Sample_2b=7:16, row.names=LETTERS[1:10]) groups <- rep

Hi everyone, I am working on analyzing spectra which requires hundreds of kruskal wallis tests to be performed. Once significant results are found post tests are performed ect. My questions are: how do I perform all these kw tests and save just the pvalues in a table? Is there anyway to have R automatically perform post tests if a result of p<0.05 is found? Currently when I run the test

Hello, Is there like the kruskal wallis test in relation to ANOVA (no restrictions on normallity and variance homogenity) something (in R) for MANOVA? thanks -- Dr.Nicolaas Busscher Universit?t GH Kassel Nordbahnhofstrasse: 1a, D-37213 Witzenhausen Phone: 0049-(0)5542-98-1715, Fax: 0049-(0)5542-98-1713

Hi all, The Kruskal-Wallis test is a generalization of the two-sample Mann-Whitney test to *k* samples. That being the case, the Kruskal-Wallis test with *k*=2 should give an identical p-value to the Mann-Whitney test, should it not? x1<-c(1:5) x2<-c(6,8,9,11) a<-wilcox.test(x1,x2,paired=FALSE) b<-kruskal.test(list(x1,x2),paired=FALSE) a$p.value [1] 0.01587302 b$p.value [1]