similar to: bandwidth estimation using bw.SJ

Hello All, I used sample(..., replace = TRUE) repeatedly to generate 100 bootstrap samples of a data that I'm working with. If I re-run my codes, it'll generate a different set of 100 bootstrap samples. Is there any way that I can get the same set of 100 bootstrap re-samples at each re-run of my codes? Something like setting the seed in a random number generator... Any help is

Dear All, Is there a way to find bandwidths for Normal, and Triangular Kernels using Sheather Jones plug-in method. If not, is there a way we could convert the bandwidth obtained by the default( I assume that the default is standard normal) in hcj(x,...) to other kernels. I appreciate your help. Thanks Durai --------------------------------- [[alternative HTML version

Dear all, I have a data matrix with 7 independent variables, and each of them is a trinary variable ( - 1, 0 , 1), and I would like to know what kinds of R package or method I should use to perform the clustering. Thanks for any comment or suggestion. -- Lam C. Tsoi (Alex) Medical University of South Carolina [[alternative HTML version deleted]]

Hello, I am doing a simple regression using lm(Y~X). As my response and my predictor seemed to be skewed and I can't meet the model assumptions. Therefore I need to transform my variables. I wanted to ask what is the preferred way to find out if predictor and/or response needs to be transformed and if yes how (log-transform?). I found a procedure in "A modern approach to Regressoin in

Does anyone know how to do a fry plot in R? I have 600-800 points per image, and I really don't want to attempt that manually. Thank you! Cassie [[alternative HTML version deleted]]

Hi Folks, I've got a question concerning the calculation of the Schwarz-Criterion (BIC) done by summary.regsubsets() of the leaps-package: Using regsubsets() to perform subset-selection I receive an regsubsets object that can be summarized by summary.regsubsets(). After this operation the resulting summary contains a vector of BIC-values representing models of size i=1,...,K. My problem

Hello, I am using the lmomco package (lmom.ub and pargev) to compute the GEV parameters (location, scale, and shape), which are used to estimate return values. I was wondering how/if I can calculate upper and lower confidence (CI_u, CI_l) intervals for each return frequency using the GEV parameters to fill-in the table below? Xi (location) = 35.396 Alpha (scale) = 1.726 Kappa (shape) =

Hola a todos, Estoy utilizando la libreria e1071 para clasificar unos datos. Me gustaría poder guardar los modelos en el disco duro y no en memoria de R. He visto que hay una función: save : que guarda el modelo en memoria y load: que carga ese modelo Ejemplo: #saving the best model save(calibrate.rf.model1, file=''bestmodel.rda'') #loading the best model

Hello everyone I hope Im not bothering you all again. I have just begun to use R and so Im not yet familiarized with it.. I ve got an assignment which consists in calculating the area of a circle given a certain radius and center using the monte carlo method, which means that I have to plot a circle given its parameters. Limit the area inside it...with as many sample points as possible...and

Everytime I see this in the logs my internet access stops for about 5 minutes. Nothing else follows in the logs. I have DSL. Any ideas? Aug 15 17:14:56 bonsai gconfd (root-6095): starting (version 2.14.0), pid 6095 user 'root' Aug 15 17:14:56 bonsai gconfd (root-6095): Resolved address "xml:readonly:/etc/gconf/gconf.xml.mandatory" to a read-only configuration source at position

Context: I'm relatively new to R and am working with very large datasets. General problem: If working on a dataset requires that I produce more than two objects of roughly the size of the dataset, R quickly uses up its available memory and slows to a virtual halt. My tentative solution: To save and remove objects as they're created, and load them when I need them. To do this I'm

R-help list, I ported some R code from MS Windows PC to sun-solaris. Both do queries on the same MySQL database. PC version of R is 1.6.1 and solaris version is 1.6.0. I use RMySQL_0.5-0 (on solaris) and RODBC 1.0-1 on PC to connect to MySQL. The PC version took about 65 mins. to run wheras the Solaris (SunOS 5.7) version took 375 mins. (> 6 hrs). R and RMySQL are resident on the same host

Hi everyone, Has anyone written a parser in Java for either the ASCII or binary format produced by save()? I need to parse a single large 2D array that is structured like this: list( "32609_1" = c(-9549.39231289146, -9574.07159324482, ... ), "32610_2" = c(-6369.12526971635, -6403.99620977124, ... ), "32618_2" = c(-2138.29095689061, -2057.9229403233, ... ),

Dear all, I am doing a simulation for my model that works when I use only the rq() command. However, since I need to use the varcov matrix for my Wald test, I need to compute summary(rq(), cov=TRUE). But the simulation does not work because of the error: tau + h > 1: error in summary.rq I tried to use: if (tau + h > 1) stop("tau + h > 1: error in summary.rq") But the

Hello, I'm currently trying to model the movement of a slope (v.obs) with a regression model. The data can be found following the given links: either http://www.sendspace.com/file/dnugwc or http://rapidshare.com/files/420569660/sel.day.txt I want to use the Box-Cox transformation to normalize the response as well as the predictor variables. The scatterplot looks like this: library(zoo)

Hi all, I am wondering how to write a 'list' object to a file. I already gone through some threads like http://mail.python.org/pipermail/python-list/2001-April/080639.html, however could not trace out any reliable solution. I tried following : > write.table(calc, file="c:/data.csv") Error in data.frame("200501" = c(-0.000387071806652095, -0.000387221689252648, :

Dear R users, Thanks to a new contributor, Luigi Cerulo, amap provides 2 more metrics, and another agglomeration linkage for distance computation, and hierarchical clustering. The new centroid linkage method implemented in the cluster3 software tool (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm). In this method the pseudoitem of a cluster is computed as the average of the items

$B!g!g!g!g!y!g!g!g!g!g!y!g!g!g!g!g!g!g!g!g!y!g!g!g!g!g!y!g!g!g!g!g(B $B!c(BNO.I don't veceive your mail$B!d"M!!(Becchi_daisuki69@yahoo.com.au $B!c:#8e!"<u?.$r5qH]$9$k>l9g$O!d"M(Bsweets85433@yahoo.com.au $B!g!g!g!g!y!g!g!g!g!g!y!g!g!g!g!g!g!g!g!g!y!g!g!g!g!g!y!g!g!g!g!g(B

Dear all, Sorry if I am posting to the wrong place, but I could not find the link for registration on the bug tracker, that?s why I am writing here. I think there is inconsistency between two R functions from the stats package, bw.nrd0 and bw.nrd. Consider the following vector: D <- c(0, 1, 1, 1, 1) bw.nrd(D) returns zero bandwidth for this object even without a warning. Considering the

I need a "complete" list of "reliable" services for repeating the signal of a station to the listeners (live radio). The situation is that a station will be doing live radio and has a little DSL connection and external bandwith will be required. I don't mind if it costs much or it is cheap, I now just want to join information to say my customer: You have all this coices,