[R] Help with yaImpute

From: Wesley Roberts <wroberts_at_csir.co.za>
Date: Mon, 28 Jul 2008 11:43:54 +0200


Dear fellow R-users

I am trying to do some imputation using K-NN with the yaImpute library. All seems to be going well until I try to use AsciiGridImpute. All my data are correctly formatted and I am able to run and view the results of yai. Below is my code:



library("yaImpute")
data <- read.csv("Data.csv", header=TRUE, sep=",") y <- subset(data, select = c(BA_PER_HA_)) x <- subset(data, select = c(IK1_mean,IK2_mean,IK3_mean,IK4_mean))

RF <- yai(x = x, y = y, method = "randomForest")

xfiles <- list(IK1_mean = "IK1_mean.asc", IK2_mean = "IK2_mean.asc", IK3_mean = "IK3_mean.asc", IK4_mean = "IK4_mean.asc")

outfiles <- list(BA_PER_HA_ = "BasalAreaPerHA.asc")

AsciiGridImpute(RF, xfiles, outfiles)


And here is the error I receive

################################################################

> source("K-NN-Basal-Area.R")

Loading required package: randomForest
randomForest 4.5-25
Type rfNews() to see new features/changes/bug fixes. Rows per dot: 10 Rows to do: 1062
ToDo: ..........................................................................................................
Done:
First six lines of predicted data for map row: 1
       object.yRefs...toKeep.
1x0001                  -9999
1x0002                  -9999
1x0003                  -9999
1x0004                  -9999
1x0005                  -9999
1x0006                  -9999

Error in AsciiGridImpute(RF, xfiles, outfiles) :   BA_PER_HA_ is not present in the predicted data
################################################################

The first six lines of the output file contain the header information for an ESRI ASCII grid and look like this ncols 1564
nrows 1062
xllcorner -77876.53335875
yllcorner -3304144.3940076
cellsize 4.000718
NODATA_value -9999

Which is exactly the same as the input ascii grids used during the AsciiGridImpute command.

I was wondering if anyone on the R mailing list had ever come across this error and if so how did you correct it.

I am using R version 2.7.1 on Ubuntu 7.10

Many thanks for all your help.
Kind regards,
Wesley Roberts

-- 
This message is subject to the CSIR's copyright terms and conditions, e-mail legal notice, and implemented Open Document Format (ODF) standard. 
The full disclaimer details can be found at http://www.csir.co.za/disclaimer.html.


and is believed to be clean.  MailScanner thanks Transtec Computers for their support.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Mon 28 Jul 2008 - 15:40:33 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 28 Jul 2008 - 16:32:36 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive