आर पैकेज "spatstat": po.gon = TRUE जब एक poisson प्रक्रिया मॉडल (फ़ंक्शन: पीपीएम) में गैर-चिकनी शर्तों के लिए आपको मानक त्रुटियां कैसे मिलती हैं?

Question

आर पैकेज में spatstat (मैं वर्तमान संस्करण, 1.31-0 का उपयोग कर रहा हूं), एक विकल्प use.gam है। जब आप इसे सत्य पर सेट करते हैं, तो आप रैखिक भविष्यवाणी में चिकनी शर्तों को शामिल कर सकते हैं, वैसे ही आप आर पैकेज mgcv के साथ करते हैं। उदाहरण के लिए,

g <- ppm(nztrees, ~1+s(x,y), use.gam=TRUE) 

अब, अगर मैं अवरोधन के लिए एक विश्वास अंतराल चाहते हैं, तो आप आमतौर पर summary या vcov उपयोग कर सकते हैं, जो काम करता है आप gam का उपयोग नहीं करते जब लेकिन विफल रहता है आप का उपयोग गम

vcov(g) 

करते हैं जब

जो त्रुटि संदेश

Error in model.frame.default(formula = fmla, data = 
    list(.mpl.W = c(7.09716796875, :invalid type (list) for variable 's(x, y)' 

मैं वाकिफ हूँ देता है कि इस मानक त्रुटि सन्निकटन यहाँ उचित क नहीं है एन आप gam उपयोग करें, लेकिन यह चेतावनी संदेश द्वारा कब्जा कर लिया गया है: - मैं उद्देश्य मैं उन्हें उपयोग कर रहा हूँ के लिए इन मानक त्रुटियों के उपयोग का औचित्य साबित करने के लिए तैयार कर रहा हूँ

In addition: Warning message: model was fitted by gam(); 
      asymptotic variance calculation ignores this 

मैं इस बारे में चिंतित नहीं हूँ - मुझे सिर्फ संख्याएं चाहिए और ऐसा करने के लिए "लेखन-मेरा-खुद" से बचना चाहूंगा।

मुझे जो त्रुटि संदेश मिला है वह उस डेटा सेट पर निर्भर नहीं है जो मैं उपयोग कर रहा हूं। मैंने nztrees उदाहरण का उपयोग किया क्योंकि मुझे पता है कि यह spatstat के साथ पूर्व लोड हो गया है। ऐसा लगता है कि यह चर के बारे में शिकायत कर रहा है, लेकिन मॉडल स्पष्ट रूप से सिंटैक्स को समझता है क्योंकि यह मॉडल (और अनुमानित मूल्य, मेरे अपने डेटासेट के लिए, काफी अच्छा दिखता है, इसलिए मुझे पता है कि यह सिर्फ कचरा पंप नहीं कर रहा है)।

क्या किसी के पास इस बारे में कोई सुझाव या अंतर्दृष्टि है? क्या यह एक बग है? मेरे आश्चर्य के लिए, मैं इस ऑनलाइन के बारे में कोई चर्चा नहीं ढूंढ पाया। किसी भी मदद या संकेत की सराहना की जाती है।

संपादित करें: हालांकि मैंने निश्चित रूप से यहां अपने स्वयं के प्रश्न का उत्तर दिया है, मैं समय के लिए अपना उत्तर स्वीकार नहीं करूंगा। इस तरह, अगर कोई दिलचस्पी लेता है और spatstat के अगले संस्करण की प्रतीक्षा किए बिना इसके लिए "वर्कअराउंड" खोजने के प्रयास में इच्छुक है, तो मैं उसे बक्षीस दे सकता हूं। अन्यथा, मैं केवल बक्षीस अवधि के अंत में अपना उत्तर स्वीकार करूंगा।

Answer 1

मैंने इस बारे में पैकेज लेखकों में से एक एड्रियन Baddeley से संपर्क किया है।उन्होंने तुरंत जवाब दिया और मुझे बताया कि यह वास्तव में सॉफ्टवेयर के साथ एक बग है और जाहिर है, मैं इसे सामना करने वाला पहला व्यक्ति हूं। सौभाग्य से, इस मुद्दे को ट्रैक करने और इसे सही करने के लिए उसे केवल थोडा समय लगा। फिक्सस्टैट, 1.31-1 की अगली रिलीज में फिक्स शामिल किया जाएगा।

संपादित करें:spatstat के अद्यतन संस्करण जारी किया गया है और इस बग अब और नहीं है:

g <- ppm(nztrees, ~1+s(x,y), use.gam=TRUE) 
sqrt(vcov(g)[1,1]) 
[1] 0.1150982 
Warning message: 
model was fitted by gam(); asymptotic variance calculation ignores this 

अन्य रिलीज नोट्स के लिए the spatstat website देखें। इस धागे में पढ़ने और भाग लेने वाले सभी लोगों के लिए धन्यवाद!

Answer 2

मुझे यकीन नहीं है कि आप इस तरह के रुझान को निर्दिष्ट कर सकते हैं जो संभवतः त्रुटि उत्पन्न कर रहा है। यह सुनिश्चित करने के लिए प्रतीत नहीं होता भावना according to the documentation:

The default formula, ~1, indicates the model is stationary and no trend is to be fitted.

लेकिन बजाय तुम इतनी तरह मॉडल निर्दिष्ट कर सकते हैं:

g <- ppm(nztrees, ~x+y, use.gam=TRUE) 
#Then to extract the coefficientss: 
>coef(g) 
(Intercept)    x    y 
-5.0346019490 0.0013582470 -0.0006416421 
#And calculate their se: 
vc <- vcov(g) 
se <- sqrt(diag(vc)) 
> se 
(Intercept)   x   y 
0.264854030 0.002244702 0.003609366 

क्या यह मेकअप भावना/अपेक्षित परिणाम? मुझे पता है कि पैकेज लेखक r-sig-geo मेलिंग lsit पर बहुत सक्रिय हैं क्योंकि उन्होंने अतीत में मेरी मदद की है। आप अपने प्रश्न को उस मेलिंग सूची में भी पोस्ट करना चाहते हैं, लेकिन जब आप ऐसा करते हैं तो आपको यहां अपने प्रश्न का संदर्भ देना चाहिए।