Displaying 1 result from an estimated 1 matches for "newpptpr".
Did you mean:
new_ptr
2011 Apr 12
1
How to set the dimension of a matrix correctly?
...kc <- krige.conv(pptd, loc=pgrid,krige=krige.control(type.krige="OK",trend.d="2nd",trend.l="2nd",cov.pars=ols[2]$cov.pars)) pvalxx <- which(kc$predict < 0) kc$predict[pvalxx] <- 0 pptpred <- kc$predict }else{ pptpred <- 0*(1:pdim[1]) }}
newpptpred <- matrix(pptpred, nrow=62, ncol=165) # need to fit the lat/lon frame
plat <- seq(37.5,42,by=0.07273)plon <- seq(-105.5,-93.5,by=0.07273)
lat2 <- which((plat> 37)&(plat< 42))lon2 <- which((plon> -106)&(plon< -93))
newpptpred <- tpptpred[lat2,pon2]
nLevel <...