Displaying 2 results from an estimated 2 matches for "ppt1".
Did you mean:
opt1
2010 Oct 12
0
Stripes in the plot
...er = FALSE) # the size is 5764(pixels) x 528(periods)
ppt <- (ppt/100)*25.4 #convert from hundredths of an inch to mm
# the range of the area
plat <- as.matrix(seq(38.16836,41.29575,by=0.07273)) #dimension is (44,1)
plon <- as.matrix(seq(-103.9752,-94.5,by=0.07273)) #dimension is (131,1)
ppt1 <- ppt[,1] #length is 5764
ppt1 <- matrix(as.numeric(ppt[,1]), nrow=131, ncol=44)
> ppt1
[1] 1.6107435 1.6964510 1.7840456 1.8703002 1.9060464 1.8331288
[7] 1.7492645 1.6657609 1.5836155 1.5029230 1.4235667 1.3453415
[13] 1.2679848 1.1911846 1.1145816 1.0377716 0...
2011 Apr 12
1
How to set the dimension of a matrix correctly?
...as.matrix(read.table('/Users/R/Code/ppt_15day.dat'))
xx <- dim(ppt) # (77,528)
plat <- seq(37.5,42,by=0.07273)plon <- seq(-105.5,-93.5,by=0.07273)pgrid <- expand.grid(x=plon,y=plat)pdim <- dim(pgrid) # (10230,2)
#plot(pgrid, cex=0.5)
lat <- coords[,1] lon <- coords[,2]
ppt1 <- ppt[,1:xx[2]] # 1:528
pptpred <- matrix(0,ncol=xx[2],nrow=1)
############# Only test one period data##################
ptemp <- ppt1[,3] ll <- which(ptemp>0)
ppt2 <- matrix(0,nrow=length(ll),ncol=3) # (lon,lat,ptemp)
ppt2[,1] <- lat[ll] # y-axis ppt2[,2] <- lon...