search for: arecibo

Dear useRs, I am trying to build my package (nonpareff) which deals with some models of data envelopment analysis. The building worked well, but checking complains when it tests examples. Zipped nonparaeff.Rcheck is attached. Following is the log. --------------------------------------------- arecibo:tmp arecibo$ R CMD build nonparaeff/ * checking for file 'nonparaeff/DESCRIPTION' ... OK * preparing 'nonparaeff': * checking DESCRIPTION meta-information ... OK * checking whether 'INDEX' is up-to-date ... NO * use '--force' to overwrite the existing 'INDEX'...

Dear useRs, I am trying to build my package (nonpareff) which deals with some models of data envelopment analysis. The building worked well, but checking complains when it tests examples. Zipped nonparaeff.Rcheck is attached. Following is the log. --------------------------------------------- arecibo:tmp arecibo$ R CMD build nonparaeff/ * checking for file 'nonparaeff/DESCRIPTION' ... OK * preparing 'nonparaeff': * checking DESCRIPTION meta-information ... OK * checking whether 'INDEX' is up-to-date ... NO * use '--force' to overwrite the existing 'INDEX'...

Dear UseRs, I would like to know the way to create a new column by naming it simultaneously. For example, in for() loop I want to create columns named as paste("test", i, sep = ""), as shown below. ---------------------------- dt <- data.frame(a = c(1, 2, 3), b = c(3, 2, 2), c = c(1, 3, 5)) for(i in 1:2){ dt$as.name(paste("test", i, sep = ""))

Hi useRs! I would like to know how to make aggregated data.frame with aggregate() tabulated. For example, I run the following command to aggregate re with respect to group1 and group2. > (aggr <- with(final, aggregate(re, group1, group2, mean))) Group.1 Group.2 x 1 1992 15 0.16392 2 1993 15 0.15467 3 1994 15 0.15456 4 1995 15 0.15391

Hi expeRts, I would like to calculate weighted mean by two factors. My code is as follows: R> tmp <- by(re$meta.sales.lkm[, c("pc", "sales")], re$meta.sales.lkm[, c("size", "yr")], function(x) weighted.mean(x[,1], x[,2])) The result is as follows: R> tmp size: micro yr: 1994 [1] 1.090

Dear UserRs, I draw filled contour plot with filled.contour function. By default, colors of filled.contour is terrain.colors which alters from green (the lowest level) to white (the highest level). I want the colors to alter from white (the lowest level) to green (the highest level). How can I set this color levels? Thanks in advance.

Dear UseRs, I would like to know the way to find classes of each column of data.frame(). Thank you in advance. ========================================================= Dong-hyun Oh Center of Excellence for Science and Innovation Studies Royal Institute or Technology, Sweden e-mail: oh.dongh at gmail.com cel: +46 73 563 45 22

Dear UseRs, I wrote following function in order to solve Data Envelopment Analysis. Reason for posting is that the function is slow when nrow(dat) is large. I wonder if other functions could substitute the for() loop in the code, such as mapply(). Can anybody help to rewrite the dea() function as efficiently as possible? The code is as follows:

Dear UseRs, I would like to know which function is the most efficient in finding convex hull of points in 3(or 2)-dimensional case? Functions for finding convex hull is the following: convex.hull (tripack), chull (grDevices), in.chull (sgeostat), convhulln (geometry), convexhull.xy (spatstat), calcConvexHull (PBSmapping). I also would like to know if there is a function that can be used

Dear expeRt, I occasionally use help.search() function to find appropriate functions in R. Whenever committing help.search() after rerunning R, it takes so much time to find help (around 2 minutes). Is there any way to make it faster in help.search()? I think indexing help database is one of the solution, but i don't know how to do. Best, Dong-hyun Oh

Dear expeRt, I would like to know whether a function similar to str_replace() in php exists in R. Looking forward to hearing from you. Best, ========================================================= Center of Excellence for Science and Innovation Studies, Royal Institute of Technology Drottning Kristinas v?g 30 100 44 Stockholm, Sweden e-mail: donghoh at infra.kth.se, doriaba2 at

Hi, my data set is data.frame(id, yr, y, l, e, k). I would like to estimate Lee and Schmidts (1993, OUP) model in R. My colleague wrote SAS code as follows: ** procedures for creating dummy variables are omitted ** ** di# and dt# are dummy variables for industry and time ** data a2; merge a1 a2 a; by id yr; proc sysnlin maxit=100 outest=beta2; endogenous y; exogenous l e k

Dear UseRs, I would like to know the created time and date of specific object. Is there any function for it? Thank you in advance. Sincerely,

Dear expeRts, I would simplify following code. --------------------------------------------- youtput <- function(x1, x2){ n <- length(x1) y <- vector(mode="numeric", length=n) for(i in 1:n){ if(x1[i] >=5 & x1[i] <= 10 & x2[i] >=5 & x2[i] <=10) y[i] <- 0.631 * x1[i]^0.55 * x2[i]^0.65 if(x1[i] >=10 & x1[i] <= 15 &

Dear UseRs, I wrote the following function to use MLE. --------------------------------------------- mlog <- function(theta, nx = 1, nz = 1, dt){ beta <- matrix(theta[1:(nx+1)], ncol = 1) delta <- matrix(theta[(nx+2):(nx+nz+1)], ncol = 1) sigma2 <- theta[nx+nz+2] gamma <- theta[nx+nz+3] y <- as.matrix(dt[, 1], ncol = 1) x <- as.matrix(data.frame(1,