Displaying 1 result from an estimated 1 matches for "thininng".
Did you mean:
thinning
2005 May 19
0
Random/systematic selection of rows in a matrix
...process
param<-c(6.63,2.24,1.82,4.36)
names(param)<-c("beta","sigmasq","tausq","varphi")
#Simulations loop - Function GaussRF from package RandomFields
for(i in 1:N.sim){
mcola<-GaussRF(x=x,y=y,param=param,grid=TRUE,model="gauss")
#Process thininng out with 'pcell': probability that the process will manifest itself
for(j in 1:540){
for(k in 1:180){
if(runif(1)<pcell) mcola[k,j]<-mcola[k,j] else mcola[k,j]<-0
}
}
file1.out<-paste("mcolasim7",i,"txt",sep=".")
write(t(mcola),file1.out,ncol=ncol(t...