similar to: xyplot of function only shows diagonal plots

Hi r-help community, This question isn't so much a syntax/coding one, but here goes: Let's say I have matrix of arbitrary dimensions and I'd like to reorder the rows in such a way that I could maximize the sum of the entries along the diagonal. For example, for this 3x3 matrix: [,1] [,2] [,3] [1,] 3 4 13 [2,] 9 1 2 [3,] 2 11 1 rearranging the rows

Dear list, Sorry for asking this question, but has something changed in the syntax for bwplot in Lattice? In an old publication, I used > bwplot( VOTMS ~gender |type * group, data=merge(vot,words,by="ord"), nint=30, horizontal=F, layout=c(3,3), box.ratio=0.8) which produced a lovelly 3x3 lattice plot with one

Hi, I am trying to find an efficient way of applying a linear regression model to different factor combinations in a data frame. I want to obtain the output with minimal or no use of loops if possible. Please let me know if this query is unclear. Thanks, Murtaza

Hi All, I have a greyscale image that I am reading in through RGDAL and placing in a pixmap object. As an example use the logo.jpg file that comes with the RGDAL package: #Read the file logo <- system.file("pictures/logo.jpg", package="rgdal")[1] x <- new("GDALReadOnlyDataset", logo) #Create the pixmap object xGrey <- getPixmapGDAL(x) Now I would

Hi, On 03-08-15 20:09, Ilia Mirkin wrote: > On Mon, Aug 3, 2015 at 1:31 PM, Hans de Goede <hdegoede at redhat.com> wrote: >> Hi, >> >> >> On 03-08-15 17:36, Ilia Mirkin wrote: >>> >>> On Mon, Aug 3, 2015 at 9:02 AM, Hans de Goede <hdegoede at redhat.com> wrote: >>>> >>>> Hi, >>>> >>>> On

Hi all, Why cov(y, y) only gives one value, and cov(t(y), t(y)) gives 3x3 NA matrix? Here my y is listed below and it is a 3x1 matrix. I am expecting that if I have a random vector y=[y1 y2 y3]', here " ' " denotes a transposition so that y is a column vector, where y1, y2, y3 are independent random variables... then cov(y, y) should be E[ y * y' ] - E[y] * E[y] ',

I want to form a 3x3 crosstabulation for the signs of two vectors (i.e. Negative, Zero, Positive). The problem is that I am simulating the data so for some iterations one of the categories is absent. Thus the resulting table shrinks to 3x2. I want it to be 3x3 with zero column corresponding to the missing category. Moreover, I have tried but failed to give the dimension names. -- Sohail Chand

useR's I have a matrix from which I want to take multiple subsets from, according to a particular scheme I will now describe. The matrix below (mat) is 5x5, and I want to take 9 subsets of it, each of dimension 3x3. The best way to explain what the result should look like is with the following: > dat <- c(3,6,1,9,12,9,2,10,6,5,3,13,1,4,8,9,4,6,10,11,2,7,3,5,10) > miss <-

Hello, I am calculating this thing with vectors (b) and matrices (G,P): b'G/sqrt(b'Pb) where the denominator is a quadratic form and therefore always a scalar. In Scilab, it is quite simple: b'*G/sqrt(b'*P*b) However, in R, the denominator is an (1x1)matrix and R claims it is non conformable and I have to use drop() or as.numeric(). Like this: > b = 1:2 > G=diag(1,2) >

Dear List, I have a couple of (~200) 3x3 transition probability matrices (ie each defines a Markov chain). They are all estimated from the same underlying process, so it ie meaningful to take their elemetwise mean and standard deviation. [1] First question: supposing that they are given in a list l, how do I get their elementwise mean and standard deviation? Fortunately, the mean of trans. prob.

Dear Friends, I have been trying to save multiple 3x3 (mfrow=c(3,3) graphics inside a loop using tiff figure format (not using PDF or savePlot functions) with no success. Could you please help? Here is a simplified example code: dat=data.frame (ID=rep(1:10,each=10),IDV=rep(seq(1:10),times=10)) dat$DV <- with(dat, 50+15*IDV) dat=dat[order(dat$ID,dat$IDV),] for(i in 1:10){ dt1 =

Hello R Users, I am new to R. I have two data sets i) TOMS aerosol optical depth(AOD) and ii) TOMS ozone(O3). > > AOD data is on 1x1 grid and O3 data is on 5x5 grid. > > First I want to read AOD and O3 as it is and then I want to regrid AOD on > 5x5 grid as O3. > > Reading is first problem. > > FIRST PROBLEM READING AOD: > > AOD data is in following format: >

Hi, In function chol2inv with the option LINPACK set to false (default), it raises an error when the matrix is 1x1 matrix (i.e. just a real) saying 'a' must be a numeric matrix This error is raised by the underlying C function (modLa_chol2inv in function Lapack.c). Everything is normal, but I wonder if we could have another behavior when we pass a 1x1 matrix. I spent time this morning

Dear colleagues, I would need to add text (some rows of information) in a matrix. For example, given this matrix 1 2 3 4 5 6 7 8 9 I would need to add this info: THIS IS AN EXAMPLE OF a 3x3 MATRIX 1 2 3 4 5 6 7 8 9 I have been looking for a function that works similar to "fopen" in matlab, but unfortunately I have not found It in R. Thank you in advance for your help! Carlos

R-help, I would like to place nine (3X3) plots per page. I am not properly implement mfrow(3,3) in the script below: jpeg("xyplot.jpg") #names output file my_args <- commandArgs() #sets up to take args from dos batch command mfcol(3,3) #set page for 3X3 TEMPS <- c(-15,25,85)#list of temps VBATS <- c(3,3.6,4.7)#list of Bats BOARDS <-

Hi All, My apologies if this is a totally newbie question. I want to calculate the probability that a particular set of genes is picked repeatedly for 3 samplings. For example, if from a total of 'N' genes, I pick 'A' number of genes in the first pick, 'B' number of genes in the second pick, and 'C' number of genes in the third pick, I would want to calculate

Hi guys, I used the following R code to generate one plot library(lattice) xyplot(Y~X1|as.factor(X2)*as.factor(X3), groups = as.factor(X4), data=mydata) Both X2 and X3 have three values. X4 has two values. I got 3x3 grids and in each grid there are two curves about y~x1 for the two X4 values. I am quite happy with the plot except that I need a different layout of the 3x3 layout. For example,

Hi all, I try to make a program for comparison between matrix cells. Here um simple example: 3 |a b a 2 |b b b 1 |a b a ------- 1 2 3 This is tabulated like this: x y sp 1 1 a 1 2 b 1 3 a 2 1 b 2 2 b 2 3 b 3 1 a 3 2 b 3 3 a Then, I need to make the follow calcule: ============================================== Change in sp levels from 1x1 area to 2x2 area:

Hi all, I try to make a program for comparison between matrix cells. Here um simple example: 3 |a b a 2 |b b b 1 |a b a ------- 1 2 3 This is tabulated like this: x y sp 1 1 a 1 2 b 1 3 a 2 1 b 2 2 b 2 3 b 3 1 a 3 2 b 3 3 a Then, I need to make the follow calcule: ============================================== Change in sp levels from 1x1 area to 2x2 area:

Hello, I'm having a bit of trouble connecting our current CentOS Openswan server with a Vyos server via IPSec. I've posted this on the VyOS forums, but haven't had many helpful responses, so I thought I would ask here. http://forum.vyos.net/showthread.php?tid=26504&pid=29703#pid29703 Basically our Openswan configuration is as follows: conn VYOS keyingtries=0