similar to: acf and na.pass

Dear all, I would like to compute power of a square non symmetric matrix. This is a part of a simulation study. Matrices are quite large (e.g., 900 by 900), and contains many 0 (more than 99 %). I have try the function mtx.exp of the Biodem package: library(Biodem) m <- matrix(0, 900, 900) i <- sample(1:900, 3000, replace = T) j <- sample(1:900, 3000, replace = T) for(x in 1:3000)

Hello, I have some questions concerning encoding and package distribution. We develop the ade4 package. For some data sets included in the package, there are accentued character (e.g. ?,?...). The data sets have been saved using latin1 encoding, but some of us use utf-8 and can not see some data sets which contains accented chracters. e.g: librarry(ade4) data(rankrock) rankrock in this case,

Dear All, I would like to plot text with a box around it. I used strwidth and strheight to compute the size of the box which is plotted with rect: z <- rnorm(10) # horizontal text works plot(rnorm(10)) x1 <- 5 y1 <- 0 label <- "Label" cha <- paste(" ", label, " ", sep = "") xh <- strwidth(cha, cex = par("cex")) yh <-

Hello list, I have found a problem (bug?) with Sweave. I hope that someone could help me. Try this little example : \documentclass[a4paper]{article} \title{toto} \author{toto} \begin{document} \maketitle <<ni1, fig=T, eval=TRUE, echo=TRUE, debug=TRUE, results=verbatim, include =FALSE, width=7, height=7>>= a<-rnorm(1)+1 a plot(1:10) @ patati <<ni2, fig=F, eval=TRUE,

A L E R T E V I R U S Notre antivirus a d?tect? le VIRUS "W32/Sobig.e at MM" dans votre courrier destin? ? : -> dray at biomserv.univ-lyon1.fr Nous refusons de le transmettre. Nous vous conseillons de d?contaminer votre poste de travail. Pour votre information, vous trouverez ci-dessous les entetes de votre message Pour des

Hello list, I would like to cross-compile R packages using R 2.4.0. I am working on Linux Debian and cross-compiled (windows binaries) without problems with older R version. I have used the doc of Yan and Rossini in the contributed section of the R documentation (same version of MinGW...). When I try to cross-compile R (make R), the procedure stopped and returns :

Hello list, I am writing a paper for Rnews. I use Sweave to do it. I did not find information about writing a paper for Rnews using Sweave and have some questions: - Is there a problem to use the environment 'Sinput' in the place of 'example' or 'smallexample'. It works fine but perhaps there are some technical/editorial problems ? - I have some long lines of code in

Hi all, I've been trying to figure out how to superimpose a histogram on a biplot that shows the relative contribution of each axis. I have been using the NIPALS function (http://biomserv.univ-lyon1.fr/~dray/ files/softwares/nipals.R) to run principal component analyses. Here is a toy example. source("http://biomserv.univ-lyon1.fr/~dray/files/softwares/nipals.R")

Hello list, I would like to know if it is possible to get output of help(foo) with Sweave. If I insert ?foo in a chunk, the console give me the help but no output is obtained in the tex file. Perhaps something to do with options ? Thanks in advance ! -- St??phane DRAY (dray at biomserv.univ-lyon1.fr ) Laboratoire BBE-CNRS-UMR-5558, Univ. C. Bernard - Lyon I 43, Bd du 11 Novembre 1918, 69622

Dear all, I am trying to compile Ralpha on my computer (OS = Debian). SVN Revision is: Revision: 42843 Last Changed Date: 2007-09-14 ./configure returns : ---------------------------------------------------- stephane at pcdray:~/Rdev/R-alpha$ ./configure checking build system type... i686-pc-linux-gnulibc1 checking host system type... i686-pc-linux-gnulibc1 loading site script

Hello, i' m looking for a function to compute geary statistic for spatial correlation. Thanks. -- St?phane DRAY --------------------------------------------------------------- Biom?trie et Biologie ?volutive - Equipe "?cologie Statistique" Universite Lyon 1 - Bat 711 - 69622 Villeurbanne CEDEX - France Tel : 04 72 43 27 56 Fax : 04 78 89 27 19 04 72 43 27 57 E-mail

Hello list, (PC Windows 2000, R 1.2.3) i want to compile C source file for R on windows 2000. I'm working with Metrowerks CodeWarrior for creating and compiling C files. I've see that R.lib file must be obtain with R.ddl by Visual C++. I don't have Visual C++ (and i don't want to buy it !!) but i want to obtain R.lib file. Is there a solution to my problem ??? Any help would

Hello list, i have a glm with family=binomial, link=logit but there is over-dispersion. So, in order to take into account for this problem i choose to do a glm with family=quasibinomial(). I'm not an expert on this subject and i ask if someone could validate my approach (i'm not sure for the tests) : quasi_glm(myformula,quasibinomial(),start=mystart) summary(quasi) # test t for

Hello List, i need help for eigen and svd functions. I have a non-symmetric square matrix. These matrix is not positive (some eigenvalues are negative). I want to diagonalise these matrix. So, I use svd and eigen and i compare the results. eigen give me the "good" eigenvalues (positive and negative). I compare with another software and the results are the same. BUT, when i use svd,

Hello, I am trying to play with RODBC library and Excel Files. In my file (doubs.xls) there are 2 spreadsheets: > library(RODBC) > connection<-odbcConnect("Excel Files") > sqlTables(connection) TABLE_CAT TABLE_SCHEM TABLE_NAME TABLE_TYPE REMARKS 1 F:\\Th?se\\R\\Doubs NA Faune$ SYSTEM TABLE NA 2 F:\\Th?se\\R\\Doubs NA Milieu$

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Dear R-help, After making a typo (reminiscent of FORTRAN 77, I guess) I found the following: > identical(all.equal(2^(-10:10), 2**(-10:10)), TRUE) [1] TRUE I have tried to find the documentation about the ** operator but I was unsuccesful this way: > sessionInfo() Version 2.3.1 (2006-06-01) powerpc-apple-darwin8.6.0 attached base packages: [1] "methods" "stats"

Dear all, I'm trying to produce two dotcharts side-by-side within a Sweave document. When I'm compiling this example: \documentclass{article} \begin{document} <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) dotchart(1:10) @ <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) hist(1:10) @

Dear all, I'm confused by the time serie parameters in data(lh) : sueoka:~ lobry$ R --vanilla --quiet > tsp(lh) [1] 1 48 1 because documentation says: QUOTE A regular time series giving the luteinizing hormone in blood samples at 10 mins intervals from a human female, 48 samples. UNQUOTE So that I would expect the time serie to end at 480 minutes or 8 hours. Shouldn't we have