search for: sheather

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 delet...

...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 R" (Sheather, 2009): There they suggest an approach with the function bctrans from alr3...but it seems that it is deprecated. So what is the best way (box-cox test) find the best transformation for predictor and response simultaneously? AFAIK boxcox from MASS is used only used for transformation of the predicto...

Hello All, When I use bw.SJ (based on Sheather & Jones, 1991) in R to estimate the bandwidth for a highly skewed data, I get the following message: "sample is too sparse to find TD". I played around with the parameters such as no. of bins (nb), lower, upper (range over which to minimize) to no avail. My sample size is 50,000....

...r of variables i ^^ This seems to match the calculation done by extractAIC but it doesn't! Maybe anyone can tell me about the reason of the variation of the BIC-values? Best regards, Jan Henckens ### Minimal Example: require(leaps) bridge <- read.table("http://www.stat.tamu.edu/~sheather/book/docs/datasets/bridge.txt", header=TRUE) fmla.full <- formula(Time ~ .) (lm.model <- summary(regsubsets(fmla.full,data=bridge,weights=NULL, intercept=TRUE, method="forward"))) lm.model$bic ### The first two models constructed via lm(): extractAIC(lm(Time~Dwgs,data=bridge)...

...nce 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 Hall-Sheather bandwidth is very high because I also vary the number of observations from 40 to 300. Is there anyone that could help me? Thanks in advance! All the best, Julia [[alternative HTML version deleted]]

...ctor variables. The scatterplot looks like this: library(zoo) library(alr3) load("sel.day.txt") sel.p1<-window(sel, start=as.POSIXct("2008-04-05"), end=as.POSIXct("2009-04-01")) pairs(~v.obs+ snow+ HH6.1+ Q.Enz+ pcpt+ qd,data=sel.p1,gap=0.4,cex.labels=1.5) In Sheather: "A Modern Approach to Regression with R" the function bctrans is used to calculate lambda for the variables. I use "yeo.johnson" since there are values=0 in the data. Doing this creates following output: 2> summary(bctrans(~v.obs+ snow+ pcpt+ Q.Enz+ qd+ HH6.1, data=sel.p1...

You could always use 40-gigabit between the two storage systems which would speed things dramatically, or back to back 56-gigabit IB. ---------------------------------------- From: zfs-discuss-request@opensolaris.org Sent: Monday, March 18, 2013 11:01 PM To: zfs-discuss@opensolaris.org Subject: zfs-discuss Digest, Vol 89, Issue 12 Send zfs-discuss mailing list submissions to