similar to: stats/debugging question hotelling t-sq

I want to install package ICSNP but I got the following message. An instructor of my class even colud not figure out what's going on. Does anybody know about this error? > install.packages("ICSNP") --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.4 Warning in

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

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

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.

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

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

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

Hello,   Anyone familiar with robust ANOVA in R? Please help.   I am conducting a between-between-within (2 X 2 X 2) ANOVA, with the focus on the significance of the three-way interaction. Because of very uneven sample sizes, and the violation of covariance homogeneity, I need to use some sort of robust ANOVA method, and I prefer bootstrap. I am aware of Rand Wilcox's robust package that can

Dear R useRs, The new contributed packages 'ICS' and 'ICSNP' are available on CRAN. Descriptions: The 'ICS' package implements the 2 scatter matrix transformation to obtain an invariant coordinate system or independent components, depending on the underlying assumptions. The result of the transformation is an object of the S4 class ics which is provided by this

Dear R useRs, The new contributed packages 'ICS' and 'ICSNP' are available on CRAN. Descriptions: The 'ICS' package implements the 2 scatter matrix transformation to obtain an invariant coordinate system or independent components, depending on the underlying assumptions. The result of the transformation is an object of the S4 class ics which is provided by this

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

On the Squirrelmail website I setup in the FAQ a quick install guid on how to get squirrelmail to work with Dovecot. Maybe someone can verify my setup procedure really quick. Maybe someone thinks this is wrong. Also if you have something to add please do so. No one has bothered to do this before so I took it upon my self. Here is the link in the FAQ its about 8 lines down

I've been using the nlme library to construct lme models and I'd like to extract an r-sq value for the model - how do I do this? With glm I just used the reported (null deviance-deviance)/(null deviance). Deviance is not reported directly for lme nor can it be extracted using the deviance() function. How do I work out the deviance or can I work out R-sq in some other way? Thanks, Yvonne

Hi, I encountered a very strange error with R (2.14.1). I tried to build an R package and eventually worked (passing all the checks etc.). Then, I tried this command "R CMD INSTALL --binary pkg", which supposedly "installs and produce a binary source archive (*.zip) for use on Windows only." Nothing happened, but I could not run "R CMD check pkg" anymore. Further,

Dear R-Users, I have a trivial problem, but extensive googling and ??'ing did not solve it: I want to obtain the sums of squares from a summary.aov() object for later use. Example: > DV <- rnorm(100) > IV1 <- as.factor(rep(c("male", "female"), times = 50)) > IV2 <- as.factor(rep(c("young", "old"), times = 50)) > >

Hi, On Tue, 2023-03-07 at 09:49 +0800, Xuan Zhuo wrote: > On Mon, 6 Mar 2023 12:58:22 -0500, "Michael S. Tsirkin" <mst at redhat.com> wrote: > > On Mon, Mar 06, 2023 at 12:15:33PM +0800, Xuan Zhuo wrote: > > > If the queue of xdp xmit is not an independent queue, then when the xdp > > > xmit used all the desc, the xmit from the __dev_queue_xmit() may

On Wed, Mar 08, 2023 at 04:13:12PM +0800, Yunsheng Lin wrote: > On 2023/3/8 15:14, Xuan Zhuo wrote: > > On Wed, 8 Mar 2023 14:59:36 +0800, Yunsheng Lin <linyunsheng at huawei.com> wrote: > >> On 2023/3/8 10:49, Xuan Zhuo wrote: > >>> If the queue of xdp xmit is not an independent queue, then when the xdp > >>> xmit used all the desc, the xmit from

Hi all: Can someone advice me on how to hold the residuals Mean sq value on a string so it can be used in other calculations. I was trying something like this: Msquare<-dfr$Mean sq but fails..Thanks dfr <- read.table(textConnection("percentQ Efficiency 1.565 0.0125 1.94 0.0213 0.876 0.003736 1.027 0.006 1.536 0.0148 1.536 0.0162 2.607 0.02 1.456 0.0157 2.16 0.0103

Hello R users, I am trying to run a three factor ANOVA on a data set with unequal sample sizes. I fit the data to a 'lm' object and used the Anova function from the 'car' package with the 'type=III' option to get type III sums of squares. I also set the contrast coding option to 'options(contrasts = c("contr.sum", "contr.poly"))' as

If the queue of xdp xmit is not an independent queue, then when the xdp xmit used all the desc, the xmit from the __dev_queue_xmit() may encounter the following error. net ens4: Unexpected TXQ (0) queue failure: -28 This patch adds a check whether sq is full in XDP Xmit. Thanks. Xuan Zhuo (2): virtio_net: separate the logic of checking whether sq is full virtio_net: add checking sq is full