Displaying 1 result from an estimated 1 matches for "dfdbtran".
Did you mean:
dfdbtrans
2008 Mar 05
2
matrix inversion using solve() and matrices containing large/small values
...965
> coeff=1.106836*10^(-8)
> echtecoeff=c(481.46,19919.23,-93.41188,0.5939589,-0.002846272,8.030726e-6
,-1.155094e-8,6.357603e-12)/10000000
> dosis=c(0,29,70,128,201,290,396)
> dfdb <-
array(c(1,1,1,1,1,1,1,dosis,dosis^2,dosis^3,dosis^4,dosis^5,dosis^6,dosis^7),dim=c(7,8))
> dfdbtrans = aperm(dfdb)
> sigerr=sqrt(coeff*dosis^macht)
> sigmadosis = c(1:7)
> for(i in 1:7){
sigmadosis[i]=ifelse(sigerr[i]<2.257786084*10^(-4),2.257786084*10^(-4),sigerr[i])
}
> omega = diag(sigmadosis)
> infomatrix = dfdbtrans%*%omega%*%dfdb
**********************
I need the inverse...