Eğri uyumu sonuçlarını bilimsel bir makalede bildirme


11

(Umarım bu soru bu siteye uyar; değilse özürlerimi kabul et).

Belirli bir simülasyon çalıştırdım ve bir zaman serisi y (t), t = 0, 1, ... 20 aldım. Bazı işlevleri denedikten sonra şunu buldum:

y(t) =~ 1 / (A t + B)

Burada A ve B, R ^ 2> 0,99 ile doğrusal regresyon kullanılarak hesaplanan I katsayılarıdır.

Bu sonuçları bilimsel bir raporda bildirmenin standart yolu nedir? özellikle:

A. Teorik bir açıklamam yok, çıktının neden böyle göründüğünü biliyorum (azaldığını ve aşağıdan sınırlandığını biliyorum, ama çok daha fazla değil). Bu sadece başarılı bir tahmindi. Denediğim diğer tüm başarısız tahminleri tanımlamalı mıyım?

B. Simülasyonu her çalıştırdığımda, A ve B'nin biraz farklı değerlerini alıyorum. Sadece rastgele bir çalışma rapor etmeli miyim yoksa simülasyonu birçok kez çalıştırmalı mıyım ve sonuçları ortalama mı? Eğer öyleyse, kaç kez yeterlidir?


Ne iletmek istersiniz? Her bir simülasyon neyi temsil eder?
Bill Barth

Arazi mülkiyeti simülasyonu. N vatandaşı ve N arsa vardır. Başlangıçta, her arsa rastgele bir vatandaşa verilir. Daha sonra, her yıl, her arazi belirli bir olasılıkla p satılır ve eğer gerçekten satılırsa, alıcı rastgele seçilir. 50 yıl sonra, şu anda arazisi olmayan bazı arazilerin asıl sahiplerine iade edildiği bir "Jübile" prosedürü uygularım. Her Jübile'den (t) sonra topraksız vatandaş sayısını (y) ölçüyorum. Kesinlikle y (t) artmaz. Tahmin edilebilir bir oranda azaldığını ve 0'a yaklaştığını göstermek istiyorum.
Erel Segal-Halevi

AB

2
xnn=0Nn

Bill: Yani A ve B'yi birçok kez hesaplamam, sonra ortalama ve std'yi raporlamam gerekir mi? Bence daha iyi bir yaklaşım, tüm simülasyonlardan tüm örneklerle tek bir doğrusal regresyon yapmaktır. Ama simülasyonu kaç kez çalıştırmalıyım?
Erel Segal-Halevi

Yanıtlar:


5

Dağıtımınıza bir güç yasası uydurmaya çalışıyorsunuz . Çok ilginç. Bunlar her zaman grafik teorisinde , sosyal ağlarda ve diğer yerlerde görülür .

Verilerinizi buraya ve buraya sığdırmaya yönelik bazı eğitimler var .

Ayrıca, A sorusuna atfen, bir kişinin arazi satın alma olasılığı zaten sahip oldukları arazi miktarına nasıl bağlıdır? Bir güç yasasının verilerinize neden uygun olduğunu açıklamak için Barbasi'nin modelini kullanabilirsiniz .

güncelleme: Bunu kullandım ve harika çalışıyor: https://pypi.python.org/pypi/powerlaw


Tüm bağlantılar için +1! Ben de güç yasasını düşündüm, ama basit bir form (y = A t ^ k) B sabiti (y = (A t + B) ^ - 1) nedeniyle bulduğum formu içermiyor. Daha genel bir form var mı?
Erel Segal-Halevi

Eğrinin şeklini tanımlamakla ilgileniyorsanız, bir güç yasasına uymadan önce çarpanlara ayırıp kaydırmalısınız. B'ye sahip olmanız eğrinin şekliyle ilgili değildir.
dranxo

Özür dilerim, seni anlamadım, "o zaman faktör ve değişmelisin" derken ne demek istiyorsun?
Erel Segal-Halevi

X = t + B / A olarak ayarlayın. Sonra (At + B) ^ {- 1} = (A * x) ^ {- 1}, linklerdeki form.
dranxo

1
tuvalu.santafe.edu/~aaronc/courses/7000/csci7000-001_2011_L3.pdf
dranxo

7

Sorunuzla ilgili birkaç düşünce:

  • Model uyumunuzu nasıl rapor ettiğiniz, kitlenize ve alanınıza bağlı olacaktır. Örneğin, benim alanımda, R ^ 2 gibi model uyum istatistikleri çok nadiren bildirilir - ne etkileyici ne de özellikle yararlı olarak kabul edilir. Bunun yerine, geldiğiniz modele nasıl ulaştığınıza ilişkin bazı kriterler açıklanma eğilimindedir ve daha sonra model sonuçlarınızı rapor edersiniz - hepimiz modele gerçekten doğru şekilde uyduğunuzu varsayıyoruz.
  • "Ben bu formda oldum" kötü bir açıklama. Gerçekten kötü biri. Penisilin veya kinin keşfi gibi kazara deha hikayelerine olan düşkünlüğüne rağmen, "kör aptal şans" güvenilir bir bilimsel süreç değildir. Örneğin, bu formun verilerinizi sığdırma konusunda iyi olduğunu gösterdiniz, ancak henüz verilerinizi sığdırma konusunda en iyi olduğunu göstermediniz . Tek başına R ^ 2, modelinizin verilere ne kadar uyduğunu değerlendirmek için yeterli bir ölçüm değildir. Bkz Anscombe en dörtlüsü .
  • @ Rcompton'un belirttiği gibi, bilmeden bir güç yasası dağıtımına uymaya çalışıyorsunuz gibi görünüyor, ancak bir güç yasasına iyi uymayı başarsanız bile, bunun bir güç yasası olduğunu düşündüğünüz bir neden bulursanız gerçekten en iyisi . Zaman içinde Y'yi çizmek, CrossValidated'a (veya istatistiklerle daha rahat bir kolej / bölüme) gitmek ve kabaca bu görünümü verebilecek dağıtımları sistematik olarak gözden geçirmek yeterli olabilir. Güç yasası dağıtımının yanı sıra size daha iyi uyum sağlayabilecek başkaları da var.

Analizler için +1. "kabaca bu görünümü verebilecek dağıtımları sistematik olarak gözden geçirin." - bunları nerede bulabilirim?
Erel Segal-Halevi

@ErelSegalHalevi CrossValidated'da başlayabilirsiniz, bu siteler kardeş-site istatistikleri ve veri analizi ile ilgilidir.
Fomite
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.