similar to: elegant way to minus on each row of a matrix

Hi guys, I hit on a problem when I use rgl. Could you try to run the code here in this link and see why the first persp3d gives a red bounding box and the second shows black? http://rafb.net/p/g1i7ur33.html (sorry for not pasting the code here directly but my previous email got filtered by this list so I suspect my code looks like spam to the spam filter) I'm expecting black color to be

I used the following to assign each row the group average, where w is the new group average variable, x is the variabale to be averaged, and g is the nominal group indicator: w <- ave(x,g) Now I want to calculate the group average, but WITHOUT each row's value of x. Is there an easy way to do this? I'm sure I'm missing something obvious here, but for the life of me I can't

Hi Guys, I'm cooking up my time series code. I want a data frame with first column as timestamp in POSIXlt format. I hit on this the problem of how to create an array/list/vector of POSIXlt objects. Code is as follows > dtt=array(dim = 2) > t=as.POSIXlt( strptime("07/12/07 13:20:01", "%m/%d/%Y %H:%M:%S",tz="GMT")) > dtt [1] NA NA > t [1]

Hi Guys, Need your wisdom on this. Say I have a time series (in zoo format) like this > x <- zoo(11:21) > x 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18 19 20 21 I want to do a "moving window sampling" of it. The result can either be a matrix or a dataframe like this my.super.moving.window(x, length=3, by=1) 11 12 13 12 13 14 13 14 15 14 15 16 .... 18 19 20

Full_Name: Bo Zhou Version: 2.6.1 (2007-11-26) OS: Windows XP Submission from: (NULL) (207.237.54.242) Hi, I found an arithmetic problem when I'm doing something with POSIXct The code to reproduce it is as follows (This is the recommended way of finding out time zone difference on R News 2004-1 Page 32 URL http://cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf) a=Sys.time()

Hi Guys, I'm using zoo package now. I found lag is not doing what I assumed. > x <- zoo(11:21) > z <- zoo(1:10, yearqtr(seq(1959.25, 1961.5, by = 0.25)), frequency = 4) > x 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18 19 20 21 > lag(x) 1 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 20 21 > z 1959 Q2 1959 Q3 1959 Q4 1960 Q1 1960 Q2 1960 Q3 1960 Q4

I need to use R to model a large number of experiments (say, 1000). Each experiment involves the random selection of 5 numbers (without replacement) from a pool of numbers ranging between 1 and 30. What I need to know is what *proportion* of those experiments contains two or more numbers that are consecutive. So, for instance, an experiment that yielded the numbers 2, 28, 31, 4, 27 would be

Hello, I am new to using R and I am having problems get boot() to work properly. Here is what I am trying to do: I have statistic called "cs". cs takes a data matrix (154 x 5) and calculates 12 different scores for me. cs outputs the data as a vector (12 x 1). cs doesn't really use weights, per se, however I have included this as one of the 2 arguments cs can take. I try

Hi All, I am illustrating a simple, two-way ANOVA using the following data and I'm having difficulty in expressing the predicted values succinctly in R. X<- data.frame(read.table(textConnection(" Machine.1 Machine.2 Machine.3 53 61 51 47 55 51 46 52 49 50

Hello Suppose I have a matrix mat=(1:16,2) [,1] [,2] [,3] [,4] [1,] 1 5 9 13 [2,] 2 6 10 14 [3,] 3 7 11 15 [4,] 4 8 12 16 there is a vector end=c(2,3,1,3) #coerce the 1st 2 numbers of the 1st column to a vector [1] 1 2 #coerce the 1st 3 numbers of the 2nd column and append it to the previous vector [1] 1 2 5 6 7 #coerce the 1st number of

This is what I tried: > num.vec <- c(12.34,56.78,90.12,34.56) > names(num.vec)<-c("first","second","third","fourth") > num.vec first second third fourth 12.34 56.78 90.12 34.56 > seq<-(1:4) > num.mat<-rbind(num.vec,seq) > num.mat first second third fourth num.vec

Hi, Given a test matrix, test <- matrix(c(1,2,3,NA,2,3,NA,NA,2), 3,3) A) How to compute the counts of each column (excluding the NA) i.e., 3, 2, 1 B) How to compute the counts of each column (excluding the NA) that are greater than the column means ? i.e., 1, 1, 0 I could write a for loop, but hope to use better alternative. [[alternative HTML version deleted]]

Dear mailing list I have a data in two columns and how can i convert it to one row . thank you in advance inpute 1 2 3 4 5 6 7 8 9 1 out put 1 2 3 4 5 6 7 8 9 1 [[alternative HTML version deleted]]

Hi all, I have problems reading Unicode (UTF-16) coded tables in R 2.8.1 under Windows Vista. Imagine the following table: a b c d X 1,2 1,3 1,4 Y 2,2 2,3 2,4 Z 3,2 3,3 3,4 Usually I would use the following code to read the table: t = read.table("test.txt", header=T, sep="\t",dec=",") This works well if I create the table

I have been RTFM/doc/www, but I'm still lost. How does one tell R CMD INSTALL and R CMD check (both) that libraries are installed in non-usual places and so -I/APPS/include (or whatever) is necessary in CPPFLAGS and -L/APPS/lib (similarly) is necessary in LDFLAGS? I know I can hardwire them in pkg/src/Makevars, but this requires hand editing of that file by each installer (yuck!) rgentlem

I have been RTFM/doc/www, but I'm still lost. How does one tell R CMD INSTALL and R CMD check (both) that libraries are installed in non-usual places and so -I/APPS/include (or whatever) is necessary in CPPFLAGS and -L/APPS/lib (similarly) is necessary in LDFLAGS? I know I can hardwire them in pkg/src/Makevars, but this requires hand editing of that file by each installer (yuck!) rgentlem

Hi, I have matrix of bits and a target vector. Is there an efficient way to search the rows of the matrix for the target? I am interested in the first row index where target is found. Example: > source("lookup.R") [,1] [,2] [,3] [,4] [,5] [1,] 1 0 1 1 0 [2,] 1 1 0 1 0 [3,] 0 0 1 0 0 [4,] 1 0 0 1 1 [5,]

there are really two related problems here I have a 2D matrix A <- matrix(1:100,nrow=20,ncol =5) S <- matrix(1:10,nrow=2,ncol =5) #I want to subtract S from A. so that S would be subtracted from the first 2 rows of #A, then the next two rows and so on. #I have a the same problem with a 3D array # where I want to subtract Q for every layer (1-10) in Z # I thought I solved this one

Dear R devolepers! Wouldn't you change src/main/devPS.c file (patch will follow). Problem is with minus sign in ps output. The original code forces the name /minus to the character number 45 (-). But /minus symbol is not defined in standard Adobe encoding and (even worse) is not defined in most type1 fonts (where the name /hyphen is used instead), in this case you will have an empty

Dear R devolepers! Wouldn't you change src/main/devPS.c file (patch will follow). Problem is with minus sign in ps output. The original code forces the name /minus to the character number 45 (-). But /minus symbol is not defined in standard Adobe encoding and (even worse) is not defined in most type1 fonts (where the name /hyphen is used instead), in this case you will have an empty