similar to: need help to find type I error rate for modified F statistic

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

Dear List, I'm new in R. I'm trying to solve a simple constrained optimization problem. Essentially, let's say I have a matrix as in the object 'mm' inside the function below. My objective function should have a matrix of parameters, one parameter for each element 'mm' (4 in this case). The problem is to select the value of the parameters to maximize the function

Dear List, I'm trying to solve this simple optimization problem in R. The parameters are the exponents to the matrix mm. The constraints specify that each row of the parameter matrix should sum to 1 and their product to 0. I don't understand why the constraints are not satisfied at the solution. I must be misinterpreting how to specify the constrains somehow. library(alabama) ff <-

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

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)

hello, I try to run simulation of levene's test to find the p-value but the error of replacement has length zero occur, could anyone help me to fix this problem? asim <- 1000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) g1 <- rnorm(20,0,2) g2 <- rnorm(20,0,2) g3 <- rnorm(20,0,2) x <- c(g1,g2,g3) group<-as.factor(c(rep(1,20),rep(2,20),rep(3,20))) library(Rcmdr)

My apologies, I know that this is not a new problem, but I'm not sure how to find the answer I want to recursively loop over an object and trim trailing white space. When I use this function on a list of data.frame I get output like this: [1] "c(\" many spaces \", \" many spaces \")" "c(\" many spaces \", \" many spaces

Dear Colleagues, I hope to know how ordinal logistic regression with a mixed model is made in R. We (My colleague and I) are studying the behavior of a beetle. The attraction of beetles to a stimulus are recorded: the response is Slow, Mid, or Fast. They are based on the time after the presentation of the stimulus to the beetles. Because we do not observe the behavior continuously but do

1. The matrices are only 3x2. In your loop you'd be assigning rows 4 and 5 in the last two iterations. Are you sure that's what you want? 2. The reason it ``didn't work'' is because assign() takes the first argument as the name of the object to create, literally, instead of evaluating it. 3. You probably ought to be working with a list. E.g.: m <- replicate(5,

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, 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've sent this question 2 days ago and got response from Sarah. Thanks for that. But unfortunately, it did not really solve our problem. The main issue is that we want to use our own (manipulated) covariance matrix in the calculation of the mahalanobis distance. Does anyone know how to vectorize the below code instead of using a loop (which slows it down)? I'd really appreciate

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

Dear all, I am trying to do weighted regression using lm function in R. However, I have a question why the results from 1) lm(formula = Y~aX, weight = w) 2) lm(formula = wY~waX) are different. Aren't they supposed to have the exactly same result? Below are the R code to see difference in regression results MatY <- c(0.15,0.42,0.31,0.22) MatX <-

Hi, all. The help file for tapply says that if simplify is true, and the result of the calculation is always a scalar, then tapply will return a vector. Nonetheless: > t1 <- tapply(runif(10), rep(1:5, 2), mean) > is.vector(t1) [1] FALSE > is.array(t1) [1] TRUE > I have found this in version 0.65.1 on an SGI running Irix 6.5, and on a Mac running Linux-PPC. I've also

Hi, As it has been brought to my attention (by several parties), my current MB does not support VT-d. Since I don''t feel like purchasing a brand new system, I have searched for replacement MB-s that would fit into my existing uATX case, support my C2D CPU, and support at least 4x2GB ram. I have narrowed it down to the following four models: SUPERMICRO MBD-C2SBM-Q-O (Q35 chipset)

hello, I am attempting to compile a Dom0 CentOS 5.2 guest from source, using XenLinux 2.6.18.8 with a Xen 3.2.1 hypervisor. I am able to compile the Dom0 kernel, but when I attempt to boot into it, I receive the following messages: mount: could not find filesystem ''/dev/root'' setuproot: moving /dev failed: No such file or directory setuproot: error mounting /proc: No such