similar to: simulation of levene's test

Hello everyone, I want to calculate type I error rate for modified F statistic for one way robust anova. I need to find the j group trimmed mean and winsorized sum of squared deviations. Here I attached my code for j=2 to make it simple. Originally I have j=4. Hope you can help. I need to run it for 1000 times My problem is: i) the value of F-test obtain from my simulation below is in negative

Hello R-helpers, I want to ask your opinion since I am not so sure how to do it. This is regarding one part of my paper project and my situation is: Stage I I have 2 groups and for each group I need to compute the following steps; i) Generate 3 random numbers from normal distribution and square them. ii) Repeat step 1 for 15 times and at the end I will get 15 random numbers. I already done

Hello All, Anyone who has experience using Digium analog card clones from any of the following: 1. Zycoo 2. CTVON 3. Chinaroby 4. Etross 5. Immediate IT (IIT) 6. Realtone and can give review which one is good quality with easy configuration and error free running. Also since some of these manufacture only analog cards, does anyone have any experience using these in a single system with digital

I am trying to use Levene's test (of package car), but I do not understand quite well how to use it. '?levene.test' does not unfortunately provide any example. My data are in a data frame and correspond to 1 factor plus response. Could someone please give me an example about how to use the command levene.test(y, group) Thanks in advance, marta -- View this message in

Hi, I run this code to get the power of the test for modified bartlett's test..but I'm not really sure that my coding is right.. #normal distribution unequal variance asim<-5000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) n1<-20 n2<-20 n3<-20 mu<-0 sd1<-sqrt(25) sd2<-sqrt(50) sd3<-sqrt(100) g1<-rnorm(n1,mu,sd1) g2<-rnorm(n2,mu,sd2)

Hi, Can you please help me with this please? What I am trying to do is call a vector from R function and used in the new function So I create 4 functions with these arguments M11 <- function(TrainData,TestData,mdat,nsam) { ls <- list() I have few statments one of them is vectx <- c(,1,2,3,4,5,6,6) vectz <- c(12,34,5,6,78,9,90) and then................ ls(vectx=vtecx,vectz=vectz)

Dear all, For my research I need to perform a Levene's test on my data.I have not large experience with this test and with R,so I have some problems to use the leveneTest function. I read some previous posts on this topic but I cannot understand what I have to do. I have a five-columns data frame organised in the following way: number of observation, number of subject, score, type of

Hi, I want to perform pass through I/O from my domU. I guess I would need an additional network card for the same. Is there any preference on the second network card I should get or any one is fine? I have a x86_64 machine with Cent OS 4.4 running. Are there any additional considerations that I cannot see - Kindly let me know. Thanks, Asim _______________________________________________

Hi there, I want to do a nonparametric covariance analysis and I have tried to use the package "sm" function "sm.ancova" but it didn't work for me because I have more then one covariates (I have 18 covariates and 3 factors). I want to analyse for one factor (who has 13 levels) where the differences for my response are, using the explanatory (covariates). (The other two

Dear Mariano, See the function leveneTest() in the car package. I hope that this helps, John ----------------------------- John Fox, Professor Emeritus McMaster University Hamilton, Ontario, Canada Web: socialsciences.mcmaster.ca/jfox/ > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Marcelo > Mariano Silva > Sent: Monday, January

Hi, I''ve trying booting Xen Hypervisor PV with a centos5 template using kernel and initrd from mirror.centos.org/centos/5/os/i386. When i start PV this go to installation mode. However this should boot with template. Please help in this regard what Am i doing wrong. -- Best Regards Asim _______________________________________________ Xen-users mailing list Xen-users@lists.xen.org

Dear All, I have been using -- shorewall-4.0.14-1 on Centos 5.2 and its been working just great. i would like to monitor shorewall traffic ( i mean the allowed and dropped traffic passing through my firewall )with a browser from any location apprecite if someone could help me or advice me which software or add on utility would do the above. really apprecite thanks nd regards simon

Hello...i want to find the empirical rate for type 1 error using fixed trimmed mean. To make it easy, i'm referring to journal given by this website http://www.academicjournals.org/ajmcsr/PDF/pdf2011/Yusof%20et%20al.pdf. I already run the programme and there is no error in it but i got zero for the empirical rate of type 1 error. The empirical rate for the type 1 error given in the journal

Hello, I haven't found errors in your code. I implemented the test in the paper (the first, fixed symetric mean) and it also gives me zero Type I errors, when alpha = 0.05. Try to see the value of min(pv) or to plot the histogram of 'pv', hist(pv) and you'll see that there are no significant p-values, at that level. Anyway I'll continue to look at it, but my first

Hi all, I have a data set with missing value. I would like to estimate those missing value by using normal ratio method. Below is part of my data: AS BL Serdang Jhr Phg Target station 0 0.0 12.8 0.0 23.7 0.0 6 0.0 81.7 0.2 0.0 NA 0 1.5 60.9 0.0 0.0 15.5 1 13.0 56.8 17.5 32.8 6.4 4 3.0 66.4

Hi, All! To calculate Levene's Test for Homogeneity of Variance I use R Commander, and this is the output: > levene.test(Dataset$age, Dataset$sex) Levene's Test for Homogeneity of Variance Df F value Pr(>F) group 1 0.8739 0.3567 33 I am not sure what means "Pr(>F)"? Can anyone explain/translate this? Regards, Iurie Malai Department of Psychology and

Dear All I am trying to use Levene's test (of package car), but I do not understand quite well how to use it. '?levene.test' does not unfortunately provide any example. My data are in a data frame and correspond to 4 factors plus response. Could someone please give me an example about how to use the command levene.test(y, group) ? Thanks in advance, Paul

Dear All, I sent already before a message concerning Levene's test in nested model, but I didn't get any answer. Optimistically I hope to get an answer this time. I also point a new question related to the whole model, because I haven't find any sure answer if I am analysing it in a suitable way or not. I really have tried to do my homework. I have response variable (y) and four

Dear All, I am willing to test homogenity of variances using Levene's test. I know how to do it with simple model, but I am lost when dealing with a little bit more complex model. I have response variable (y) and four factors (a, b, c, d). One of these four factors (d) is nested within another factor (c). In addition, I would like to take into account only 2nd degree interactions in my

Hello, I notice that when I do Levene's test to test equality of variances across levels of a factor, I get different answers in R and SPSS 16. e.g.: For the chickwts data, in R, levene.test(weight, feed) gives F=0.7493, p=0.5896. SPSS 16 gives F=0.987, p=0.432 Why this difference? Which one should I believe? (I would like to believe R :) Ravi -- View this message in context: