«computational-statistics» etiketlenmiş sorular

İstatistik ve hesaplama arayüzünü ifade eder; istatistiksel amaçlar için algoritma ve yazılım kullanımı.

1
Veri matrisi diyagonal olduğunda kement problemine kapalı form çözümü
Sorunumuz\newcommand{\diag}{\operatorname{diag}} var: minw∈Rd(1n∑i=1n(⟨w,xi⟩−yi)2+2λ||w||1),minw∈Rd(1n∑i=1n(⟨w,xi⟩−yi)2+2λ||w||1),\min_{w\in\mathbb{R}^{d}}\left( \frac{1}{n}\sum_{i=1}^{n} \left( \langle w,x_{i}\rangle-y_{i} \right)^{2} +2\lambda||w||_1\right), şu varsayımla: ∑i=1nxixTi=diag(σ21,...,σ2d).∑i=1nxixiT=diag⁡(σ12,...,σd2).\sum_{i=1}^nx_ix_i^T=\diag(\sigma_1^2,...,\sigma_d^2). Bu durumda kapalı bir çözüm var mı? Ben var: (XTX)−1=diag(σ−21,...,σ−2d),(XTX)−1=diag⁡(σ1−2,...,σd−2),(X^TX)^{-1}=\diag\left(\sigma_1^{-2},...,\sigma_d^{-2}\right), ve bu yüzden cevap olduğunu düşünüyorum : wj=yjmax{0,1−λn|yj|},wj=yjmax{0,1−λn|yj|},w\,^j=y\,^j\max\left\{0,1-\lambda \frac{n}{|y^j|}\right\}, için yj=∑i=1nyixijσ2iyj=∑i=1nyixijσi2y\,^j=\displaystyle\sum_{i=1}^n\frac{y_ix_i\,^j}{\sigma_i^2} , ancak emin değilim.

7
İstatistik teorisi ve uygulamalarından anlam çıkarma
Kısa bir süre önce mühendislik matematiği ile birlikte tıbbi ve biyolojik modelleme üzerine yüksek lisansımı bitirdim. Eğitim programım oldukça yüksek notlarla yönettiğim matematik istatistikleri üzerine önemli miktarda ders içermesine rağmen (bir liste için aşağıya bakın), genellikle istatistik teorisine ve uygulamalarına bakarken tamamen kayboldum. Şunu söylemeliyim ki, "saf" matematiğe kıyasla, istatistikler …

1
Büyük bir veri setine tekrar tekrar karmaşık bir model takarken hesaplama verimliliğini nasıl optimize edebilirim?
MCMCglmmKarışık efektler modeli çalıştırmak için R paketini kullanarak performans sorunları yaşıyorum . Kod şöyle görünür: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn, verbose=F) Verilerde yaklaşık 20.000 gözlem vardır ve bunlar yaklaşık 200 okulda kümelenmiştir. Kullanılmamış tüm değişkenleri veri çerçevesinden bıraktım ve çalıştırmadan önce diğer tüm nesneleri bellekten kaldırdım. …

4
Zaman Serilerinde Aykırı Tespit: Yanlış pozitifler nasıl azaltılır?
Zaman serilerinde aykırı algılamayı otomatikleştirmeye çalışıyorum ve burada Rob Hyndman tarafından önerilen çözümde bir değişiklik kullandım . Diyelim ki, çeşitli ülkelerden bir web sitesine günlük ziyaretleri ölçüyorum. Günlük ziyaretlerin birkaç hundrend veya binlerce olduğu bazı ülkelerde, yöntemim makul çalışıyor gibi görünüyor. Bununla birlikte, bir ülkenin günde sadece 1 veya 2 …

3
Lisansüstü düzeyde istatistiksel kavramları daha iyi anlamak için bilgisayar simülasyonlarını kullanma
Merhaba, İstatistik alanında lisansüstü bir ders alıyorum ve Test istatistiklerini ve diğer kavramları ele alıyoruz. Bununla birlikte, formülleri sık sık uygulayabiliyorum ve işlerin nasıl çalıştığına dair bir tür sezgi geliştirebiliyorum, ancak belki de çalışmamı simüle edilmiş deneylerle desteklersem, eldeki sorunlara daha iyi sezgi geliştireceğim hissine kapılıyorum. . Bu yüzden, sınıfta …

1
Bu “maksimum korelasyon katsayısı” nedir?
Tipik bir görüntü işleme istatistiği 14 olan Haralick doku özelliklerinin kullanılmasıdır. Bu özelliklerin 14'ünü merak ediyorum: Bir bitişik harita verildiğinde (iki tamsayı ampirik dağılımını görebiliriz ), şu şekilde tanımlanır: ikinci özdeğerinin kare kökü , burada :PPPi,j&lt;256i,j&lt;256i,j < 256QQQQQQ Qij=∑kP(i,k)P(j,k)[∑xP(x,i)][∑yP(k,y)]Qij=∑kP(i,k)P(j,k)[∑xP(x,i)][∑yP(k,y)]Q_{ij} = \sum_k \frac{ P(i,k) P(j,k)}{ [\sum_x P(x,i)] [\sum_y P(k, y)] } …

1
Bir tedavi grubu için karşılaştırılabilir bir Kontrol grubu mu arıyorsunuz?
Matematik destek yazılımı kullanan 30 kişilik bir tedavi grubum var (Kaliforniya'da 30 okul). Basit bir analizde, öğrencilerin tedavi grubumuzla karşılaştırılabilir bir kontrol grubu arasındaki ortalama Matematik büyümesini karşılaştırmak istiyorum. CA'da yazılımı kullanmayan birçok okul var. Kontrol grubunun benzer performans gösteren okulları içermesini istiyorum (temel puanları makul bir hata payı olan …


3
R'de (veya genel olarak) regresyon katsayılarını belirli bir işaret olmaya zorlamak mümkün müdür?
Bazı gerçek dünya verileri ile çalışıyorum ve regresyon modelleri bazı mantıksız sonuçlar veriyor. Normalde istatistiklere güvenirim ama gerçekte bunlardan bazıları doğru olamaz. Gördüğüm temel sorun, bir değişkendeki artışın, gerçekte, negatif olarak ilişkilendirilmeleri gerektiğinde, yanıtta bir artışa neden olmasıdır. Regresyon katsayılarının her biri için belirli bir işareti zorlamanın bir yolu var …

1
Düşük seviyeli lineer sistemin hızlı hesaplanması / tahmini
Doğrusal denklem sistemleri, hesaplama istatistiklerinde yaygındır. Karşılaştığım özel bir sistem (örneğin, faktör analizinde) sistemdir A x = bbirx=bAx=b burada Burada D , kesinlikle pozitif bir diyagonal olan bir n × n diyagonal matristir, Ω bir m × m ( m ≪ n ile ) simetrik pozitif yarı tanımlanmış matristir ve …

4
İstatistiksel yazılımı test etme
İstatistiksel yazılımın testinde hangi teknikler / yaklaşımlar yararlıdır? Özellikle maksimum olasılık kullanarak parametrik tahmin yapan programlarla ilgileniyorum. Başka programlardan veya yayınlanmış kaynaklardan elde edilen sonuçları karşılaştırmak her zaman mümkün değildir, çünkü çoğu zaman kendi başıma bir program yazarken ihtiyacım olan hesaplama zaten mevcut bir sistemde uygulanmamıştır. Doğruluğu garanti edebilecek yöntemlerde …

4
Degrade iniş neden gereklidir?
Maliyet fonksiyonunu farklılaştırabildiğimizde ve her parametreye göre kısmi farklılaşma yoluyla elde edilen denklemleri çözerek ve maliyet fonksiyonunun minimum nerede olduğunu bulabildiğimizde parametreleri bulabiliriz. Ayrıca, türevlerin sıfır olduğu birden fazla yer bulmak mümkün, böylelikle tüm bu yerleri kontrol edebilir ve küresel minima bulabiliriz. bunun yerine gradyan inişi neden yapılır?

2
Negatif olmayan tamsayılarda ayrık dağılımdan nasıl örnek alınır?
Aşağıdaki bağımsız dağıtım var, burada bilinen sabitler:α , βα,β\alpha,\beta p ( x ; α , β) =Beta ( α + 1 , β+ x )Beta ( α , β)için x = 0 , 1 , 2 , ...p(x;α,β)=Beta(α+1,β+x)Beta(α,β)için x=0,1,2,... p(x;\alpha,\beta) = \frac{\text{Beta}(\alpha+1, \beta+x)}{\text{Beta}(\alpha,\beta)} \;\;\;\;\text{for } x = 0,1,2,\dots Bu dağıtımdan …

3
Python ile Zaman Serisi Anomali Tespiti
Birkaç zaman serisi veri kümesinde anomali tespiti uygulamam gerekiyor. Bunu daha önce hiç yapmadım ve tavsiye almayı umuyordum. Ben python ile çok rahat, bu yüzden çözüm uygulanmasını tercih ederdim (benim kod çoğu işimin diğer bölümleri için python olduğunu). Verilerin açıklaması: Son 2 yılda toplanmaya başlanan aylık zaman serisi verileri (yani …

1
örnekleme maliyeti
Aşağıdaki simülasyon sorun geldi: verilen bir dizi bilinen gerçek sayılar, bir dağıtım ile tanımlanır burada , pozitif kısmını belirtir . Bu dağılımı hedefleyen bir Metropolis-Hastings örnekleyicisini düşünürken, algoritmanın sırasını den ya düşürmek için çok sayıda sıfır olasılıktan yararlanarak verimli bir doğrudan örnekleyici olup olmadığını merak ediyorum. .{ω1,…,ωd}{ω1,…,ωd}\{\omega_1,\ldots,\omega_d\}{−1,1}d{−1,1}d\{-1,1\}^dP(X=(x1,…,xd))∝(x1ω1+…+xdωd)+P(X=(x1,…,xd))∝(x1ω1+…+xdωd)+\mathbb{P}(X=(x_1,\ldots,x_d))\propto (x_1\omega_1+\ldots+x_d\omega_d)_+(z)+(z)+(z)_+zzzO(2d)O(2d)O(2^d)O(d)O(d)O(d)

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.