random forest - How can I subsample a SpatialPointsDataFrame in R -


i working on running randomforest. i've imported point data representing used , unused sites , created raster stack raster gis layers. i've created spatialpointdataframe of used , unused points underlying raster values attached. 

require(sp) require(rgdal) require(raster)  #my raster stack xvariables <- stack(rlist)  #rlist = list of raster layers     # reading in spatial used , unused points. ldata <- readogr(dsn=paste(path, "data", sep="/"), layer=used_avail) str(ldata@data)   #attach raster values point data. v <- as.data.frame(extract(xvariables, ldata)) ldata@data = data.frame(ldata@data, v[match(rownames(ldata@data), rownames(v)),])

next plan run random forest using data. problem is, have large data set (over 40,000 data points). need sub sample data having hard time figuring out how this. i've tried using sample() function think because have spatialpointsdatafram wont work? i'm new r , appreciate ideas.

thanks!

subsetting spatial*dataframe object simple; use

spsubset <- spobject[<sample_criterion>,]

here's example load colleges in nj, grab random sample of size=20. there's example load states , grab "new jersey".

library(rgdal) set.seed(1) # random sample of nj colleges... samplesize=20 sppoints <- readogr(dsn=".",layer="nj_college_univ_nad83njsp") spsample <- sppoints[sample(1:length(sppoints),samplesize),]  # extract nj states tiger/line file states   <- readogr(dsn=".",layer="tl_2013_us_state") nj       <- states[states$name=="new jersey",] nj       <- sptransform(nj,crs=crs(proj4string(spsample)))  # render map nj.df    <- fortify(nj) library(ggplot2) ggplot() +   geom_path(data=nj.df, aes(x=long,y=lat, group=group))+   geom_point(data=as.data.frame(coordinates(sppoints)),               aes(x=coords.x1,y=coords.x2),colour="blue", size=3)+   geom_point(data=as.data.frame(coordinates(spsample)),               aes(x=coords.x1,y=coords.x2),colour="red", size=3)+   coord_fixed() + labs(x="", y="") + theme(axis.text=element_blank())

the tiger/line file of states can found here. nj college shapefile here.


-

Comments

Post a Comment

Popular posts from this blog

php - regexp cyrillic filename not matches -

i take list of files php glob(."*.pdf") . files have cyrillic names. example: ООО «Рога и копыта» then load company names database , try find them in file names: preg_match("/$firm_name/ui", $file_name, $match); i have problem: firms if $firm_name visually similar $file_name , preg_match says there no match. mb_detect_encoding both variables return utf-8. there problem in file name, or maybe there encoding or something, can not understand, ask help. upd: $file_name='ООО «Рога и копыта»'; preg_match("/... .Рога и копыта./ui", $file_name, $match); as result, ОО «Рога и копыта» , 1 О missing
Read more

c# - OpenXML hanging while writing elements -

Image
i have program pulls data database, caches file , exports data multiple formats (excel, excel 2003, csv). i'm using openxml sdk 2.0 excel work. these export processes run in parallel (using parallel.foreach ), , amount of data can pretty large - e.g. csvs 800mb. during these larger exports, i've noticed writing of xml documents hang. instance, if have 8 exporting in parallel, @ point "pause". hang around same point: //this.writer openxmlwriter created worksheetpart. this.writer.writeelement(new cell() { cellvalue = new cellvalue(value), datatype = cellvalues.string }); when happens, pause debugger (vs2013 in case) , notice threads blocking around same portion of code - bit deeper in openxml sdk - stem call openxmlwriter.writeelement . i dug through source using justdecompile didn't find answers. appears there intermediary stream in use writing isolated storage , is, reason, blocking. underlying stream each of these filestream . here screenshot ...
Read more

sql - Select Query has unexpected multiple records (MS Access) -

i may tired, can't seem figure out happening. have 2 queries. query 1 has 10 fields. query 2 has 5 fields. fields 1 through 3 same in both queries. trying write select query , simple add field 5 query2 query1 result should query1 additional field (query2.field5). i join first 3 fields both queries , choose select records in query1 , in query2 match. don't know sql, use query designer. expect same number of records have query 1 instead have 3 times or more. can tell me doing wrong. tried reversing joins still seem same records. thanks in advance. here query: the fields in common divisionname , supplierid = vendorid , commodityname . select [part revenue exposed query p1].divisionname, [part revenue exposed query p1].commodityname, [part revenue exposed query p1].suppliername, [part revenue exposed query p1].partnumber, [part revenue exposed query p1].partdescription, [part revenue exposed query p1].bucode, ...
Read more