2017-08-07 7 views
0

CRCのspatstatブック(第17章)にある小さなlinnetオブジェクトのlinearK推定値を学習しようとしていますが、linearK関数を使用するとspatstatがエラーをスローします。私は以下のrコードのコメントのプロセスを文書化しました。エラーは次のとおりです。seqのlinearKエラー。 default()をNAにすることはできません。NaN

Error in seq.default(from = 0, to = right, length.out = npos + 1L) : 'to' cannot be NA, NaN or infinite 

これを解決する方法がわかりません。私も、私はLPPオブジェクトに問題があると感じlinearK(d1_lpp)

Error in seq.default(from = 0, to = right, length.out = npos + 1L) : 'to' cannot be NA, NaN or infinite 

上のエラーを取得し、その結果、その後

# I have data of points for each data of the week 
# d1 is district 1 of the city. 
# I did the step below otherwise it was giving me tbl class 
d1_data=lapply(split(d1, d1$openDatefactor),as.data.frame) 

# I previously create a linnet and divided it into districts of the city 
d1_linnet = districts_linnet[["d1"]] 

# I create point pattern for each day 
d1_ppp = lapply(d1_data, function(x) as.ppp(x, W=Window(d1_linnet))) 
plot(d1_ppp[[1]], which.marks="type") 

# I am then converting the point pattern to a point pattern on linear network 
d1_lpp <- as.lpp(d1_ppp[[1]], L=d1_linnet, W=Window(d1_linnet)) 


d1_lpp 
Point pattern on linear network 
3 points 
15 columns of marks: ‘status’, ‘number_of_’, ‘zip’, ‘ward’, 
‘police_dis’, ‘community_’, ‘type’, ‘days’, ‘NAME’, 
‘DISTRICT’, ‘openDatefactor’, ‘OpenDate’, ‘coseDatefactor’, 
‘closeDate’ and ‘instance’ 
Linear network with 4286 vertices and 6183 lines 
Enclosing window: polygonal boundary 
enclosing rectangle: [441140.9, 448217.7] x [4640080, 4652557] units 

# the errors start from plotting this lpp object 
plot(d1_lpp) 

"show.all" is not a graphical parameter 
Show Traceback 
Error in plot.window(...) : need finite 'xlim' values 

coords(d1_lpp) 


x y seg tp 
441649.2 4649853 5426 0.5774863 
445716.9 4648692 5250 0.5435492 
444724.6 4646320 677  0.9189631 
3 rows 

そして、私はそれがハード・エラーとどのように解釈するために見つける:私はこのプロセスを以下のいそれらを解決する。誰かが私を案内してくれますか?

おかげ

+0

組み込みのデータセットまたは人工生成データを使用して、データにアクセスするか、再現可能な例が必要です。 –

+0

私はこのリンクにデータを追加しました。 [データ](https://drive.google.com/open?id=0B_2R6KanT3YIQmxFTGVUaVAxYkk)。ありがとうございました。 – BKS

答えて

0

私は線形ネットワーク上のマークされた点パターンをプロットしようとしたときにバグがplot.lppであり確認することができます。それは間もなく修正されることがうまくいくでしょう。あなたは、私がlinearKで問題を再現することはできません

plot(unmark(d1_lpp)) 

を使用してマークされていない点パターンをプロットすることができます。 spatstatのどのバージョンを実行していますか?私のラップトップの開発版ではspatstat_1.51-0.073がすべて動作します。このコードは最近変更されているため、これは開発バージョン(https://github.com/spatstat/spatstatを参照)に更新することで解決される可能性があります。

+0

私は通常のクレーンバージョンを使用しています。私はgithubからインストールして、linearKが動作するようです。私はもっ​​と実験します。ありがとうございました。 – BKS

関連する問題