r - idw() or krige() Error: dimensions do not match when missing values -


functions idw() , krige() gstat package keep reporting errors when either response or predictor variable contains missing values (na), when na.action set na.omit:

require(gstat) data(meuse) coordinates(meuse) = ~x+y data(meuse.grid) gridded(meuse.grid) = ~x+y  meuse2 <- as.data.frame(meuse) meuse2[1, 'zinc'] <- na meuse2 <- spatialpointsdataframe(spatialpoints(meuse), meuse2)  # idw response var int <- idw(zinc ~ 1, meuse2, meuse.grid, na.action = na.omit) # error: dimensions not match: locations 310 , data 154  # krige response var m <- vgm(.59, "sph", 874, .04) int <- krige(zinc ~ 1, meuse2, meuse.grid, model = m, na.action = na.omit) # error: dimensions not match: locations 310 , data 154  # krige predictor var meuse3 <- as.data.frame(meuse) meuse3[1, 'dist'] <- na meuse3 <- spatialpointsdataframe(spatialpoints(meuse), meuse3) int <- krige(zinc ~ dist, meuse3, meuse.grid, model = m, na.action = na.omit) # error: dimensions not match: locations 310 , data 154

is bug? have filter our data manually , merge results original data frames? isn't there easier solution? why there na.action option then?

the na.action argument deals missing values within newdata (not locations or data).

this stated in ?idw / ?krige / ?predict.gstat

function determining should done missing values in 'newdata'. default predict 'na'. missing values in coordinates , predictors both dealt with.

there no method deal na values within locations or data (and hence error saying there 2 more values in locations data (ie. x , y coordinate of missing data point)

you can work removing location missing value

int <- idw(zinc ~ 1, meuse2[!is.na(meuse2$zinc),],newdata= meuse.grid)

Comments

Post a Comment

Popular posts from this blog

basic authentication with http post params android -

i doing basic authentication passing username , password , using basicnamevaluepair sending post params response service. my method: public stringbuilder callservicehttppost(string username, string password, string type) { // create new httpclient , post header httpclient httpclient = new defaulthttpclient(); httppost httppost = new httppost(webservice + type); httpresponse response = null; stringbuilder total = new stringbuilder(); try { url url = new url(webservice + type); /*string base64encodedcredentials = base64.encodetostring((username + ":" + password).getbytes(), base64.url_safe | base64.no_wrap);*/ string base64encodedcredentials = "basic " + base64.encodetostring( (username + ":" + password).getbytes(), base64.no_wrap); httppost.setheader("auth...
Read more

vb.net - Virtual Keyboard commands -

i wondering difference between virtual these keyboard commands is: keyeventf_extendedkey , keyeventf_keyup is. everywhere have looked gives me description based off integers , not want know each of them does. you've tagged question vb.net, these have nothing @ vb.net. they're constants defined in windows header files, use win32 api functions. as far difference, can't tell looking @ values. individual values not particularly important, that's why named identifiers used. what's important used , documentation functions tells mean. the first one, keyeventf_extendedkey , used keybdinput structure (which used along e.g. sendinput function) pass information synthesized keyboard input. if flag used, means scan code should interpreted extended key. technically, means scan code preceded prefix byte value 224 (&he0 in hexadecimal notation). the second one, keyeventf_keyup , 1 of flags available use structure. means key being released (going up)...
Read more

c++ - End of file on pipe magic during open -

i have c++ application in starting process(wireshark) following. if (fp == null){ fp = popen(processpath, "r"); //processpath process want start if (!fp){ throw std::invalid_argument("cannot start process"); } fprintf(fp, d_msg);//d_msg input want provide process } else if(fp != null){ fprintf(fp, d_msg); } the problem when execute c++ application, start wireshark error end of file on pipe magic during open what should avoid that? also tried using mkfifo create named pipe , execute it. used this: if (fp == null){ system("mkfifo /tmp/mine.pcap"); fp = popen("wireshark -k -i /tmp/mine.pcap", "r"); if (!fp){ dout << "cannot start wireshark"<<std::endl; throw std::invalid_argument("cannot start wireshark"); } input = fopen("/tmp/mine.pcap", "wb"); fprintf(input , d_msg); fclose(input)...
Read more