similar to: test for equality of two data sets with multidimensional variables

I would like to test the hypothesis that the difference between pairs, for several variables, is zero. This is easily done separately for each variable with: lm(Y ~ rep(0, nrow(Y))) where Y is a matrix whose columns are the differences for each variable between pair members. However, I would like to get an overall probability across all variables from a Wilks or Pillai-Bartlett statistic as in

# Peter B. Mandeville kindly offered to send me code for Hotelling's T^2 # Test. Unfortunately there seems to be no route to his machine. # So i'm trying to reach him via the Mailing List. ------------------------------------------------------------------------ Sir, this morning i recieved your message about the availability of the code for Hotelling's Test. I hurried to find out

I don't find rbind, cbind or other bind that bind in other dimesions for multidimension array while keeping the result of the same dimension. For example, say, dim(x)=c(2,3,4), dim(y)=c(3,3,4). I want to have the result z, where dim(z)=c(2+3,3,4), when I bind x and y in the first dimenion. The bind can also be in the 2nd and 3rd dimesion. Could somebody let me know if this family of

Hi All, I am trying to reproduce this using R instead. [image: Full-size image (38 K)] I tried using the following code *SChla <- read.csv("SM_Chla_data.csv")* *Atlantis <- SChla[16:66,]* *head(Atlantis)* * * Seamount Station Depth Pico Nano Micro Total_Ch dbar Latitude Longitud 16 Atlantis 1217 Surface 0.0639 0.1560 0.0398 0.2597 2.082 -32.71450 57.29733

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

R-help, I have the following file I want to import to R (some lines removed) Calibrated CTD data for station:00280001 Calibrated:23/8 2001, Salinity Unsmoothed, Fluorescence Uncalibrated Maximum observed depth: 36 m QUAL has one digit for each of pressure, temp., sal. and fluor. QUAL=1:Uncal., QUAL=2:OK, QUAL=6:Interp., QUAL=9:No data

Dear All, I am running a simulation to obtain coverage probability of Wald type confidence intervals for my parameter d in a function of two parameters (mu,d). I am optimizing it using "optim" method "L-BFGS-B" to obtain MLE. As, I want to invert the Hessian matrix to get Standard errors of the two parameter estimates. However, my Hessian matrix at times becomes

I am doing a cluster analysis with hclust. I want to get hclust to output the Hotelling's T squared statistic for each cluster so I can evaluate is data points should be in a cluster or not. My research to answer this question has been unsuccessful. Does anyone know how to get hclust to output the Hotelling's T squared statistic for each cluster? Mike [[alternative HTML version

This is a proposal for either adding a new, or updating an existing command line option such that a diagnostic can optionally be produced for unused -D macros. Long-lived large projects with thousands of files and many contributors have a tendency to accumulate build options over time. As time passes some build options like macros become replaced, obsolete, or simply no longer used. At the same

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, In the Unix environment, I open a window by x11(). May I specify the position of this window by specifying the position of the top left of the window as in Windows environment? Or some other parameters can be used to do that? Thank you, Shengzhe

R 1.8.0 on Windows XP Professional. A huge THANK YOU to the R Team for this marvelous software. I am making lattice plots of oceanographic data. The usual layout does not conform to plotting conventions that marine scientists use when depth is the independent variable. Under those conventions, plots are made with the origin at the upper left, depth on the vertical axis (increasing as it

Dear R users, I am doing MCMC using Metropolis-Hastings. My model is bivariate-log-normal and the prior for variance-covariance is wishart distribution. I am wondering if there are some simple codes about how to get the density of Wishart distribution in my case ? Thanks in advance. Meihua [[alternative HTML version deleted]]

Hello, When I open a new device by "windows()", how to hide the 3 buttons on the top-right corner of this window, since I want to make these buttons not work, and the window can be closed only by "dev.off()" or "graphics.off()". Thank you, Shengzhe

Hello, I have a data set with 15 variables (first one is the response) and 1200 observations. Now I use pls package to do the plsr as below. trainSet = as.data.frame(scale(trainSet, center = T, scale = T)) trainSet.plsr = mvr(formula, ncomp = 14, data = trainSet, method = "kernelpls", model = TRUE, x = TRUE, y = TRUE) from the model, I wish to know the

I am adding -D k=v to yaml2obj, similar to clang -D. This makes it easy to generate {32-bit,64-bit} x {big-endian,little-endian} tests. --- !ELF FileHeader: Class: ELFCLASS[[BITS]] Data: ELFDATA2[[ENCODE]] Type: ET_DYN Machine: EM_X86_64 # RUN: yaml2obj -D BITS=32 -D ENCODE=LSB %s -o %t.32le # RUN: yaml2obj -D BITS=32 -D ENCODE=MSB %s -o %t.32le # RUN: yaml2obj

Hello, I wish to plot some figures in a window in turn, but the size of these figures is different, so how can I change the size of the window by resetting the parameters before each plotting? Thank you, Shengzhe

Hello, When I use "windows()" to open several windows, e.g. 4 windows, device 2, device 3, device 4 and device 5. Normally using "dev.set()" to specify the current active window, if there is a way to specify the active window by mouse click? When I click device 3, device 3 becomes active (shown on device header), and when I click device 4, device 4 becomes active, so on....

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

?The idea itself is indeed good. Regarding to escaping: I think we should have it. Imagine the following example (I've took it from D73828). --- !ELF FileHeader: Class: ELFCLASS[[BITS]] Data: ELFDATA2LSB Type: ET_EXEC Machine: EM_386 # RUN: yaml2obj %s --docnum=4 -D BITS=32 -o %t-32bit.o # RUN: yaml2obj %s --docnum=4 -D BITS=64 -o %t-64bit.o Without escaping it would