similar to: vector normal to a plane

Dear All, I have a problem that I think could be solved much more efficiently, but I don't have a clue how to accomplish this. I have a matrix W with dimensions k*(p+1): Let's say W is 5*4 and looks like this: > W [,1] [,2] [,3] [,4] [1,] 1 -1 -1 1 [2,] 1 1 1 1 [3,] 1 2 -2 -1 [4,] 1 0 -1 -1 [5,] 1 -2 -1 0 I want to take

Dear All, I have written the following programs to find a non-singular (10*10) covariance matrix. Here is the program: nitems <- 10 x <- array(rnorm(5*nitems,3,3), c(5,nitems)) sigma <- t(x)%*%x inverse <- try(solve(sigma), TRUE) while(inherits(inverse, "try-error")) { x <- array(rnorm(5*nitems,3,3), c(5,nitems)) sigma <- t(x)%*%x inverse <-

I'm trying to get R running on AIX 5.1 with the native AIX compilers, VisualAge C++ 5.0.2.0 and XL Fortran Compiler 7.1.0.0. R compiled fine, but in running the test cases, the base-Ex.R tests fail in the Choleski decomposition with the following: > x <- matrix(c(1:5, (1:5)^2), 5, 2) > m <- crossprod(x) > Q <- chol(m) Error in chol(m) : NA/NaN/Inf in foreign function call

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

I am pretty new to R. So this may be an easy question for most of you. ? I would like to calculate the squared distances of a large set (let's say 20000) of vectors (let's say dimension of 5) to a fixed vector. ? Say I have a data frame MY_VECTORS with 20000 rows and 5 columns, and one 5x1 vector y. I would like to efficiently calculate the squared distances?between each of the 20000

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

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 guys, I am still learning R, and not well familiar with all the apply functions. I am trying to find faster alternatives to replace the for cycle. Is there a faster way to do the example below? nm <- 1000 b <- matrix (rnorm (5000, 0, 1), nrow = 500, ncol = nm) a <- matrix (0, nm, nm) for (i in 1 : nm) { for (j in 1 : nm) { if ( j == i) { next } a[i, j] <- t (b [, i]) %*% b[, j] } }

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

Dear R helpers, I have a 3D scatter plot that I have generated from scatterplot3d (which looks great- thanks!) and I can see that the points in my graph fall in a plane. Following the example 5 from 3D scatter plot (below) I have fitted a regression plane. Now what I would like to do is a rotation so that my new co-ordinate system is about the fitted plane (by finding the normal to the plane

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"