Displaying 2 results from an estimated 2 matches for "plon".
Did you mean:
plan
2011 Apr 12
1
How to set the dimension of a matrix correctly?
...for help,Peiling
------------------
library(geoR) #functions for geostatistical data analysis
coords<- as.matrix(read.table('/Users/R/Code/stncoords.dat'))ppt<- 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#############...
2010 Oct 12
0
Stripes in the plot
...ot;)
ppt <- read.table ('/Users/R/Code/LULC/pptnewlulc.dat', header = 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.502923...