similar to: how to upload [forwarded]

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

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, 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)

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

> -----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 >

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

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

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

Dear all, I would like to ask you the following: I have data about different manuscripts, together with data about the presence/absence of copying errors, in the days when manuscript were really manuscripts. I would ideally like to use the data to draw a phylogenetic tree, so that I can infer which ms was copied from which. The error presence/absence is coded in binary format. The plan is to use

Dear all, I am putting together a business case for R and I would like to know if there exist statistics on the distribution of R in academia, research and commercial institutions. The same information about S-plus would be welcome. Please reply only if you have consistent and certain information. Thank you very much, giovanni

Together with Alessio Boattini of the University of Bologna we have created a package called Biodem. Biodem provides a number of functions for Biodemographycal analysis, and we hope it will be useful to the anthropological community. Because Biodem contains all the functions found in Malmig (a package I maintain), I would like to orphan it, or, even better, have it removed from CRAN. Finally,

Together with Alessio Boattini of the University of Bologna we have created a package called Biodem. Biodem provides a number of functions for Biodemographycal analysis, and we hope it will be useful to the anthropological community. Because Biodem contains all the functions found in Malmig (a package I maintain), I would like to orphan it, or, even better, have it removed from CRAN. Finally,

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

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,

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>

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 am counting the number of occurrences of the terms listed in one vector in another vector. My code runs: for( i in 1:length(vector3)){ vector3[i] = sum(1*is.element(vector2, vector1[i])) } where vector1 = vector containing the terms whose occurrences I want to count vector2 = made up of a number of repetitions of all the elements of vector1 vector3 = a vector of NAs that is