Displaying 1 result from an estimated 1 matches for "pvalxx".
Did you mean:
pval
2011 Apr 12
1
How to set the dimension of a matrix correctly?
...# plot(bin2) # fig2
ini1 <- max(bin1$v) ols <- variofit(bin1, fix.nugget = F,weights="cressie",ini.cov.pars=c(ini1,4)) 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 <...