similar to: Sparse Matrices in R

Hi, I am writing some basic smoothers in R for cleaning some spectral data. I wanted to see if I could get close to matlab for speed, so I was trying to compare SparseM with Matrix to see which could do the choleski decomposition the fastest. Here is the function using SparseM difsm <- function(y, lambda, d){ # Smoothing with a finite difference penalty # y: signal to be smoothed #

Full_Name: Julien Gagneur Version: 2.2.1 OS: Linux (Suse 9.3) Submission from: (NULL) (194.94.44.4) Using the SparseM library (SparseM_0.66) and the e1071 library (e1071_1.5-12) I fail using svm method with a sparse matrix. Here is a sample example. I experienced the same problem under Windows. > library(SparseM) [1] "SparseM library loaded" > library("e1071")

Art (and group), I'm doing this as a form of missing value analysis. Approximately 30% of the cases are missing data for one variable. To impute values for those cases, I'd like to match those cases that are missing the variable to all other cases and then take an average of those to infill. I realize there are many methods for imputing data. I'm not well versed on any in

Hi, hres <- hclust(smatr,method="single") hresd<-as.dendrogram(hres) as.dendrogram(hres) `dendrogram' with 2 branches and 380 members total, at height 2514.513 plot(hresd,leaflab="none") #<-error here. #the plotted dendrogram is incomplete. The x axis is not drawn. #The interested reader can download the save(hresd,file="hres.rda") #from the

Hi! I am installing R in non standard directory. ./configure --prefix=/non/standard/directory The installation works fine. But after starting R i get .libPaths() /usr/lib/R/library but they cant be there of course! Have I to set some additional switches during config? Eryk Dipl. bio-chem. Eryk Witold Wolski @ MPI-Moleculare Genetic Ihnestrasse 63-73 14195 Berlin 'v' tel:

hello I am returning to some libraries that I had previously customised by amongst other things adding additional functions. I had simply typed these new functions into the file in library/R/thelibrary. These however do not seem to be loaded now as they previously were under older versions of R. What's changed and what do i need to do? Thanks Stephen

Hi! Would like to initialize a huge matrix.csc (Pacakge SparseM) with all elements 0 and afterwards set a few alements nonzero. The matrix which I like to allocate is so huge that I can not use A <- matrix(a,n1,p) before: A.csr <- as.matrix.csc(A) because I can not allocate such a huge matrix A. But I believe that the much more memmory efficient model in case of csc matrix should do it for

This seems like a bug to me. Can someone verify this? First we define a function f that returns its second argument and lapply it to 1:2 using 9 as the second argument and all seems well. Note that "{" as a function does the same thing as f, as illustrated with f(1,9) and "{"(1,9); however, when we attempt to use "{" in the very same way we used f in lapply,

Hi! I encountered some things trying to build own packages under Windows 2000 for which i found no answer in the documentation. First. 1. I do not found the rcmd.exe in the bin directory of R_HOME. I installed R using the setup.exe 2. I updated the R intatllation with one of the zip files. (rw1041sp.zip) now i had the rcmd.exe in the bin. 3. I do not found any hint in the docu that perl is

Hi! I guess that it cant work. but maybee I am wrong. I would like to define a function rbind. > rbind function (..., deparse.level = 1) { which works only for objects of class Myclass Is it possible to use something like setMethod("rbind",signature(...="Myclass",deparse.level="numeric") This gives an error. Or should I use rbind.Myclass Eryk. Dipl. bio-chem.

Hi! I would like to give my objects functions like "subset", "union", et cetera and name it also in this way. But this functions names are already used in the base package and they are not generic. I am right in that if they would be generic then i still can use the neat names with my objects? I have the impression that there are more and more packages are going to be

Hello, I'm working with a very large, very sparse X matrix. Let csr.X <- * as.matrix.csr*(X) as described by the SparseM package. The documentation says that "Indexing .... work just like they do on dense matrices". To me this says that I should be able to perform operations on the rows of csr.X in the same way I would on X itself. E.g. f <- function(x){ for (i in 1:n){

Hello All: This function obviously fails x <- function(z) paste("go", z, sep = ".") <- 10 x("now") But is there a way to define the name of a variable through passing a parameter in a function call? Thanks, ANDREW

Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document figure the quality is bad. Anyone knows a way to reduce the size but keep the quality? /E -- Dipl. bio-chem. Witold Eryk Wolski MPI-Moleculare

Hi! I have installed the new R on windows. I wanted to reinstall the XML package. I am not able to find the XML.zip anymore. I am quite shure that they where a windows binary version. Has anyone old XML windows binary? Eryk Dipl. bio-chem. Eryk Witold Wolski @ MPI-MG Dep. Vertebrate Genomics Ihnestrasse 73 14195 Berlin 'v' tel: 0049-30-84131285 / \ mail:

Hi! Would like to plot two image plots with color scale (like levelplot (lattice) implements) in one graphic device. Tried to plot 2 levelplot's (lattice) (the data has not much in common) in one graphics device by setting par(mfrwo=c(2,1)) which does not work. Is there a way to force levelplot to use only a half of the graphic device? Is there a function like image which provides the color

Hi! >library(cluster) In this lib a function called silhoutte.default is defined than on the R prompt it type >silhouette.default Error: Object "silhouette.default" not found R1.9.0 The same error are at R1.8.1 And I knew that a function silhoutte.default are present. But the same piece of code works in R1.6.2 ??? Eryk Dipl. bio-chem. Eryk Witold Wolski @

I have been struggling with this problem for a while and I hope someone could help me. Or if someone could point me to a section in the manual I would be grateful. x <- "my" plot(1:10, main=expression(paste( x, Delta, "values" ))) Q : How do I get the title to say "my (triangle symbol) values" ? The following trial-and-error produced mainly errors :

Hi, Executing the following code causes the prompt to enter browser() mode. >dd<-as(Npatrick1,"list") Called from: as(Npatrick1, "list") Browse[1]> R : Copyright 2004, The R Foundation for Statistical Computing Version 2.0.0 beta (2004-09-20), ISBN 3-900051-07-0 Some addtitional info about Npatrick1 > is(Npatrick1) [1] "Massvectorlist"

Hi! How to generate a empty string of a precise length without using a loop? I need a empty string (char **) to pass it with .C function to a c dll. The c routine writes the result into the string. Eryk. Dipl. bio-chem. Eryk Witold Wolski @ MPI-MG Dep. Vertebrate Genomics Ihnestrasse 73 14195 Berlin 'v' tel: 0049-30-84131285 / \ mail: