Displaying 2 results from an estimated 2 matches for "termwise".
Did you mean:
termlist
2008 May 14
4
Newbie question about vector matrix multiplication
...887114
Y 0.0002765438 0.0002570286 0.0002117336
Z 0.0007887114 0.0002117336 0.0009168750
And a weight vector generated by
w<- read.table("c:/r.weights")
X Y Z
1 0.5818416 0.2158531 0.2023053
I want to compute the product of the matrix and vectors termwise to
generate a 3x3 matrix, where m[i,j]=w[i]*co[i,j]*w[j].
0.000423773 7.47216E-08 4.41255E-08
7.47216E-08 1.96566E-11 4.29229E-11
4.41255E-08 4.29229E-11 4.11045E-11
Is this possible without writing explicit loops?
Thank you,
Dan Stanger
Eaton Vance Management
200 State Street
Boston, MA...
2003 Jan 30
2
mgcv, gam
...tried to obtain this with trigonometric terms, aca:
> Rot.gam2 <- gam(cbind(Rotavirus,Total)~ s( I(sin((MesN/12)*2*pi)),
bs="cr" )+
+ s( I(cos((MesN/12)*2*pi)), bs="cr" ),
data=na.omit(Rot),
+ family=binomial)
Warning messages:
1: Termwise estimate degrees of freedom are unreliable
but this does not seem to be wholy satisfactory. Especially, the
result of plot.gam shows smooths with confidence bounds
which cannot be discriminated visually from the smooth itself!,
and the warning messages about unreliable estimated
degrees of fre...