search for: rowvar

I am using the R 2.2.1 in a Windows XP environment. I have a dataframe with 12 columns and 1,000 rows. (Some of the rows have 1 or fewer values.) I am trying to use rowVars to calculate the variance of each row. I am getting the following message: ?Error in na.remove.default(x) : length of 'dimnames' [1] not equal to array extent? Is there a good work-around?

I've posted a sample file for estimating loglinear models for square tables (mobility models) at http://www.xs4all.nl/~jhckx/mcl/R/ Comments and suggestions are welcome. John Hendrickx

Mark Leeds <mleeds at mlp.com> wrote (to S-News): > does anyone know the R equivalent of the SPlus rowVars function ? Andy Liaw <andy_liaw at merck.com> replied: > More seriously, I seem to recall David Brahms at one time had created an R > package with these dimensional summary statistics, using C code. (And I > pointed him to the `two-pass' algorithm for variance.) Here are the...

Hi, I have what I think is an easy question. I have a data frame, called stockdata, of stock prices that looks like this: date ticker close 1 01/02/1998 GE 24.667 2 01/05/1998 GE 25.104 3 01/06/1998 GE 24.771 4 01/07/1998 GE 24.979 5 01/08/1998 GE 24.750 6 01/02/1998 HIT 71.125 7 01/05/1998 HIT 72.313

I need to calculate a large number of t statistics, and would like to do so via matrix operations. So far I have figured out a way to calculate the mean of each row of the matrix: d <- matrix(runif(100000,1,10), 1000, 10) # some test data s <- rep(1,ncol(d)) # a sum vector to use for matrix multiplication means <- (d%*%s)/ncol(d) This is at least 1 order of magnitude faster than

Dear R friends I wanted a function to make a simple percent table that would be easy for students to use. The goal originally was to have a simple thing people would call like this pctable(rowvar, colvar, data) and the things "rowvar" and "colvar" might be names of variables in data. I wanted to avoid the usage of "with" (as we now see in the table help). Then some people wanted more features, and I agreed with the suggestion to create a formula interface tha...

...test stable REvolution R on OS X if that is pertinent. Here's the contents of a simple test2.csv CSV file: #,Status,Project 5842,New,Test > snortalerts = read.table( "/Users/lcox/Documents/test2.csv", header=TRUE, > sep=",", row.names="#") Error in data[[rowvar]] : attempt to select less than one element I can't see how it could get any more simple, yet it doesn't work. I'm obviously missing something basic, but based on the error, I can't see what it is. Landon -- View this message in context: http://www.nabble.com/Simple-question-abo...

Yesterday I asked for help about read.table with a CSV file. I received the following help. As always many thanks for the prompt responses. Now I load my csv file in a text editor ( pfe) and delete all spaces. My original mail is at the end. From: Peter Dalgaard BSA <p.dalgaard at biostat.ku.dk> Do you have spaces before the commas in your file? > Yes, number , space,

Yesterday I asked for help about read.table with a CSV file. I received the following help. As always many thanks for the prompt responses. Now I load my csv file in a text editor ( pfe) and delete all spaces. My original mail is at the end. From: Peter Dalgaard BSA <p.dalgaard at biostat.ku.dk> Do you have spaces before the commas in your file? > Yes, number , space,

I've uploaded a package colSums_1.0.tar.gz to CRAN /src/contrib/Devel. It contains functions colSums, colMeans, colVars, colStdevs, rowSums, rowMeans, rowVars, and rowStdevs. These do simple, fast arithmetic on columns/rows of a matrix, or more generally across dimensions of an array, e.g. colSums(m) = apply(m, 2, sum) but faster. They should be compatible with the corresponding S-Plus functions. The core code was written by Doug Bates <bates at s...

Hi people, If there is a fn "colMeans" why isn't there a "colSDs" etc Thanks, Phil. -- Philip Rhoades Pricom Pty Limited (ACN 003 252 275) GPO Box 3411 Sydney NSW 2001 Australia Mobile: +61:0411-185-652 Fax: +61:2:8923-5363 E-mail: pri at chu.com.au -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all, I am trying to run a two sample t-test on a matrix which is a 196002*22 matrix. I want to run the t-test, row-wise, with the first 11 columns being a part of the first group and columns 12-22 being a part of the second group. I tried running something like (temp.matrix being my 196002*22 matrix) t.test(temp.matrix[,1:11],temp.matrix[,12:22],paired=TRUE) or somthing like

...is summary(svytable(~dietperception+dietstatus,dudleyls1rake,na.rm=TRUE),"C hisq") It is throwing up an error message as follows: Error in NCOL(y) : only 0's may be mixed with negative subscripts In addition: Warning messages: 1: is.na() applied to non-(list or vector) in: is.na(rowvar) 2: is.na() applied to non-(list or vector) in: is.na(colvar) The dietperception data set does have some NA's where as there are none in dietstatus. The table is svytable(~dietperception+dietstatus,dudleyls1rake) dietstatus dietperception Good...

...lease. Of the functions in R packages that are loaded by default, hist.default is the only one I found with a "freq" argument. -------------------------------------------------- These functions have both "weights" and "freq" arguments in S-PLUS: hist rowMeans rowSums rowVars colMeans colSums colVars colQuantiles cor, var, stdev density skewness, kurtosis qqline qqnorm.default quantile Some others have "weights" but not (yet) "freq", including mean median ppoints tabulate. Other functions like lm have always had a weights argument. --------------...

Hi, Currently, the "base" has colSums, colMeans. It seems that it would be useful to extend this to also include colMin, colMax (of course, rowMin and rowMax, as well) in order to facilitate faster computations for large vectors (compared to using apply). Has this been considered before? Please forgive me if this has already been discussed before. Thanks, Ravi Ravi Varadhan, Ph.D.

Hi together, I have some data in a matrix structure - say 1000 rows with 10 columns. And I like to do some calculations (like max, avg or min) on row level. The only solution I found so fare was using a loop like for (i in 1:1000) { max[i] <- max(matrix[I,]) } It looks like that this is not very fast. Does an other way exists? Kind regards Ulrich