similar to: T2 hoteling

Dear ALL I cannot run this line stat.obs <- apply (GS, 2, function(z) Hott2(t(DATA[which(z==1),]), cl)) Error in Hott2 (t(DATA[which(z == 1), ]), cl) : subscript out of bounds I will be glade if you guide me. ******************************************************************************* *GS is a matrix 1857*200 *DATA is a matrix 1857*79 *cl <- as.factor(y)

We were looking at conducting a pilot program with Dell to virtualize our lab environments. We are specifically looking at Dell/Citrix provisioning servers. The OS in this setup is streamed from the server on every boot. I was curious if the R Project had any licensing stipulations in this type of setup. This is in the beginning phases and we wanted to answer some questions about licensing.

Dear All I can not download R for Linux. I do not know which file I should install? Best Regards, Soheila

I'm sure this has been discussed many times, but I have a question about hoteling. My understanding would be this: A phone sitting on a desk. A user hits 9000 and it asks what extension you'd like to become. You type "1001" and then it asks for your password. You type 1234, and it says you're "logged in". You now are accepting calls at your phone and you're

Dear All Is there any function for "median test" in R? Best Regards, Soheila [[alternative HTML version deleted]]

Hi, I want to do a 2 sample hotelling's test but i can't figure out how. When i type T2.test it says there is no such test and when i tried library(rrcov) it says there is no such program. Cheers. -- View this message in context: http://r.789695.n4.nabble.com/How-to-do-Hotelling-s-t2-test-tp4128748p4128748.html Sent from the R help mailing list archive at Nabble.com.

I want to calculate the probability that a group will include a particular point using the squared Mahalanobis distance to the centroid. I understand that the squared Mahalanobis distance is distributed as chi-squared but that for a small number of random samples from a multivariate normal population the Hotellings T2 (T squared) distribution should be used. I cannot find a function for

Hi all, I've got the following matrix : ? mat <- matrix(rnorm(700), ncol=5, dimnames=list( paste("f", c(1:140), sep="_"), c("A", "B", "C", "D", "E"))) I can see that currently most of the multivariate Hotelling T2 tests are limited for application on two groups/samples. I wud appreciate if someone can provide me a

Dear All How can I transform 1 column to 2 columns in R? 211217_s_at GO:0005249 211217_s_at GO:0005251 211217_s_at GO:0005515 211217_s_at GO:0015271 211217_s_at GO:0030955 211217_s_at GO:0005249 211217_s_at GO:0005251 211217_s_at GO:0005515 211217_s_at GO:0015271 211217_s_at GO:0030955 Best Wishes, Soheila [[alternative HTML version deleted]]

var1 var2 var3 var4 var5 var6 var7 var8 var9 var10 gold 2 3 1 2 4 0 1 4 4 3 2 2 4 2 4 3 4 2 4 4 4 2 3 3 0 0 4 1 0 2 4 4 2 1 4 0 3 2 0 0 2 4 4 2 3 4 0 2 2 0 0 0 3 4 2 2 2 3 2 2 0 0 0 2 4 2 2 4 1 1 2 0 0 3 3 3 2 3 4 1 4 0 0 0 0 3 4 2 3 1 0 2 2 1 0 2 3 3 2 0 3 1 1 1 1 2 1 2 3 2 1

Thank you Matthieu for your helpful suggestions. Unfortunately I still have problems. I have tried to compile it via your suggestion. " this is strange... you should have the usual summary... I have on my machine library(urca) test2<-ur.df(nottem, type="none", lags=1) summary(test2) " When I type in "library(urca)" I receive the following. "Error in

Dear R-users I've simulated a longitudinal multivariate normal data set from which I've simulated missing-patterns such as MCAR MAR and a simple kind of non-MAR. I've imputated the values so I now have 'complete' data sets. I'm trying to perform a T2-test as done in the multivariate case under th enormal assumption. Is there something I've to think about when performing

I've been working with the ARM backend's disassembler and working it out pretty heavily over the past several months, and have run into a bug involving the LDR instruction and all variants which use the T2I_ld multiclass. Codegen doesn't seem to trigger this issue, because either the correct variant is always selected, or certain combinations of parameters never occur, I haven't

Dear All I want to chose just spacial columns in R. (read table) data1<- read.table("/home/Documents/data.txt",header=F,sep = "\t", as.is =F) data.2<-data1[sub.data[,1],] The dimension of data.2 is correct but my data are transformed NA head(data.2) V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 NA <NA> NA NA NA

Hi All, I'm new to R and am trying to run a unit root test on the vector "y" (a time series of inflation (i.e. changes in the Consumer Price Index quarter on quarter)). I've run the Augmented-Dickey-Fuller Test below (R's URCA package). It gives me an error that it cannot find the function ur.df unless I comment out the third last line of code (see below). I try to call

Hi I spent hours looking over my formula. Somehow I cant find the reason why it gives me different answer. help appreciated. x = as.matrix(read.table("http://www.niehs.nih.gov/research/atniehs/core/microarrays/docs/heinloth.txt",1)) x = t(x) #now rows are subjects, cols are genes x = x[order(rownames(x)),] #order by treatment group oxygen, ultra-violet, gamma radiation y =

Dear All I want to arrange our output in 3 columns in R. Example: N 04320 03010 00514 00604 00563 03430 03022 V 22.000000 0.000000e+00 15.0000000 15.0000 13.0000 1.300000e+01 51.0000000 p 0.640625 2.910383e-11 0.4257812 0.4375 0.6875 3.356934e-04 0.2434998 N 03060 04140 03050 00562 00565 V 20.0000000 62.00000000 273.0000000

Dear All I have 54000 plots in R, How can I observe them? If I‌ have to save them one-by-one? Soheila [[alternative HTML version deleted]]

By accident, I left out the lines defining n1 and n2. Here it is as a function. Peter B. hotelling <- function(d1,d2){ k <- ncol(d1) n1 <- nrow(d1) n2 <- nrow(d2) xbar1 <- apply(d1,2,mean) xbar2 <- apply(d2,2,mean) dbar <- xbar2-xbar1 v <- ((n1-1)*var(d1)+(n2-1)*var(d2))/(n1+n2-2) t2 <- n1*n2*dbar%*%solve(v)%*%dbar/(n1+n2) f <-

Unable to resolve error in seq.default(which(text == "")[1] + 1, length(text), 1): wrong sign in 'by' argument I am trying to run this code get.msg <- function(path) { con <- file(path, open = "rt", encoding = "latin1") text <- readLines(con) # The message always begins after the first full line break msg <- text[seq(which(text ==