search for: tblackw

I am calling plot() with argument pch as a vector of numeric symbol codes, the same length as x and y. Is there some code which produces no symbol - a blank - so that I can come back with a second call to points() and fill in these locations using a different fill color and a different symbol size ? There's always a work-around, but both x and y are generated on the fly by calls

Hi dear R-listers, I'm trying to fit a 3-level model using lme in R. My sample size is about 2965 and 3 factors: year (5 levels), ssize (4 levels), condition (2 levels). When I issue the following command: > lme(var~year*ssize*condition,random=~ssize+condition|subject,data=smp,method ="ML") I got the following error: Error in logLik.lmeStructInt(lmeSt, lmePars) :

...haracter rather than factor and still reads in the numeric columns as numbers. You can also specify specific columns such as as.is=3:4 if you want columns 3 and 4 to be character but other character columns to be factors. --- Date: Mon, 15 Mar 2004 11:04:08 -0500 (EST) From: Tom Blackwell <tblackw at umich.edu> To: Janet Gannon <jangann at mote.org> Cc: <r-help at stat.math.ethz.ch> Subject: Re: [R] Simple numeric "as.is" question Janet - Try x2 <- as.numeric(as.character(x)) hist(x2) I'm not a Windows user, so I can't test this before sendin...

Hi R maintainers: when I use update.packages() And I try to update the package "maps" maps : Version 2.0-9 in D:/rw1080/library Version 2.0-10 on CRAN Update (y/N)? y The following message appears: trying URL `http://cran.r-project.org/bin/windows/contrib/1.8/maps_2.0-10.zip' Error in download.file(url, destfile, method, mode = "wb") : cannot open URL

...or 'outer' l.w <- length(w) l.x <- length(x) z <- abs( rep( x , l.w ) - rep( w, each = l.x ) ) dim( z ) <- c( l.x, l.w ) 3) find s as follows: s <- pmax( abs( w - min(x) ), abs( w - max(x) ) ) 4) use colSums() instead of apply(, , sum) Thomas W Blackwell <tblackw at umich.edu>: Matrix z is potentially enormous. On large matrices, it's a familiar truism that as.vector(z %*% rep(1, dim(z)[2])) is MUCH faster than apply(z, 1, sum). The equivalent "inline" code for "max" is not so nice. So here's an alternative version...

I have used PCA for data classification by visual examination of the 3D scatter plot of the first 3 principal components. I now want to use the results to predict the class for new data. I have used predict.princomp to predict the scores and then visualise the results on a 3D scatter plot using the rgl library. However, is there an R function that will fit the new data to the class assignments

I have checked that section already. Sorry, I should have mentioned that. Memory limit increase does not work. Installtion of msvcrt.dll does not work either. Thank you. -----Original Message----- From: ripley at stats.ox.ac.uk [mailto:ripley at stats.ox.ac.uk] Sent: Wednesday, March 05, 2003 2:44 PM To: Igor Oleinik Cc: r-help at stat.math.ethz.ch Subject: Re: [R] problem with cclust[er]

I have a factor (with "n" observations and "k" levels), with only "nobs" < n of the observations not missing. I would like to produce a (n x k) model matrix with treatment contrasts for this factor, with rows of NAs placeholding the missing observations. If I use model.matrix() I get back a (nobs x k) matrix. Is there an easy way to get the (n x k) without

I have the 'breaks' for the histogram ('hist') but I want weight the cells instead of using actual observations. I thought that using freq=FALSE implied that the numbers in 'x' were weights but this turned out to be wrong. Any help and/or comment is very much appreciated. Regards, M?rten M?rten Bjellerup Doctoral Student in Economics School of Management and Economics

Greetings, I am trying the get the standard errors of multiple comparisons using Tukey's HSD. These are not reported by the function TukeyHSD. When I apply the following code to the data, which I store as PROLE4.TXT, several unexpected things happen. First, the function TukeyHSD works for all the comparisons but the function simint doesn't. Second, after the application of na.omit

Dear r-help - I just noticed that in my R-1.7.1 on i386-pc-linux-gnu, the page displayed by help("print") ends with the line " ## Printing of factors illustrated for ex" and then no more. It looks as though something got truncated here. I think this is an R that I compiled from source off of CRAN, but I can't quite remember. - tom blackwell - u michigan

Dear list - I just discovered to my surprise that I cannot define a function with an argument named 'break' or 'while'! 'breaks' is okay. Maybe this is no surprise to the R developers. R-1.7.1, 2003-06-16, i686-pc-linux-gnu. - tom blackwell - u michigan medical school - ann arbor -