search for: eilers

OK, I searched through the archives and didn't find anything so I must be doing something stupid. I just compiled 2.6.3 and am trying to get it to run on linux 2.6.8.1 kernel with a more or less fedora core2 environment. I was trying between two systems but I have narrowed it down to what I believe is a pretty simple case. Here is my setup/test: Last login: Mon Oct 4 19:30:49 2004

...the choleski decomposition the fastest. Here is the function using SparseM difsm <- function(y, lambda, d){ # Smoothing with a finite difference penalty # y: signal to be smoothed # lambda: smoothing parameter # d: order of differences in penalty (generally 2) # based on code by Paul Eilers 2002 require(SparseM) m <- length(y) E <- as(m,"matrix.diag.csr") D <- diff(E,differences=d) B <- E + (lambda * t(D)%*%D) z <- solve(B,y) z } To do this in Matrix, I am not sure how to get an Identity matrix of dimension m. From looking at the documentation I...

As the request for the Savitzky-Golay Algorithm in R has come up several times, I here include my implementation based on code written for Matlab. Savitzky-Golay uses the pseudo-inverse pinv() of a matrix. There is an 'generalized inverse' ginv() in the MASS package, but I use a simpler form because I didn't want to 'require' MASS any time I apply Savitzky-Golay.

Apologies if this is this too obscure for R-help. In package splines, ns(x,,knots,intercept=TRUE) produces an n by K+2 matrix N, the values of K+2 basis functions for the natural splines with K (internal) knots, evaluated at x. It does this by first generating an n by K+4 matrix B of unconstrained splines, then postmultiplying B by H, a K+4 by K+2 representation of the nullspace of C (2 by K+4),

...urbed the lmtp process). -- J?rgen Gmach . juergen.gmach at apis.de . +49 9482 941545 APIS Informationstechnologien GmbH . http://www.apis.de Gewerbepark A 13, 93086 W?rth/Donau . Deutschland Sitz der GmbH: W?rth/Donau, Amtsgericht Regensburg (HRB 6684) Gesch?ftsf?hrer: Julia Anna Dietz, J?rgen Eilers, Peter Rosenbeck

...ariables (used for varying coefficient models) can now be factor variables, to enable easy conditioning of smooths on factors. * The default p-values for smooth terms have been substantially improved. * see ?gam.models and ?summary.gam for further details. New smoothers ---------------------- * Eilers and Marx style P-splines are now built in, along with a cyclic version. See ?p.spline. * An adaptive smoother class has been added. See ?adaptive.smooth. * The interface for adding user defined smooths has been modified and simplified. See ?smooth.construct. A fuller list of changes is at http:/...

...ariables (used for varying coefficient models) can now be factor variables, to enable easy conditioning of smooths on factors. * The default p-values for smooth terms have been substantially improved. * see ?gam.models and ?summary.gam for further details. New smoothers ---------------------- * Eilers and Marx style P-splines are now built in, along with a cyclic version. See ?p.spline. * An adaptive smoother class has been added. See ?adaptive.smooth. * The interface for adding user defined smooths has been modified and simplified. See ?smooth.construct. A fuller list of changes is at http:/...

On examining non-linearity of Cox coefficients with penalized splines - I have not been able to dig up a completely clear description of the test performed in R or S-plus. >From the Therneau and Grambsch book (2000 - page 126) I gather that the test reported for "linear" has as its null hypothesis that the spline coefficient is the same at the center of basis. Thus, in the example

...at the coefficients of the spline terms fall in a linear series. > That > is, a linear trend test on the coefficients. This is what coxph does. As > with > the LR test, the chi-square dist is conservative. I have not worked at > putting > in the more correct distribution. See Eilers and Marx, Statistical Science > 1986. > > > And what is the null for the non-linear test? > > The linear test is "is a linear better than nothing", the non-linear one is > a > sequential test "is the non-linear better than the linear". The second > t...

more specific on what your problem is? cheers, jerry --------------------------------------------------------------------- http://www.valinux.com/ VA Linux Systems gcarter@valinux.com http://www.samba.org/ SAMBA Team jerry@samba.org http://www.plainjoe.org/ jerry@plainjoe.org --"I never saved anything for the swim back."

Author: setje Repository: /hg/zfs-crypto/gate Revision: a1ec8b0f5c3c118ca5324dcc2c364e5cdb95648f Log message: 6270888 DHCP setup fails on nvidia built-in NICs on AMD64 machines and intel built-in NICs on Dimension2350s Files: update: usr/src/uts/common/io/strplumb.c

Author: setje Repository: /hg/zfs-crypto/gate Revision: a8b6ee3e951dd078ba5402d1fdc1f5fda4d4af29 Log message: 6308072 need to update boot image with one supplied by brand artist 6308073 need to suppress splashimage on cd/dvd media Files: update: usr/src/grub/install_menu update: usr/src/grub/splash.xpm.gz

Author: setje Repository: /hg/zfs-crypto/gate Revision: cc6c7f4cdc6315564dc66bc978503b27d8a4499f Log message: 6341178 bootadm refuses to check read-only alternate root Files: update: usr/src/cmd/boot/bootadm/bootadm.c

Author: setje Repository: /hg/zfs-crypto/gate Revision: bb9f8eccd7d8aef91c3bf4f6b393f812a8d0ce16 Log message: 6362643 GRUB graphics still not in line with unified "coolstart" branding 6363028 progress indicator needs to move and turn green to accommodate "coolstart" branding 6363031 boot progress graphics need to be updated with unified "coolstart" branding 6372042

Author: setje Repository: /hg/zfs-crypto/gate Revision: 538dbfddda692163b723b8101962a74a070a8b5c Log message: 6414238 no menu on x86 cd Files: update: usr/src/cmd/boot/scripts/root_archive.ksh