search for: takezawa

...mand 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 something else before this procedure? Is perl required for this? ******** E-mail: takezawa at affrc.go.jp ******** ***** http://cse.naro.affrc.go.jp/takezawa/patent-e.html ***** -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[u...

...oks 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 equation. In my humble opinion, the '^4' in the first term of the second line of this equation should be '^3'. K. Takezawa -- ***** r.otasuke@gmail.com ***** http://cse.naro.affrc.go.jp/takezawa/intro.html [[alternative HTML version deleted]]

...ble(p), work = double(2 * p), PACKAGE = "base") It may indicate that QR decomposition is used to derive an additive model instead of backfitting. I am wondering if my guess is correct, or this "the _backfitting algorithm" has another meaning. ******** E-mail: takezawa@affrc.go.jp ******** ***** http://cse.naro.affrc.go.jp/takezawa/intro.html ***** [[alternative HTML version deleted]]

...00 -1.0475699 0.2672612 I think that "my_df1" is "dffits"( http://en.wikipedia.org/wiki/DFFITS ), but in R language, "my_df2" gives the difinition of "dffits". Please let me know why. -- ***** r.otasuke@gmail.com ***** http://cse.naro.affrc.go.jp/takezawa/intro.html [[alternative HTML version deleted]]

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 derivative, s...

...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/WileyTitle...

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-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[u...

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

...E = "base") It may indicate that QR decomposition is used to derive an additive model instead of backfitting. I am wondering if my guess is correct, or this "the _backfitting algorithm" has another meaning. -- ***** r.otasuke@gmail.com ***** http://cse.naro.affrc.go.jp/takezawa/intro.html [[alternative HTML version deleted]]

...1 2 3 1.00000000 0.08891993 0.03258845 "re2" 0.72777002 0.16019126 0.02373425 11.11111 14.32524 13.12762 Apparently, "re1" is not close to "re2". I would like to know why. -- ***** r.otasuke@gmail.com ***** http://cse.naro.affrc.go.jp/takezawa/intro.html [[alternative HTML version deleted]]

Hello, I'm using mvpart option xv="1se" to compute a regression tree of good size with the 1-SE rule. To better understand 1-SE rule, I took a look on its coding in mvpart, which is : Let z be a rpart object , xerror <- z$cptable[, 4] xstd <- z$cptable[, 5] splt <- min(seq(along = xerror)[xerror <= min(xerror) + xvse * xstd]) I interprete this as following: the

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),