similar to: How to specify multiple regular expressions for pattern argument

Dear all, In a plot command like xyplot(Y~X|ID*PERIOD,data=...) xyplot will generate all the possible ID*PERIOD combinations. But not all of them have data in there. So I have a lot of empty plots. How can I suppress those empty plots and ask xyplot only to generate plots actually with data. Thanks. Jun [[alternative HTML version deleted]]

Dear all, This question was asked a few years ago. Back then, the answer was NO. Just wonder if the package has been updated to make it possible. An example (Theoph is a dataset coming with R) ggplot(data=Theoph, aes(x=Time, y=conc)) + geom_point() + facet_wrap(~Subject) gives me all 12 subjects on one page. Can I break it down to several pages? say 2x2 on each page? Thanks. Jun

I'm working with some data, and am trying to generate it in the following format. state city zipcode I like pizza 0 0 0 I live in Denver 0 1 0 All the fun stuff is in Alaska 1 0 0 he lives in 66062

Hi all, > diff_operator <- "\\(" > strsplit(cond, diff_operator) [[1]] [1] "andsin" "log_angle_1_4)" > diff_operator <- "\\sin(" > strsplit(cond, diff_operator) Error in strsplit(cond, diff_operator) : invalid regular expression '\sin(', reason 'Missing ')'' When I am going to split with "("

Hello all: Is there any way to add points to a point pattern, while keeping a given minimum distance amongst this new points and the pre-existing points? Thanks, Frederico [[alternative HTML version deleted]]

Hi, I was looking for a function like "which()" but only returns the first argument. Compare: x <- c(1,2,3,4,5,6) y <- 4 which(x>y) returns: 5,6 which(x>y)[1] returns: 5 which(x>y)[1] is exactly what i need. I did use this but the dataset is too big (~18 mio. Points). That's why i need a more effective way to get the first element of a vector which is

Hi, everyone, I was trying to sort a data frame by two columns, one increasing, the other decreasing and got an error. "Error in FUN(left) : invalid argument to unary operator", The command is "BA[order(BA[1],-BA[2]),]". BA is the data frame. It was working if I used increasing on both columns. Why the decreasing symbol "-" is not working here? Thanks. -- Jun

Dear R-help members, I'm processing a large amount of MODIS data where quality assessment information is stored as an integer value for each pixel. I have to converted this number to an 8 digit binary flag to get access to the stored quality code (e.g. in2bin(165,8) = 1 0 1 0 0 1 0 1). Unfortunately, I did not manage to find a package providing a fast function to do so. I need to run

Are they the same with .RData being the newer format?? Thanks, ...Tao

Dear all, I have an R script that uses Rcpp, and I have been trying to parallelize it using mclapply (I tried with the multicore and the parallel library) Sometimes (not always, interestingly), the CPU use for each core drops, usually so that the total over all cores reaches 100%, i.e., as fast as if using just one single core fully. I tried my code directly from within emacs, and also using a

Hi, I'm trying to user lmer function from lmerTest package because, if I understood correectly, it allows to make better inference than lmer method from lme4 package. However, whatever I do I keep getting this error: Error in lme4::lFormula(formula = mark ~ ssCount + sTime+ : rank of X = 1660 < ncol(X) = 1895 any ideas what could be a problem? thanks, Srecko [[alternative HTML

Dear list, We have three time course profiles with very different scales, and we want to show them in one plot. Is it possible to have three y axis? I guess not, then what would be other options? something like two 2-y axis plots on a three dimensional view? Appreciate any comment. Jun Shen [[alternative HTML version deleted]]

Hello, there!   Basically my problem is very clear. I would like to take a (numerical) integration of a function f(x,y) which can be quite complex of x and y, over a disk (x-a)^2+(y-b)^2<= r^2 (with r constant). However, after some search in R, I just cannot find a function in R that suits my purpose. Function Integrate applies to one dimensional, and adaptIntegrate to rectangle. In my case, it

R-help, I have the following file I want to import to R (some lines removed) Calibrated CTD data for station:00280001 Calibrated:23/8 2001, Salinity Unsmoothed, Fluorescence Uncalibrated Maximum observed depth: 36 m QUAL has one digit for each of pressure, temp., sal. and fluor. QUAL=1:Uncal., QUAL=2:OK, QUAL=6:Interp., QUAL=9:No data

Hello r-help, I've been banging my head against the computer in an attempt to learn how to divide my matrix into segments by rows. I want to be able to return each segment as a newly named object. I've tried looking at the apply functions and creating a for loop but brain no work. Here's the basic starting objects that I believe would be needed to separate the matrix. mat <-

Hi, I have run a groups of t test and obtained a list of the t-test results (about 30). How do I extract p-values from the 30 t-test results and make a table output? Thanks Jun Example t-test output: ============================ $Lambda_z Welch Two Sample t-test data: x[PA1.BS.DN$Onecode == 2] and x[PA1.BS.DN$Onecode == 3] t = -0.1128, df = 33.635, p-value = 0.9109 alternative

Dear R users, I am running simulations (1000), and in my simulation I am looking at specific sums. For example, if the sum is >=4 then count this, if say <3, then don't count, if the sum=3, then generate a random number from uniform distribution, if this number is say less than 0.5, then count this sum, if greater than 0.5, then don't count. I am having trouble with introducing this

Hi all: My OS is 32bit winxp,but I wanna install 64bit R2.14.1. >From the following website,it says "You can also go back and add 64-bit components to a 32-bit install, or vice versa" http://cran.r-project.org/bin/windows/rw-FAQ.html#Can-both-32_002d-and-64_002dbit-R-be-installed-on-the-same-machine_003f Does it mean that I can install and run 64bit R2.14.1 under 32bit

Hi, everyone, Try the following command to see if you get TRUE or FALSE. I get FALSE on a unix platform but TRUE on Windows. Any comment? all(paste(Theoph[1],Theoph[2])==paste(Theoph[[1]],Theoph[[2]])) -- Jun Shen PhD PK/PD Scientist BioPharma Services Millipore Corporation 15 Research Park Dr. St Charles, MO 63304 Direct: 636-720-1589 [[alternative HTML version deleted]]

Hi, all, I am thinking to compare a group of slopes from regression lines to see if they are different overall, and then make specific comparisons between groups. How can I achieve that in R? I searched the archives and there are only discussions about comparing two lines a time. Thanks. A sample data set is like the following. I would like to compare the regression slopes between the five