Displaying 8 results from an estimated 8 matches for "otasuke".
2008 Oct 19
2
definition of "dffits"
...3 4 5
-1.3333333 0.4082483 0.6000000 -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]]
2007 Dec 18
2
"gam()" in "gam" package
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 *
2007 Dec 18
1
R-users
...qraux = double(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.
--
***** r.otasuke@gmail.com *****
http://cse.naro.affrc.go.jp/takezawa/intro.html
[[alternative HTML version deleted]]
2007 Dec 26
1
Cubic splines in package "mgcv"
...ith 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]]
2008 Oct 01
0
xpred.rpart() in library(mvpart)
...print(re2)
}
The result is:
"re1"
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]]
2008 Sep 16
1
1-SE rule in mvpart
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
2008 Jul 26
0
gam() of package "mgcv" and anova()
...of smooth terms:
edf Ref.df F p-value
s(x1) 3.367 3.867 1.992 0.104
I do not understand why the value of "Ref.df" is larger than "edf"
by 0.5. And I failed to find out an equation to derive the F value(1.992).
Any help would be greatly appreciated.
***** r.otasuke@gmail.com *****
[[alternative HTML version deleted]]
2012 Jul 06
1
Definition of AIC (Akaike information criterion) for normal error models
Dear R users (r-help@r-project.org),
The definition of AIC (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