similar to: Warning: matrix by vector division

Hi, I have a matrix, say m=matrix(c( 983,679,134, 383,416,84, 2892,2625,570 ),nrow=3 ) i can find its row/col sum by rowSums(m) colSums(m) How do I divide each row/column by its rowSum/colSums and still return in the matrix form? (i.e. the new rowSums/colSums =1) Thanks. Casper -- View this message in context:

Dear R users, Is there the possibitily in R under X11 to get (after typing help(...) command) separate help window, as it is in Windows version? Best wishes, ================================= Dr. Alexey B. Shipunov Section of Molecular Systematics Jodrell Laboratory Royal Botanic Gardens, Kew, Richmond, Surrey, TW9 3DS, U.K. e-mail: a.shipunov at rbgkew.org.uk

Dear R experts, I have a matrix (from some sort of classification) like this: object group [1,] 1 1 [2,] 2 2 [3,] 3 1 [4,] 4 1 [5,] 5 3 And I need something like this: [,1] [,2] [,3] [,4] [,5] [1,] 1 0 1 1 0 [2,] 0 1 0 0 0 [3,] 1 0 1 1 0 [4,] 1 0 1 1 0 [5,] 0 0 0 0 1 where all

Dear all, I have two data frames- x1 and y1 with same row names and column names(actually the names of the patients). x1 a b c d e a 1.0000000 0.4730679 0.6226994 0.6036036 0.6433333 b 0.4730679 1.0000000 0.6227273 0.6303855 0.5730858 c 0.6226994 0.6227273 1.0000000 0.7290503 0.6900585 d 0.6036036 0.6303855 0.7290503 1.0000000

Hi, I would like to construct a transition matrix from a data frame with annual transitions of marked plants. plant<-c(1:6) class<-c("seed","seed", "seed", "veg", "rep", "rep") fate<-c("dead", "veg","veg","rep", "rep", "veg") trans<-data.frame(plant, class, fate)

Hello R-Users, its again me with a question. I´m using R 2.5.0 on Mac Power Book running Mac OS X 10.4.10 I try to calculate distances betweeen two data tables looking like this V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 V25 V26 V27 V28 V29 V30 V31 V32 V33 V34 V35 V36 1 1 0 0 0 1 1 0 1 0 0 0 1 0 0 0 0 0 1 0 0

Hi, xtabs() creates a table of counts. I want a table of proportions -- that is, I want to divide every vector (along a particular dimension) by its sum. The tiny example below does that. The call to xtabs() creates a matrix "A" with dimensions ("x1","x2","y"). I transform "A" using aperm() and aaply() to get the matrix "B". The

Hello, might be rather easy for R pros, but I've been searching to the dead end to ... twsource.area <- table(twsource, area, useNA="ifany") gives me a nice cross tabulation of frequencies of two factors, but now I want to convert to pecentages of those absolute values. In addition I'd like an extra column and an extra row with absolute sums. I know, Excel or the

Dear list, Im trying to do the following operation but im not able to do it This is my table: 1 2 3 1 0 7 4 2 0 2 0 3 0 1 3 4 0 3 4 what i would like to do is divide each row values with the corresponding column' sum,namely: 1 2 3 1 0 0.54 0.36 2 0 0.15 0 3 0 0.08

Hello all! I have a problem with ggplot2 library. I want to do an heat map and the y variables are the year months. If I use the following code, he y values are in alphabetical order, but I want it in month order. The code is: library(reshape) library(ggplot2) library(scales) p <- ggplot(data.m, aes(variable, Month)) + geom_tile(aes(fill = value),

Dear all, I would like to ask one question related to statistics, for specifically on defining dummy variables. As of now, I have come across 3 different kind of dummy variables (assuming I am working with Seasonal dummy, and number of season is 4): > dummy1 <- diag(4) > for(i in 1:3) dummy1 <- rbind(dummy1, diag(4)) > dummy1 <- dummy1[,-4] > > dummy2 <- dummy1 >

Dear all, I am highly disappointed of the following decision: === ... Starting with R 2.3.0, CRAN binaries support Mac OS X 10.4 (Tiger) and higher only. === Of course, === It is, however, possible to compile binaries for earlied OS X versions from sources. === But the process, described on http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html#Building-R-from-sources consists of more than

Dear all, let say I have following data: dat <- structure(list(V1 = structure(c(1L, 4L, 5L, 3L, 3L, 5L, 6L, 6L, 4L, 3L, 5L, 6L, 5L, 5L, 4L, 4L, 6L, 2L, 3L, 4L, 3L, 3L, 2L, 5L, 3L, 6L, 3L, 3L, 6L, 3L, 6L, 1L, 6L, 5L, 2L, 2L), .Label = c("C", "G", "I", "O", "R", "T"), class = "factor"), V2 = c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L,

Dear R helpers: I have tried many times to find some way to read a list into R. But I faid. Here is an example: I have a file 'List.txt' which includes data as follows: [[1]] [1] 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.5 0.0 0.0 0.0 [19] 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 [[2]] [1] 0.0000000 0.5000000 0.0000000 0.0000000

Hello, I am trying to get a linear model of y ~ log(x). *> lm (y~log(x))* However, I always get an error report: /Error in model.frame.default(formula = y ~ log(x), drop.unused.levels = TRUE) : variable lengths differ (found for 'log(x)')/ *Here was my y:* > y [1] 0.4500000 0.0500000 0.5000000 0.4000000 0.0000000 0.5000000 0.4000000 [8] 0.0500000

Hi, I'm trying to find an implementation of monotonic regression in R and I haven't been able to find anything that's really related to this. isoMDS in the MASS package uses monotonic regression, however, I was wondering if there is any standalone function for monotonic regression? Basically what I'm trying to do is implement monotonic regression where I can see not just the

Hello everyone I'm trying to use Scheffe's complete cubic model (mixture design). In the bibliographies, they indicate that the term is of the type: A * B * (A-B). But I see that trying to adjust the three cubic terms results in singularities. I know this implies not having the inverse matrix: solve (t (X)% *% X) does not exist. The bibliographies show all three cubic terms. So my

I have two lists (and possibly more): col.prop <- structure(list(B = structure(c(0.5, 0.5, 0.6, 0.4, 0.555555555555556, 0.444444444444444), class = c("cast_matrix", "matrix"), .Dim = 2:3, .Dimnames = list( NULL, NULL)), D = structure(c(1, 0, 0.666666666666667, 0.333333333333333, 0.8, 0.2), class = c("cast_matrix", "matrix"), .Dim = 2:3, .Dimnames =

Dear R wizards: [a] I believe that the pdf device does not yet fully support TeXencoding. (under R-2.2.0, the pdf file created with Textext as font encoding still dies when post-processed by ghostscript.) are there any workarounds, or are there utilities that would allow a TeXencoded font to be re-encoded/converted into ISOLatin, perhaps, which R could then handle beautifully? [b] is there a

While building a package, I see the following: * checking R code for possible problems ... NOTE cheat.fit: no visible binding for global variable 'Zobs' plot.jml: no visible binding for global variable 'Var1' I see the issue has come up before, but I'm having a hard time discerning how solutions applied elsewhere would apply here. The entire code for both functions is below,