similar to: My question on creating binaries

How do i unsubscribe??

Hi I am running a Savitzky-Golay smoothing filter (http://tolstoy.newcastle.edu.au/R/help/04/02/0385.html) for variables in my dataset, dim (272:90). I managed to run the code for individual variables in the dataset and then combine the results into a single dataset. My novice attempt at this task is shown below csg<-NULL for (i in 1:ncol(data.all)) {

Issue: Three WinXP Pro boxes using Offline Files disconnect from all Samba shares when a user either performs a synch or accesses a file on the share with offline files enabled. About half the time, being able to reconnect for any length of time requires a reboot of the WinXP machines. Disabling the Offline Files in Windows alleviates the problem. Background and Network

Hi I'm seeing some "odd" behaviour with cbind(). My code is: > cat <- read.table("cogs_category.txt", sep="\t", header=TRUE, quote=NULL, colClasses="character") > colnames(cat) [1] "Code" "Description" > is.factor(cat$Code) [1] FALSE > is.factor(cat$Description) [1] FALSE > is.factor(rainbow(nrow(cat))) [1]

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

I have some data coming from SQL sources that I wish to relate in various ways. For reasons only known to our IT people, this can't be done in SQL at present. I am looking for an R'ish technique for looking up new columns on a data frame. As a simple, hardwired example I have tried the following: # This gives me two columns, one the lookup value and the second one # the result column,

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 List, I would like to orthonormalize vectors contained in a matrix X taking into account row weights (matrix diagonal D). ie, I want to obtain Z=XA with t(Z)%*%D%*%Z=diag(1) I can do the Gram-Schmidt orthogonalization with subsequent weighted regressions. I know that in the case of uniform weights, qr can do the trick. I wonder if there is a way to do it in the case of non uniform

Hello, I would like to have some information about acf with missing values. Let us consider this example: x=rnorm(100) x2=x x2[sample(100,10)]=NA acf1=acf(x) acf2=acf(x2,na.action=na.pass) The computation of the acf is different for the two data sets. Looking at the the code reveals that with missing values, the computation of acf takes into account the number of pairs of non-NA values (i.e.

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

Hello list, Sorry, these questions are not directly linked to R. If I consider an indefinte real matrix, I would like to know if the symmetry of the matrix is sufficient to say that their eigenvectors are real ? And what is the conditions to ensure that eigenvectors are real in the case of an asymmetric matrix (if some conditions exist)? Thanks in Advance, St?phane DRAY

Greetings. I am working on setting up an asterisk server (SIP only) and am running into a few issues getting RTP working correctly. Here is our setup: SIP Client (Public IP) <---> Asterisk Server (Public IP/Private IP) <--> Nortel CSG (Internal IP) <--> PSTN So far we have SIP to SIP working through Asterisk without any problems (using various sip clients). When I call from

Hello list, I have a problem ... and do not know how to solve it. I would like create a function that estimate the quality of the fit of different functions with different values of parameters. This problem is related to the following one: I would like to create a function "evaluatemyfunction" which evaluate a function f for different values of parameters. For instance, if

Hi Thanks once again for your help, I do appreciate it..... however.... Here is what I get with your test.... (under tcsh - i normally use bash, but I will keep everything the same) users/mwatson> env http_proxy=http://wwwcache.bbsrc.ac.uk:8080/ R >options(internet.info=0) >update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' unable to connect to

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

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 used scan function to import data into R. I have done some analysis and find strange results. I have found my problem : when importing data with scan, this can slightly modify the data : > write(c(0.251,3.399,-0.481,0.266),"essai.txt") > scan("essai.txt") Read 4 items [1] 0.251 3.399 -0.481 0.266 > print(scan("essai.txt"),17) Read

Hello list, my question is related to svd of a matrix: b=matrix(rnorm(50),10,5) mysvd=svd(b) I would like to compute each xi where xi = di* ui %*% t(vi). I do it by : xlist=sapply(1:ncol(b), function(x1,y) y$d[x1]*y$u[,x1]%*%t(y$v[,x1]),y=mysvd,simplify=F) # result is a list xarray=array(sapply(1:ncol(b), function(x1,y) y$d[x1]*y$u[,x1]%*%t(y$v[,x1]),y=mysvd),c(nrow(b),ncol(b),ncol(b))) #

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,