search for: pavlicova

...ersion 1.7.0 produces very different picture which I believe is wrong (compare files version161.with_at.jpg and version170.with_at.jpg). WHY is that? Am I missing something? I have attached all 4 plots and also the commands I used to create the small example. Thank you for all your help. Martina Pavlicova --------------------------------------------------------------------- library(lattice) x <- row(matrix(NA,11,11))-6 y <- col(matrix(NA,11,11))-6 z <- x*y jpeg("version161.without_at.jpg") foo <- levelplot(z~x*y,contour=T, cuts=5, ##at=c(-10,-5,0,5,10),...

without having to check the vector element by element? Thanks a lot! Jason ===== Jason G. Liao, Ph.D. Division of Biometrics University of Medicine and Dentistry of New Jersey 335 George Street, Suite 2200 New Brunswick, NJ 08903-2688 phone (732) 235-8611, fax (732) 235-9777 http://www.geocities.com/jg_liao

...llow the normal distribution. However I know that the variable pumps is distributed along Poisson distribution. How can I model that? I would like to use LOG or SQRT transformation, but I don't know how. ***************************************** Thank you very much for all your help. Martina Pavlicova -- Department of Biostatistics Columbia University 722 W. 168th Street, 6th floor New York, NY 10032 Phone: (212) 305-9405 Fax: (212) 305-9408 Email: mp2370 at columbia.edu

...o1)~c(1:20), xlim=as.character(seq(.01, .2, by=.01)), horizontal=F, bty=n, ylim=c(-3,11), panel=function(x,y,axes=F,...){ panel.bwplot(x,y,bty=n,...) panel.bwplot(x, foo2, bty=n, col="red", ...) } ) Thank you very much for your help. Martina Pavlicova -------------------------------------------------------------------------- Department of Statistics Office Phone: (614) 292-1567 1958 Neil Avenue, 304E Cockins Hall FAX: (614) 292-2096 The Ohio State University E-mail: pavlicov at stat.ohio-state.edu Columbus, OH 43210-1247...

...: duplicate case value Rstreams.c:278: this is the first entry for that value make: *** [Rstreams.o] Error 1 ERROR: Compilation failed for package `Rstreams' > [558] -------------------------- Does anybody know what is wrong or what should I do? Thank you very much for your help. Martina Pavlicova -------------------------------------------------------------------------- Department of Statistics Office Phone: (614) 292-1567 1958 Neil Avenue, 304E Cockins Hall FAX: (614) 292-2096 The Ohio State University E-mail: pavlicov at stat.ohio-state.edu Columbus, OH 43210-1247...

Hello everybody, this is a bit off topic, so maybe you can just reply to me personally... What is the typical notation for 'skewness', 'kurtosis' and maybe 'excess kurtosis'? Thank you, Martina -------------------------------------------------------------------------- Department of Statistics Office Phone: (614) 292-1567 1958 Neil Avenue, 304E Cockins Hall

The sapply function is refusing to return a result for what seem to me to be mysterious reasons. Here is a toy example: set.seed(111) X <- list(x=runif(20),y=runif(20)) rvec <- seq(0.01,0.15,length=42) foo <- function(x,X,cc) { mean((X$x)^x + (X$y)^cc) } bar <- function(x,a,b){a+b*x} try.b <- sapply(rvec,bar,a=1,b=2) # This runs without a problem and