search for: y12

...code into R, and gauss has a matrix product function called the horizontal direct product (*~), which is some sort of variant on the Kronecker product. For example if x is 2x2 and y is 2x2 the horizontal direct product, z, of x and y is defined (in the Gauss manual) as: row 1 = x11*y11 x11*y12 x12*y11 x12*y12 row 2 = x21*y21 x21*y22 x22*y21 x22*y22 Or in R code if: x <- matrix(seq(1,4,by=1),2,2, byrow=TRUE) y <- matrix(seq(5,8,by=1),2,2, byrow=TRUE) The resulting matrix, if I had an operator, would be the following matrix z, here formed in a contrived manner: z.1 <- c(5, 6...

...ements all of which have the same group indices (thus the same x values in the errbar plot). In other words, my data looks like this: X1 (group index), Y11 (measurement), error of Y11 (error in measurement) X2, Y21, error of Y21 ... Xn, Yn1, error of Yn1 // End of the first set of measurements X1, Y12, error of Y12, X2, Y22, error of Y22, ... Xn, Yn2, error of Yn2 // End of the second set of measurements X1, Y13, error of Y13, X2, Y23, error of Y23, ... Xn, Yn3, error of Yn3 // end of the third set of measurements I need to put one errbar plot that include all three cases with different symbol...

Hi, We got a question about interpretating R-suqared. The actual outputs for a test dataset is X=(x1,x2, ..., xn). model 1 predicted the outputs as Y1=(y11,y12,..., y1n) model n predicted the outputs as Y2=(y21,y22,..., y2n) ... model m predicted the outputs as Ym=(ym1,ym2,..., ymn) Now we have two ways to calculate R squared to evaluate the average performance of committee model. (a) Calculate R squared between (X, Y1), (X, Y2), ..., (X,Ym), and t...

Hi all, I have happened to work on MS .NET for sometime now, and I found that this language offers RNG what is called as Donald E. Knuth's subtractive random number generator algorithm (found here: http://msdn.microsoft.com/en-us/library/system.random.aspx#Y12). ? Here I was wondering whether R also have same RNG in it's inventory, so looked at ?set.seed. There I found 2 related RNGs namely 'Knuth-TAOCP-2002', 'Knuth-TAOCP'. Can somebody guide me what is the most related RNG with .Net's? ? I also want to have R to draw random numb...

...pattern: x1 y11 x2 y21 y22 y23 x3 y31 y32 ... I.e. I have an x followed by a few y's. What I would like to do is turn this wide format into a tall format with two columns: "x", "y". The structure is that xi needs to be associated with yij (e.g. x1 should next to y11 and y12, x2 should be next to y21, y22, and y23, etc.). x y x1 y11 x2 y21 x2 y22 x2 y23 x3 y31 x3 y32 ... I have looked at ?reshape but I didn't see how it could work with this structure. I have a solution using nested for loops (see below), but it's slow and not very efficient. I would lik...

Hi all, I have a multiple-linear regression problem. There are 13 columns of data, the whole data matrix is: n x 13, where n is the number of samples. Now I want to regress EACH of the first 12 columns onto the 13th column, with 2-parameter linear model y_i = b0 + b1 * x_i, where i goes from 1 to n, and b0 is the intercept. How do I create a design matrix to do the 12-column regression

...2 1 2 X115 X116 X21 X22 X23 X24 X25 X26 X27 X28 X29 X210 X211 X212 1 2 0 1 0 1 1 1 0 1 0 1 0 1 X213 X214 X215 X216 Y1 Y2 Y3 Y4 Y5 Y6 Y7 Y8 Y9 Y10 1 1 0 1 1 2 3 4 5 6 7 8 1 2 Y11 Y12 Y13 Y14 Y15 Y16 3 4 5 6 7 8 So, e.g. for a 3 column matrix: > matrix(do.call(c,imp), ncol=3) [,1] [,2] [,3] [1,] 1 0 1 [2,] 2 1 2 [3,] 1 0 3 [4,] 2 1 4 [5,] 1 1 5 [6,] 2 1 6 [7,] 1 0 7 [8,] 2...

Hello everyone, I want to transform a data.frame into an array (lets call it mydata), where: mydata[[1]] is the first imputed dataset...and for each mydata[[d]], the first p columns are covariates X, and the last one is the outcome Y. Lets assume a simple data.frame: Imputed = data.frame( X1 = c(1,2,1,2,1,2,1,2, 1,2,1,2,1,2,1,2), X2 =

...2 1 2 X115 X116 X21 X22 X23 X24 X25 X26 X27 X28 X29 X210 X211 X212 1 2 0 1 0 1 1 1 0 1 0 1 0 1 X213 X214 X215 X216 Y1 Y2 Y3 Y4 Y5 Y6 Y7 Y8 Y9 Y10 1 1 0 1 1 2 3 4 5 6 7 8 1 2 Y11 Y12 Y13 Y14 Y15 Y16 3 4 5 6 7 8 So, e.g. for a 3 column matrix: > matrix(do.call(c,imp), ncol=3) [,1] [,2] [,3] [1,] 1 0 1 [2,] 2 1 2 [3,] 1 0 3 [4,] 2 1 4 [5,] 1 1 5 [6,] 2 1 6 [7,] 1 0 7 [8,] 2...

...for gene expression divergence and 14 matrices which stands for gene sequence divergence. I have tried joining them by using the concatanation function giving SequenceDivergence <- c(X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14) ExpressionDivergence <- c(Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14) where X1,X2..X14 are the expression matrices containing r-values and Y1,Y2..Y14 are the ones with patristic distances Now, I want to plot SequenceDivergence vs. Expression Divergence Tried doing that using plot (Sequence Divergence vs. Expression Divergence) But then getting the erro...

...for gene expression divergence and 14 matrices which stands for gene sequence divergence. I have tried joining them by using the concatanation function giving SequenceDivergence <- c(X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14) ExpressionDivergence <- c(Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14) where X1,X2..X14 are the expression matrices containing r-values and Y1,Y2..Y14 are the ones with patristic distances Now, I want to plot SequenceDivergence vs. Expression Divergence Tried doing that using plot (Sequence Divergence vs. Expression Divergence) But then getting the erro...

...Drafts: >>> y10 SELECT Drafts (QRESYNC (1381746443 1)) <<< * FLAGS (\Answered \Flagged \Deleted \Seen \Draft) <<< * OK [PERMANENTFLAGS (\Answered \Flagged \Deleted \Seen \Draft \*)] Flags permitted. Flag is \Draft instead of \Drafts? When I select Sent >>> y12 SELECT Sent (QRESYNC (1381746442 7 (2 2))) <<< * OK [CLOSED] Previous mailbox closed.?? <<< * FLAGS (\Answered \Flagged \Deleted \Seen \Draft $SENT) <<< * OK [PERMANENTFLAGS (\Answered \Flagged \Deleted \Seen \Draft $SENT \*)] Flags permitted. Why is the flag $SENT ? ht...

...?9.2?44.3?0.02 Y6?5.9?18.6?37.4?14.5?0.3?36.9?8?9.5?32.9?0.10 Y7?8.0?16.1?88.6?24.1?0.1?34.6?2?8.7?11.1?0.02 Y8?13.6?21.1?56.3?19.0?0.7?33.3?6?9..5?10.8?0.06? Y9?11.2?20.4?40.7?12.9?1.1?40.3?3?12.2?6.5?0.04 Y10?7.6?18.3?27.5?8.1?2.3?41.9?2?5.9?2.9?1.00 Y11?8.8?22.2?33.3?8.8?0.6?44.4?4?6.6?55.5?0.09 Y12?9.4?16.5?35.6?16.2?0.7?50.2?5?8..8?31.4?0.07 ? (These files have been given as attachments as well) The data file has 10 variables in all, divided into?two groups: ?G1 consisting of the first 4 variables and G2?of the remaining 6 variables. We wish to iteratively pick one variable from the group G...

After searching in R- Introduction, FAQ, help... I don't understand this: I write a function in a file (.R): tt <- function(mc) { date() mc<-read.csv2("machines.txt",na.strings="") date() } I source it in R and I type tt(). The answer is > tt() [1] "Mon Aug 14 11:18:25 2000" > The instructions following the first "date()" are ignored. Why?

Hi All, I am making a serious effort to try to learn R, but one hurdle I am facing is that I need to "see" the data as I walk through the examples in the packages. For instance, many examples on the web start by a command like data("wines"). How can I actually view what the dataset looks like prior to transformations and analysis? I have tried to use edit() , print, and

...6 > "y6",150,3,3,0,11,9,6,0,0,20,10 > "y7",340,12,8,6,26,19,0,0,0,20,25 > "y8",300,5,3,0,16,20,11,0,0,10,8 > "y9",380,8,11,5,10,25,9,2,40,30,10 > "y10",40,5,0,0,2,1,2,0,0,0,0 > "y11",230,12,9,0,11,13,4,0,0,0,10 > "y12",180,0,0,0,0,14,0,6,210,4,2 > > --------------060007010402070602040701-- > > > --__--__-- > > Message: 4 > Date: Sat, 7 Dec 2002 16:01:04 -0500 (EST) > From: Naoki Takebayashi <ntakebay at bio.indiana.edu> > To: r-help at stat.math.ethz.ch > Subject: [R]...