search for: malmig

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, Biodem has been build on a Debian Linux box and is therefore not yet available for Windows/OS X; Biodem was built with R 2.0.1 because I am using Debian 'testing' and the new R 2.1.0 is...

Hello all, To calculate the power of a matrix, I used the command "mtx.exp(X, n)", but there is an error saying "Error: couldn't find function "mtx.exp"". How can I deal with this problem? Jing

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, Biodem has been build on a Debian Linux box and is therefore not yet available for Windows/OS X; Biodem was built with R 2.0.1 because I am using Debian 'testing' and the new R 2.1.0 is...

Dear all, I have a population with three age-classes, at time t=0 the population is: n.zero <- c(1,0,0) I have a transition matrix A which denotes "fertility" and "survival": A <- matrix(c(0,1,5, 0.3,0,0, 0,0.5,0), ncol=3, byrow=TRUE) To obtain the population at t=1, I calculate: A %*% n.zero To obtain the population t=2, I calculate: A %*% (A %*% n.zero) ... and so

...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 migration model and associated functions. I am meant, according to R-exts, to upload the tar.gz at: ftp://ftp.ci.tuwien.ac.at/incoming How? never uploaded anything in my life. Cheers, Federico Calboli ================================= Federico C. F. Calboli Dipartimento...

list.files() (and dir()) don't appear to return names of directories when one uses the recursive=T argument. E.g., > dir(file.path(R.home(),"library"), pattern="^R$", recursive=T) [1] "Malmig/help/R" but the unix find commmand finds lots of R directories > z <- system(paste("find", file.path(R.home(),"library"), "-name R"), intern=T) > length(z) [1] 665 > file.info(z[1:3])[,1:3] size isdir...

..."CGIwithR","CDNmoney","clac","clim.pact","compositions","cyclones","hapassoc","haplo.score","haplo.stats","hapsim","httpRequest", "labdsv","kza","LMGene","Malmig","magic","negenes","oz","papply","spe","wavethresh","waveslim","tdthap") failPackages3 <- c("rcom","Rlsf") #put the 3 sets of rejects together failPackages <- union(failPackages1,union(...

...t X[m,1] changes with successive powers: > res_mult[1,1] [1] 0 > res_eigen[1,1] [1] 1.942362e-16 > res_mult[n,1] [1] 0 > res_eigen[n,1] [1] 1.469650e+132 Notice how poorly the off-diagonal element is approximated by the eigen-decomposition!! I am aware of the mtx.exp() function in the Malmig package, although repeatedly calling this function is not efficient as I need to extract the m-th row of the matrix raised to each power from 1:n. Any suggestions? Could I possibly exploit the fact that X is tridiagonal in my application? Thanks! Miguel -- Miguel Lacerda Department of Statistic...

Hello to everyone! I have built CRAN task view for genetics. For now I have not submit it to CRAN yet and it can be accessible from: http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.html http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.ctv I have not submitted it to CRAN, since I would like first some opinion about it. Genetics is really so broad field that I belive one person

Hi, Is there a function for raising a matrix to a power? For example if you like to compute A%*%A%*%A, is there an abbreviation similar to A^3? Atte Tenkanen > A=rbind(c(1,1),c(-1,-2)) > A [,1] [,2] [1,] 1 1 [2,] -1 -2 > A^3 [,1] [,2] [1,] 1 1 [2,] -1 -8 But: > A%*%A%*%A [,1] [,2] [1,] 1 2 [2,] -2 -5

Hi I have been wondering if there one can speed up calculating small powers of numbers such as x^8 using multiplication. In addition, one can be a bit clever and calculate x^8 using only 3 multiplies. look at this: > f1 <- function(x){x*x*x*x*x*x*x*x} > f2 <- function(x){x^8} > f3 <- function(x){x2 <- x*x;x4 <- x2*x2;return(x4*x4)} [so f1() and f2() and f3() are