similar to: Hotelling T2 test extension for multigroup data

Dear R fans, I have got a difficult sounding problem. For fitting a linear model using continuous response and then for re-fitting the model after excluding every single variable, the following functions can be used. library(MASS) model = lm(perf ~ syct + mmin + mmax + cach + chmin + chmax, data = cpus) dropterm(model, test = "F") But I am not sure whether any similar functions is

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 friends, I have a matrix with following format. group var1 var2 .......varN c1 group1 1.2399 1.4990....-1.4829 c2 group4 0.8989 0.7849.....1.8933 ... ... c100 group10 ..... I want to draw a profile plot of each condition c1 to c100, which rows in above matrix and each line representing a row should be uniquely colored according to the group(1 to 10). I think this is simple task

Hi naive question. It is possible to get R command for omitting rows or cols with missing values present. But if i want to omit rows or cols with i.e . >20% missing values, I could´t find any package-based command, probably because it is too simple for anyone to do that manually, though not for me. Can anyone please help me ? - vickie [[alternative HTML version

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.

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 <-

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 =

Good morning I've made an PCA and I'd like to plot a confidence region based on Hotelling T2? Does anyone know how to compute it? Thank you -- View this message in context: http://r.789695.n4.nabble.com/hotelling-and-confidence-region-tp3524204p3524204.html Sent from the R help mailing list archive at Nabble.com.

Hello R-users, I've been looking for a function performing one and two sample Hotelling test for testing equality of mean vectors. Has anyone implemented such a function in R? thanks a lot, Bill ============== Bill Donner Statistician

Hi My setup : Trixbox 2.6.1 & TE410P running well ..I've 2 design issues to consider : 1. I need to store the CallerId of the PSTN caller with his language preference so that next time he is played the prompt in his language that he chose the first time.What would be better - storing his number in the Asterisk DB and using Dbput and DBget ? or storing it in MySQL from the dial plan and

Dear All It is very kind of you to guide me. When I want to run this line, I see this error stat.obs <- apply(GS, 2, function(z) Hott2(t(DATA[which(z==1),]), cl)) Error in colSums(w * x) : 'x' must be an array of at least two dimensions cl <- as.factor(y) GS: a matrix with 0 or 1 GS: gene sets -> a data matrix with rows=genes, columns= gene sets, GS[i,j]=1 if gene i in

Hi, I am a graduate student at Stanford University and I have a general statistics question. What exactly is the difference between doing a two-factor repeated measures ANOVA and a Hotelling T-squared test for a paired comparison of mean vectors? Given: Anova: repeated measures on both factors, 1st factor = two different treatments, 2nd factor = 4 time points, where you are measuring the blood

Hello everyone. In my last post I did not explained my problem quite well. I made a principal component analysis and took the 2 first principal components. I made ​​a chart of my points based on the score of the 2 PC. I would like to add on this graph a 95% confidence region. To do this I used the ellipse function as follows: pcsref=PC$score[data[,1]==ref,1:2] #matrix containing the scores

Friends, I installed MS office 2k7 with wine 1.1.15. I am able to open the applications and work on it. However if i try to copy a figure from powerpoint and paste is on the word, it is not working. The reverse is also not working. In addition, i am not able to move the figures in the work file to the appropriate positions. I dnt understand what is going wrong. Any suggestion on how to resolve

Hi I have a couple of functions that work on the object created by another R command and then print out or summarise the results of this work. The main function is defined as: hotelling.t <- function(obj) { #internal commands } I then have print.hotelling.t() that takes the list returned by hotelling.t and prints it with some extra significance calculations, formatting, etc. I want to

Hello Everybody I am using Trixbox 2.4 (with TE420P & PRI lines) .. my setup is like Calls -->Asterisk-->legacy pbx--->analog extensions(agents). Whenver a call comes in , asterisk dials the ACD number of the legacy pbx which in turn decides to route to appropriate agent.. for ex : s,1,Dial(ZAP/g4/5432) [g4 is the 4th span and 5432 is the ACD number of legacy pbx under which agents

Hi, Me a newbie in RoR - switching from .NET. i want to write an application that has a newsletter functionality where I can create an interface where in I can have options to send newsletters to multiple people in the database. i need to make a form which has from hardcoded, to is selected from te database maybe a listbox next to it and then a subject and text area. how wud i go about doing

Hi friends, I have got a list where each element might have variable number of members. $`4213` [1] "214077_x_at" $`164832` [1] "225996_at" "235977_at" $`339010` [1] NA $`23410` [1] "221562_s_at" "221913_at" "49327_at" $`285386` [1] "229764_at" $`2099` [1] "205225_at" "211233_x_at"

Hello, I am using the SEM package in R to fit a multigroup latent variable model and ran into some difficulties. I have 2 questions: 1. First, I am getting the following error message and wondering what to do to fix it: Error in solve.default((N[g] - 1) * robustVcov(mod.g, adj.obj = adj.objects[[g]])) : system is computationally singular: reciprocal condition number = 4.52055e-23

Dear R users; Is there any function to plot the confidence ellipse (Hotelling's T^2) in an score plot from a PCA? I have the scores off course , say scores for PC1 and PC2, and the value of the Hotelling's T^2 statistic. Thanks for any hint PM