similar to: Drawing boxplots in pairs function

Dear R-help, I am trying to rename the variables in a dataframe, called 'T1A' here. Seems renaming was successful, but when I call one of the variable I got error message and I wanted to know why. The data frame contains 365 rows and 49 columns. I would like to name the first column `DATE` and the others T0.5, T1, T1.5,...,T24 (as this is a set of data collected every half hour for a

I'd like to model the relationship between m responses Y1, ..., Ym and a single set of predictor variables X1, ..., Xr. Each response is assumed to follow its own regression model, and the error terms in each model can be correlated. My understanding is that although lm() handles vector Y's on the left-hand side of the model formula, it really just fits m separate lm models. What should

Hi Everyone, I have some data that predicts both a nominal and ordinal response variable. I was wondering what packages in R would help me analyze the data? I was also curious if anyone could recomend me some textbooks that would help with the analysis of such data? I have the 5th edition of "Applied Multivariate Statistical Analysis" by Richard A. Johnson and Dean W. Wichern

Dear helpeRs, the following data set comes from Johnson/Wichern: Applied Multivariate Statistical Analysis, 6th ed, pp. 304-306. /X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3, 2))) Y <- structure(c(0, 2, 4, 0), .Dim = as.integer(c(2, 2))) Z <- structure(c(3, 1, 2, 8, 9, 7), .Dim = as.integer(c(3, 2)))/ I would like to compute Wilk's Lambda in R, which I know is 0.0385.

On Nov 12, 2008, at 3:37 PM, David Greene wrote: >>> Unfortunately, calling FunctionType::refineAbstractType(opaque, void >>> (...)) >>> doesn't work because RefineAbstractType doesn't recurse down into >>> the >>> pointee types. >> >> What types are you starting out with? If you had the equivalent of: >> >> %t1 = type

SUPERMIN_KERNEL and SUPERMIN_MODULES don't work to guestfish. Since guestfish sets --if-newer parameter to supermin, so the environment variables are not used under the following conditions. - the output directory exists and, - the dates of both input files and package database are older than the output To solve that, rebuild the output it when SUPERMIN_KERNEL or SUPERMIN_MODULES are

On Wednesday 12 November 2008 17:29, Chris Lattner wrote: > In fact, I strongly suspect that RefineAbstractType doesn't work if > you give it something that isn't Opaque. I'd recommend only using it > on Opaque, not random abstract types. Ok. > Sure, LLVM is more general than C. In LLVM, you can have a pointer > that points to itself, for example: > > %t =

On Wednesday 12 November 2008 18:58, Chris Lattner wrote: > On Nov 12, 2008, at 3:37 PM, David Greene wrote: > >>> Unfortunately, calling FunctionType::refineAbstractType(opaque, void > >>> (...)) > >>> doesn't work because RefineAbstractType doesn't recurse down into > >>> the > >>> pointee types. > >> > >>

I am learning glm function, but how do you fit data using exponential distribution with glm? In the help file, under "Family Objects for Models", no ready made option seems available for the distribution as well as for other distributions satisfying GLM requirements not listed there.

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

I use 1.7.0 version under Windows XP. Problem: When graphics device history recording function turned on, suppose I source a file containing lines: plot(graph1) plot(graph2) I see the graph2 in the graphic device, assuming graph2 is the last plot in the file. Now using "Previous" under history menu, I get graph1 as expected. But after that, using "Next" does not show graph2.

In MDI mode of Windows, is it possible to save window size and other parameter settings of a graphics device when it is activated, much like Rconsole file? By default the bottom part of a graphic device (windows()) is too big such that it is hidden under the horizontal scrolling and status bar of Rgui. (Of course I can use the vertical scrolling bar of Rgui to see the hidden part.) So I select

I'd like to sample multiple response values from a multivariate regression fit. For example, suppose I have m=2 responses (y1,y2) and a single set of predictor variables (z1,z2). Each response is assumed to follow its own regression model, and the error terms in each model can be correlated (as in example 7.10 of section 7.7 of Johnson/Wichern): > ex7.10 <- + data.frame(y1 =

Hi, we''re playing around with XEN (kernel 2.6.32) and REMUS. For us, sees this combination of very interesting and we would like to use the productive. However, the configuration is a great effort and often leads to problems. Until the first test position ran, it took several tries. As the further development of REMUS is planned? I''ve read that it will soon be over

Hello there, I have two data sets with 14 variables each, and wish to do the test for equality of their covariance matrices and mean vectors. Normally these tests should be done by chi square test (box provided) and Hotelling's T square test respectively. Which R functions could do this kind of test? I just find some functions could do for one dimension, but no for multidimension. Some one

I have written a few different summary functions. I want to calculate the statistics by groups and I am having trouble getting the output as a dataframe. I have attached one example with a small dataset that calculates summary stats and percentiles, I have others that calculate upper confidence limits etc. I would like the output to be converted to a dataframe with one of the columns as the

Hello, I am trying to plot two sets of boxplots together. These are estimates of two experiments and?seven?factors. The results of the two experiments I want to plot as boxplots stacked to each other. Therefore I plot first the results of the first experiment; and next with the add option the second set of boxplots. The boxplots are plotted at 'at = 1:7 - 0.15 for the first experiment and

Dear all - With the following code I get Boxplots and Stripcharts in one Plot: with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3)) with(InsectSprays, stripchart(count ~ spray, col = "red", vertical = TRUE, add = TRUE)) But the dots from the stripchart are plotted over the Boxes. Is there any possibility to have Stripchart and Boxplots side-by-side, i.e. to move boxplots and/or

Hii, I created some boxplots with this commands: x <-read.table(file="test.txt") x$group <- rep(1:8, each=5) boxplot(V3~gruppe, data=x) Now, I will connect the boxplots to each other to the min, max and median values. Can anybody help me how to do it ? greetings, J -- View this message in context: http://www.nabble.com/connecting-boxplots-tp21405459p21405459.html Sent from

I tried to draw a point on all boxplots for their means, I did: boxplot( Daten$weight~interaction(Daten$Dosis,Daten$sex, drop=TRUE)) means<-tapply( Daten$weight, Daten$Dosis, mean) points(means, pch=5, col="red", lwd=5) but only the boxplots for male got that point on them, its really weird because I don't think that I separated the sex in the codes above [[alternative HTML