similar to: R-business case

> -----Original Message----- > From: Federico Calboli > Sent: Tuesday, January 20, 2004 5:40 PM > To: r-help > Subject: [R] matrix exponential: M^0 > > I would like to ask why the zeroeth power of a matrix gives me a matrix > of ones rather than the identity matrix: > > > D<-rbind(c(0,0,0),c(0,0,0),c(0,0,0)) > > D<-as.matrix(D) > > D >

Check out textplot in the gregmisc package ... Federico Calboli <f.calboli@ucl.ac.uk> Sent by: r-help-bounces@stat.math.ethz.ch 07/13/2004 12:06 PM Please respond to f.calboli To: r-help <r-help@stat.math.ethz.ch> cc: Subject: [R] plotting a table together with graphs Dear All, I would like to ask how to add a table to a

[ This was posted to the R-packages list (which I moderate), but definitely doesn't belong there. Martin Maechler ] ------- start of forwarded message ------- From: Federico Calboli <f.calboli at ucl.ac.uk> To: r-packages at stat.math.ethz.ch Subject: how to upload Date: 15 Jul 2004 18:30:26 +0100 Dear All, I just finished a add on lib called Malmig, to calculate Malecot

Dear All, I would like to ask how to get the Sum of Squares from fitted lme model. I appreciate that lme maximises the likelihood (or REML) and uses likelihood ratio tests, but I just fail why I could not get the SS if I want them. I could use lm to calculate them, but it looks quite pointless to fit a second (and wrong in a way) model for something so trivial. I know that the issue has been

Dear All, I apologize if the following has been already asked, but I could not find anything in the archives. I am running Debian Sid (on an Athlon XP), with R 1.8.1.cvs.20040307-1 (it describes itself as "R 1.9.0 experimental" when I fire R up) and I am unable to load the library "nlme" (I apt-getted r-recommended, so the library is on the system) >library(nlme) Error in

Dear all, We compared the gene expression level in two conditions: high (H) and low (L), using real-time PCR. We had 3 samples of each condition, and we quantified each sample twice (i.e. a technical replication). Our data looks like: H1 H1 H2 H2 H3 H3 L1 L1 L2 L2 L3 L3 We want to test if the level in H vs L is significantly different. I would like to use ANOVA and I guess a nested model with

Dear All, how do I compute the left eigenvector of a matrix? I gather that "eigen" computes the right eigenvectors... Regards, Federico Calboli -- ================================= Federico C. F. Calboli PLEASE NOTE NEW ADDRESS Dipartimento di Biologia Via Selmi 3 40126 Bologna Italy tel (+39) 051 209 4187 fax (+39) 051 251 208 f.calboli at ucl.ac.uk

Dear All, I would like to ask how to customize the graph corresponding to a procrustes analysis. I have to distance matrices, that I transform to two set of coordinates by means of muti dimensional scaling: library(mva) c1<-cmdscale(mat.dist1) c2<-cmdscale(mat.dist2) I vant to rotate c2 on c1, and I use the "procrustes" analysis from the {vegan} library. library(vegan)

I use linux debian testing for which the latest debianized version of R is 1.8.1. Therefore I installed: r-base-core_1.8.1-0.cran.1_i386.deb r-base-dev_1.8.1-0.cran.1_all.deb r-base-html_1.8.1-0.cran.1_all.deb r-base-latex_1.8.1-0.cran.1_all.deb r-base_1.8.1-0.cran.1_all.deb r-doc-html_1.8.1-0.cran.1_all.deb r-doc-info_1.8.1-0.cran.1_all.deb

I am using a two-factor ANOVA model with random factor effects including the interaction, i.e. the factors are crossed. I would like to be able to generate all four variance components along with approximate confidence intervals using the NLME package. However, I do not know how to specify the random option because of two problems. First, I do not know how to enter the interaction term into the

Dear All, I am in the somewhat unfortunate position of having to reproduce the results previously obtained from (non-metric?) MDS on a "kinship" matrix using Statistica. A kinship matrix measures affinity between groups, and has its maximum values on the diagonal. Apparently, starting with a nxn kinship matrix, all it was needed to do was to feed it to Statistica flagging that the

Hello everybody, I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled with gcc 3.2.2. The p-value calculated from the chisq.test function is incorrect for some input values: > chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE) Pearson's Chi-squared test with simulated p-value (based on 2000 replicates) data: matrix(c(0, 1, 1,

Hi All, I have the following adjacency matrix for a directed graph: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0 0 0 0 0 0 0 0 [2,] 0 0 0 0 0 0 0 0 [3,] 1 0 0 0 0 0 0 0 [4,] 0 0 1 0 0 0 0 0 [5,] 0 0 1 0 0 0 0 0 [6,] 1 1 0 0 0 0 0 0 [7,] 0 0

How about this: > expm function(x,pow=2) { xd <- diag((eigen(x)$values^pow)) xm <- eigen(x)$vector %*% xd %*% t(eigen(x)$vector) return(xm) } > xa [,1] [,2] [1,] 2 1 [2,] 1 3 > expm(xa,pow=3) [,1] [,2] [1,] 15 20 [2,] 20 35 > xa %*% xa %*% xa [,1] [,2] [1,] 15 20 [2,] 20 35 > From: Federico Calboli <f.calboli at ucl.ac.uk>

Dear All, I would like to ask if there is any package that would calculate Malecot migration matrix (or the Imaizumi or Pena implementations of the MMM). I did a site-wide search, and went trough the table of contents of a number of packages that seemed promising, but with no luck. Regards, Federico Calboli -- ================================= Federico C. F. Calboli PLEASE NOTE NEW

Dear all, Actually, I want to create my own test-bed based on xen hypervisor with centos prefered. However, all I have is one desktop (with a good hardware performance). So it possible to have a dual-boot (a linux for daily use and the other for the test-bed) ? Regards -- Jihed M’SELMI http://about.me/jihed.mselmi _______________________________________________ Xen-users mailing list

Hi All, just to check, is cran2deb available for Karmic Koala (ideally also for a amd_64)? Cheers, Fede -- Federico C. F. Calboli Department of Epidemiology and Public Health Imperial College, St. Mary's Campus Norfolk Place, London W2 1PG Tel +44 (0)20 75941602 Fax +44 (0)20 75943193 f.calboli [.a.t] imperial.ac.uk f.calboli [.a.t] gmail.com

Hi All, I have a list of matrices: > x [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 > y [,1] [,2] [,3] [,4] [,5] [,6] [1,] 18 21 24 27 30 33 [2,] 19 22 25 28 31 34 [3,] 20 23 26 29 32 35 > z =list(x,y) I want to create a second list that is has a subset each matrix in the list subsetting so I get the 2nd and 3rd row of each (and

Dear All, I saved by mistake the environment I was working in after typing q(), and now I get the annoying message: [Previously saved workspace restored] I have already deleted all the objects in the environment, saving it as an empty environment, so it's just a matter of nitpicking I suppose. The message does not appear if I start R from any other place in the directory tree. I am

HI All, I have written and succesfully uploaded a new package. The licence it is under is 'GPL' --no version. My assumption is, since all the code is written in R the licence R used for R would affect the code (hence my "GPL" stands for "whatever version of the GPL R is under") I am happy with the licencing I used, but I'd like to ask if there is any transitive