similar to: assign point values as labels in spplot

I am trying to write a loop to forecast realized volatility over successive days for the purpose of VaR prediction using the HAR-RV-CJ model which is as follows: log(RV_t+1) = ?_0 + ?_CD log(CV_t) + ?_CW log(CV_t-5) + ?_CM log(CV_t-22) + ?_JD log(J_t + 1) + ?_JW log(J_t-5 + 1) + ?_JM log(J_t-22 + 1) + e_t where RV is realized volatility, CV is continuous volatility and J is the jump which is

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",

Dear All, Is there anybody has the experience in using variogram(gstat) ? Please kindly give me some hints about the results. I used variogram() to build a semivariogram plot as: tr.var=variogram(Incr~1,loc=~X+Y,data=TRI2TU,width=5) then fir the variogram to get the parameters as: v.fit = fit.variogram(tr.var,vgm(0.5,"Exp",300,1)) v.fit model psill range 1 Nug 1.484879

In the plot below, there are some grid cells that have values below 10, which is the lowest "cut" value I have specified. Is there a way, without adjusting the number of cuts, to tell R to fill in those cells with the lowest possible color (in this case greeen)? There is a white "hole" in the image about a quarter of the way in from the left side, this is what I would like to

Hi, please help!! I've been trying to fit variogram models using the REML method in the gstat package. Every time the Windows GUI crashes. For example library(gstat) data(meuse) x <- variogram(zinc ~ 1, ~x + y, meuse) v <- vgm(140000, "Sph", 800, nug = 10000) plot(x, model = fit.variogram(x, model = v, fit.method=5)) Other fit methods are non problematic (eg. fit.method=7

Dear R-helpers, I'm having trouble compiling R packages from source on Win2K. I installed R 1.7.1beta [1] on my win2k machine [2], downloaded the fields package as source [3] and tried but failed to install the package [4]. I get the same problem with other packages, so it doesn't appear to be limited to fields. Installation of precompiled packages seems to work fine. I've been

I am returning to using R after a little time away. Saw some great new functions in the fields package and decided to upgrade. Everything seems to go fine until I tried to actually load the package, when I get an error: >library(fields) Error in parse(file, n, text, prompt) : syntax error on line 227 And (of course) I am unable to access any of the functions in the fields package. Below is

Dear list, I get the following error message when installing the package GRASS. missing link(s): interp.new missing link(s): vgram Krig predict.Krig predict.se.Krig See also below for the full installation message. What does this mean exactly? Regards, Ulrich ------------------------------------------------------------------ > R CMD INSTALL

Hi, I have a 4 column file (long/lat/elev/variable) and I tried to fit the values of my variable to the XYZ space using Tps and I keep getting the following message: Warning messages: 1: GCV search gives a minumum at the endpoints of the grid search in: Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, 2: GCV search gives a minumum at the endpoints of the grid search in:

FreeBSD 5.2.1 Upgraded to dovecot 0.99.12.1 After authentication failed, I also upgraded OpenLDAP to 2.2.19 Dovecot and OpenLDAP are now the most recent versions in the ports. I can access the LDAP info through the OpenLDAP client without a problem, so slapd seems to be working fine. Both dovecot and OpenLDAP worked fine before the upgrade. Any help will be appreciated, since several people

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints

hello I try to build DEM using Krig function of fields package. And I get this error message. here is the procedure I followed: > dt<-read.table("/usr/local/bartin/stat/topostat1",header=F,sep="|") > names(dt) <-c("x","y","z") > coord<-cbind(dt$x,dt$y) > elevation<-cbind(dt$z) >

I'm new at this, so please bear with me. I am trying to get Krig (actually tps) results into a data frame at a set of points. Frankly, all I really need is to query the Krig on a specific XY pair and return a result. I can carry the rest. I am using tps from the fields package. Here's the code I use to create the Krig object: =========== Snip! =============================== dCurr

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Mike J Smith > Sent: 09 June 2005 09:58 > To: r-help at stat.math.ethz.ch > Subject: [R] krig.image help > > > Hi > > I have recently been experimenting with the use of kriging, primarily > through Goldensoftware's Surfer.

Hi all, >str(df) `data.frame': 31837 obs. of 3 variables: $ x : num 410683 410700 410720 410740 410324 ... $ y : num 43136 43126 43123 43125 42709 ... $ wz: num -101.1 -94.9 -93.3 -94.5 30.8 ... >library(gstat) >g<-gstat(id="rv",form=wz~1,loc=~x+y,data=df,model=mat,nmax=500,set=list(average=1)) >str(g) List of 3 $ data :List of 1 ..$ rv:List of 10 ..

Servus Philipp, Unfortunately the traceroute(6) results are both more or less random. Sometimes traceroute "hangs" a while, wherever, sometimes traceroute6. Sometimes traceroute is faster, sometimes traceroute6. Not reliable. Your MTU question, tried as adviced: Maximum size for IPv4 is 1466, and max size for IPv6 is 1444. Exceeding these values leads to a "ping: local error:

Hi R community I need to interpolate precipitation data for a natural park. I have precipitation data from some climate stationts. (I know the table is not complete but I only need to show you X,Y, Altitude and PrepJul) X Y Altitude PrepJan PrepFeb PrepMar PrepAp PrepMay PrepJun PrepJul 597706 4093438 41 0 9 77,8 63,1 17,5 0 2,6 597535

I recently installed 4.8-stable. When running lynx, I can only access sites outside my LAN by launching lynx as root. I can't access external sites if I launch lynx using sudo or after su-ing to root. I've searched archives and googled, and haven't found anything on this. It happened immediately after the upgrade, before any post-installation configuration, i.e. before setting up

Thanks for all the helpful responses. I include the data file and the syntax file for reference. Again, if I use the fields function, as is, I get the message: Error in svd(tempM) : error 159 in dsvdc using traceback, I get: > traceback() 4: stop(paste("error ", z$info, " in dsvdc")) 3: svd(tempM) 2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp,