Displaying 6 results from an estimated 6 matches for "lzm".
Did you mean:
lvm
2009 Jan 29
0
assign point values as labels in spplot
...es(my.dat)=~X+Y
pred.dat<-read.table("C:/temp/Predict_Location_XY.txt", header = FALSE)
names(pred.dat)<-c("x","y","p")
my.pred.dat<-data.frame(pred.dat)
coordinates(my.pred.dat)<-~x+y
gridded(my.pred.dat) = TRUE
#Set up 2 exponential variograms
lzm.vgm.exp.noNug<-vgm(1.1,"Exp",2000,0)
lzm.vgm.exp.Nug<-vgm(0.7,"Exp",3000,0.4)
#Krige the exponential variograms
lzm.krig.exp.noNug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.noNug)
lzm.krig.exp.Nug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.N...
2012 Jul 12
0
HAR-RV-CJ Moedel
...3), ]
tim = strptime(dollar[,1], "%d/%m/%Y %H:%M");
dollar$day = format(tim, "%Y-%m-%d");
##ddx <- row index for each day
ddx = split(seq.int(nrow(dollar)), dollar$day);
dollarbns = data.frame(day=as.Date(names(ddx)), m=NA, rtn=NA, rv=NA, bv=NA,
zg=NA, lzg=NA, zj=NA, zmj=NA, lzm=NA, zgtp=NA, lzgtp=NA, zmjtp=NA, zjtp=NA,
lzmtp=NA);
for (i in seq_along(ddx)) {
##quotes for day i
x = dollar$mid[ddx[[i]]];
n = length(x);
y = log(x[-1]/x[-n]);#log price difference
y2 = y^2
m = length(y);
z999 = (qnorm(1 - .001))
if (m>=260) {
##raw moments
rtn = sum...
2010 Sep 27
1
spplot cuts
Hello-
After looking through ?spplot, I would expect that I could specify the
values of the cuts:
"...?cuts? number of cuts or the actual cuts to use..."
So in the following command,
spplot(lzm.krige.dir["var1.pred"], scales=list(draw=TRUE),
xlab="Easting",ylab="Northing",
cuts=seq(0.0,0.4,by=0.01),key.space="right",cex=1.1,col.regions=terrain.colors(41),
main="Specific Yield Layer 1",sp.layout=list(pts))
I get the following warning m...
2009 Sep 10
1
sppolot: fill below minimum legend value
...quot;hole" in the
image about a quarter of the way in from the left side, this is what I would
like to correct. Thanks...Eric
The code:
pts<-list("sp.points",K.dat,pch=3,col="black")
cuts<-c(10,20,30,40,50,60,70,80,90,100,200,300,400,500,600,700,800,900,1000)
spplot(lzm.krige.dir["var1.pred"],at=cuts,colorkey=list(at=log10(cuts),at=log10(cuts),labels=as.character(cuts)),scales=list(draw=TRUE),
xlab="Easting",ylab="Northing",key.space="right",cex=1.1,col.regions=terrain.colors(30),main="Hydraulic
Conductivity of Layer 2&...
2009 Sep 11
3
For sending my R package as part of R-project
...quot;hole" in the
image about a quarter of the way in from the left side, this is what I would
like to correct. Thanks...Eric
The code:
pts<-list("sp.points",K.dat,pch=3,col="black")
cuts<-c(10,20,30,40,50,60,70,80,90,100,200,300,400,500,600,700,800,900,1000)
spplot(lzm.krige.dir["var1.pred"],at=cuts,colorkey=list(at=log10(cuts),at=log10(cuts),labels=as.character(cuts)),scales=list(draw=TRUE),
xlab="Easting",ylab="Northing",key.space="right",cex=1.1,col.regions=terrain.colors(30),main="Hydraulic
Conductivity of Layer 2&...
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...mp;DY{q}{=nS)w4O7RCQYlbf4neW
zSd&cnm*^5nQoIeKJDh4$7?cY=>#UZ>`$wWgEo_OHIfbA*O-B48Nil4PUtxTyUdo}q
z>WSzwXsgyz-M<`jBDY)ePo#e=)l;Ia{7aGFQXPFhx=BkGJj|FmT&IL9$eG547PC_C
zMa-Ea5s`Zh^|BSZ2zZME8d1PXotniQfq)2)ZqOonfvs9o8=<b7y&{II?XrCx)s^S7
z-^c|soJLdIuFe;<D$U7WXdha1l at 9!TmlZE~wLZM8L*6wad<9TEZijm)Yk)Ewb+cH~
zY#&Tv{KwTK!8!>E7G103I}uo?lMk{(V7)HW)7DCC6(q1hR}+DaVs2X^`cF)_pR1~(
z&`Z8k!)V6(F2+5iM`JJ*aI#WzQM4jhFTqWE(yEjKn{+uzL#8&jF*U;`J!S0`p=5-N
z)|(xzs|$WD1-_Oj!=++=yb>8^=;>s=R!M;>fw4^JGfO1l=7~42$<*NY2o0Xpa`GgN
zmI|TphlWeR9dgMInZ<^aC^svm at MgV1+qDYO...