search for: kunio

R-users E-mail: r-help@r-project.org Hi! R-users. I am just wondering what the definition of "dffits" in R language is. Let me show you an simple example. function() { library(MASS) xx <- c(1,2,3,4,5) yy <- c(1,3,4,2,4) data1 <- data.frame(x=xx, y=yy) lm.out <- lm(y~., data=data1, x=T) lev1 <- lm.influence(lm.out)$hat sig1 <-

Hi All, What do you consider to be the best book(reference) on nonparametric regression? I am currently reading the book of Kunio Takezawa(2006): "Introduction to nonparametric regression". Is the book of Hardle(1990): "Applied nonparametric regression" better? or maybe another book? This is off topic, but most of the books is using R or S-plus. Thanks Hennie -----------------------------------------...

re: Cubic splines in package "mgcv" I don't have access to Gu (2002) but clearly the function R(x,z) defined on p126 of Simon Wood's book is piecewise quartic, not piecewise cubic. Like Kunio Takezawa (below) I was puzzled by the word "cubic" on p126. As Simon Wood writes, this basis is not actually used by mgcv when specifying bs="cr". Maybe the point is that at the knot, this continuous function has continuous 1st and 2nd derivatives, but a discontinuous 3rd deri...

...Akaike information criterion) for normal error models has just been changed. Please refer to the paper below on this matter. Eq.(22) is the new definition. The essential part is RSS(n+q+1)/(n-q-3); it is close to GCV. The paper is temporarily available at the "Papers In Press" place. Kunio Takezawa(2012): A Revision of AIC for Normal Error Models Open Journal of Statistics, Vol.2. No.3 http://www.scirp.org/journal/ojs/ **Out now!** Guidebook to R Graphics Using Microsoft Windows ISBN: 978-1-118-02639-7 Click here to read more and to buy your print copy http://as.wiley.com/WileyCDA/W...

R-users E-mail: r-help at stat.math.ethz.ch Hi! I would like to produce DLL files to be linked to R objects on Windows98SE. The source files are written in Fortran77. I input the command below on R console. Rcmd SHLIB aaa.f The result is: Error: syntax error Does this mean that "Rcmd SHLIB aaa.f" contains symtax error, or "aaa.f" contains it? Or do I need to do

R-users E-mail: r-help@r-project.org I have a quenstion on "gam()" in "gam" package. The help of gam() says: 'gam' uses the _backfitting algorithm_ to combine different smoothing or fitting methods. On the other hand, lm.wfit(), which is a routine of gam.fit() contains: z <- .Fortran("dqrls", qr = x * wts, n = n, p = p, y = y *

R-users E-mail: r-help@r-project.org My understanding is that package "mgcv" is based on "Generalized Additive Models: An Introduction with R (by Simon N. Wood)". On the page 126 of this book, eq(3.4) looks a quartic equation with respect to "x", not a cubic equation. I am wondering if all routines which uses cubic splines in mgcv are based on this quartic

Hello R-list. I am a "long time listener - first time caller" who has been using R in research and graduate teaching for over 5 years. I hope that my question is simple but not too foolish. I've looked through the FAQ and searched the R site mail list with some close hits but no direct answers, so... I would like to estimate QAIC (and QAICc) for a glm fit using the

R users E-mail: r-help at stat.math.ethz.ch I really appreciate information from Dr. Ligges and Dr. Wang. I managed to create DLL files by MinGW and use them as subroutines on R. Thank you very much again. ******** E-mail: takezawa at affrc.go.jp ******** ***** http://cse.naro.affrc.go.jp/takezawa/patent-e.html *****

R-users E-mail: r-help@r-project.org Hi! R-users. http://finzi.psych.upenn.edu/R/library/mvpart/html/xpred.rpart.html says: data(car.test.frame) fit <- rpart(Mileage ~ Weight, car.test.frame) xmat <- xpred.rpart(fit) xerr <- (xmat - car.test.frame$Mileage)^2 apply(xerr, 2, sum) # cross-validated error estimate # approx same result as rel. error from printcp(fit) apply(xerr, 2,

R-users E-mail: r-help@r-project.org Hi! R-users. A simple object as below was created to see how gam() of package "mgcv" and anova() work. function() { library(mgcv) set.seed(12) nd <- 100 xx1 <- runif(nd, min=1, max=10) xx1 <- sort(xx1) yy <- sin(xx1)+rnorm(nd, mean=5, sd=5) data1 <- data.frame(x1=xx1, y=yy) fit1 <- gam(y~s(x1, k=5),

R-users E-mail: r-help@r-project.org I have a quenstion on "gam()" in "gam" package. The help of gam() says: 'gam' uses the _backfitting algorithm_ to combine different smoothing or fitting methods. On the other hand, lm.wfit(), which is a routine of gam.fit() contains: z <- .Fortran("dqrls", qr = x * wts, n = n, p = p, y = y *