«computational-statistics» 태그된 질문

통계 및 컴퓨팅의 인터페이스를 나타냅니다. 통계 목적으로 알고리즘 및 소프트웨어 사용.

1
데이터 행렬이 대각선 일 때 올가미 문제에 대한 닫힌 양식 솔루션
\newcommand{\diag}{\operatorname{diag}} 문제가 있습니다 : 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), : ∑i=1nxixTi=diag(σ21,...,σ2d).∑i=1nxixiT=diag⁡(σ12,...,σd2).\sum_{i=1}^nx_ix_i^T=\diag(\sigma_1^2,...,\sigma_d^2). 이 경우 폐쇄 형 솔루션이 있습니까? 나는 그것을 가지고있다 : (XTX)−1=diag(σ−21,...,σ−2d),(XTX)−1=diag⁡(σ1−2,...,σd−2),(X^TX)^{-1}=\diag\left(\sigma_1^{-2},...,\sigma_d^{-2}\right), 그래서 대답은 : 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\}, 에 대한 yj=∑i=1nyixijσ2iyj=∑i=1nyixijσi2y\,^j=\displaystyle\sum_{i=1}^n\frac{y_ix_i\,^j}{\sigma_i^2} 하지만 확실하지 않습니다.

7
통계 이론 및 응용 프로그램에서 이해하기
저는 최근 공학 및 수학을 배경으로 의학 및 생물학 모델링 석사 학위를 취득했습니다. 교육 프로그램에 상당히 많은 수의 수학 통계 (아래 목록 참조) 과정이 포함되어 있었지만 통계와 이론의 적용에 대해 완전히 잃어버린 경우가 종종 있습니다. "순수한"수학과 비교할 때 통계는 실제로 의미가 없습니다. 특히 대부분의 통계 학자 (내 과거 강사를 포함)가 …

1
복잡한 모델을 대규모 데이터 세트에 반복적으로 피팅 할 때 계산 효율성을 어떻게 최적화 할 수 있습니까?
MCMCglmm혼합 효과 모델을 실행하기 위해 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) 데이터에는 약 20,000 개의 관측치가 있으며 약 200 개의 학교에 모여 있습니다. 실행하기 전에 데이터 프레임에서 사용하지 않는 모든 변수를 삭제하고 메모리에서 다른 모든 개체를 …

4
시계열에서의 이상치 탐지 : 오 탐지를 줄이는 방법은 무엇입니까?
시계열에서 이상치 탐지를 자동화하려고하는데 여기에서 Rob Hyndman이 제안한 솔루션 수정을 사용했습니다 . 여러 국가의 웹 사이트 방문을 매일 측정합니다. 일일 방문이 몇 번 또는 몇 천 번인 일부 국가의 경우 내 방법이 합리적으로 작동하는 것 같습니다. 그러나 국가가 하루에 1 ~ 2 회 방문으로 이어지는 경우 알고리즘의 한계가 매우 좁아서 …

3
컴퓨터 시뮬레이션을 사용하여 대학원 수준의 통계 개념을 더 잘 이해
안녕하세요, 저는 통계학 대학원 과정을 수강하고 있으며 테스트 통계 및 기타 개념을 다루었습니다. 그러나 나는 종종 공식을 적용하고 물건이 어떻게 작동하는지에 대한 일종의 직관을 개발할 수 있지만, 아마도 모의 실험으로 연구를 백업하면 당면한 문제에 대해 더 나은 직관을 개발할 것이라는 느낌이 종종 듭니다. . 그래서, 우리가 수업에서 논의하는 일부 개념을 …

1
이 "최대 상관 계수"는 무엇입니까?
일반적인 이미지 처리 통계는 Haralick 텍스처 기능 ( 14)을 사용하는 것입니다. 이러한 기능 중 14 번째 기능에 대해 궁금합니다. 인접 맵 (두 정수 의 경험적 분포를 간단히 볼 수 있음)를 고려하면 다음과 같이 정의됩니다. 의 두 번째 고유 값의 제곱근 , 여기서, 이다 :피피Pi , j &lt; 256나는,제이&lt;256i,j < 256큐큐Q큐큐Q …

1
치료 그룹에 필적하는 대조 그룹을 찾으십니까?
수학 보조 소프트웨어를 사용하는 30 명 (캘리포니아의 30 개 학교)의 치료 그룹이 있습니다. 간단한 분석에서, 우리의 치료 그룹과 비슷한 대조군 사이의 학생들의 평균 수학 성장을 비교하고 싶습니다. CA에는 소프트웨어를 사용하지 않은 많은 학교가 있습니다. 통제 그룹에 유사한 수행 학교를 포함 시키길 원합니다 (그들의 기준 점수는 합리적인 오차 범위를 가진 치료 …


3
R에서 (또는 일반적으로) 회귀 계수를 특정 부호로 만들 수 있습니까?
실제 데이터를 사용하고 있으며 회귀 모델이 반 직관적 인 결과를 산출합니다. 일반적으로 나는 통계를 신뢰하지만 실제로 이러한 것들 중 일부는 사실이 될 수 없습니다. 내가보고있는 주요 문제는 실제로 하나의 변수가 증가하면 실제로는 음의 상관 관계가 있어야 응답이 증가한다는 것입니다. 각 회귀 계수에 대해 특정 부호를 강제하는 방법이 있습니까? 이 작업을 …


4
통계 소프트웨어 테스트
통계 소프트웨어를 테스트하는 데 유용한 기술 / 접근법은 무엇입니까? 특히 최대 가능성을 사용하여 파라 메트릭 추정을 수행하는 프로그램에 관심이 있습니다. 내 자신의 프로그램을 작성할 때 대부분의 경우 필요한 계산이 기존 시스템에서 아직 구현되지 않았기 때문에 다른 프로그램이나 게시 된 소스의 결과를 항상 비교할 수있는 것은 아닙니다. 정확성을 보장 할 수있는 …

4
경사 하강이 필요한 이유는 무엇입니까?
비용 함수를 차별화하고 모든 매개 변수에 대해 부분 미분을 통해 얻은 방정식을 풀어 매개 변수를 찾고 비용 함수가 최소 인 위치를 찾을 수 있습니다. 또한 파생 상품이 0 인 여러 장소를 찾을 수 있다고 생각합니다. 따라서 모든 장소를 확인할 수 있고 글로벌 최소값을 찾을 수 있습니다 왜 경사 하강이 대신 …

2
음수가 아닌 정수에 대한 이산 분포에서 표본 추출하는 방법은 무엇입니까?
는 상수로 알려진 다음과 같은 불연속 분포 가 있습니다.α , βα,β\alpha,\beta p ( x ; α , β) =베타 ( α + 1 , β+ x )베타 ( α , β)위한 X = 0 , 1 , 2 , ...p(x;α,β)=Beta(α+1,β+x)Beta(α,β)for 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 …

3
파이썬을 이용한 시계열 이상 탐지
여러 시계열 데이터 세트에서 이상 감지를 구현해야합니다. 나는 전에 이것을 한 적이 없으며 조언을 기대하고있었습니다. 나는 파이썬에 매우 익숙하므로 솔루션을 구현하는 것을 선호합니다 (대부분의 코드는 내 작업의 다른 부분에 대한 파이썬입니다). 데이터 설명 : 지난 2 년 정도 전에 수집 된 월별 시계열 데이터입니다 (즉, 24-36 시간 만). 기본적으로 여러 …

1
샘플링 비용
다음과 같은 시뮬레이션 문제가 발생했습니다. {ω1,…,ωd}{ω1,…,ωd}\{\omega_1,\ldots,\omega_d\} 알려진 실수의 분포 {−1,1}d{−1,1}d\{-1,1\}^d 에 의해 정의 P(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)_+ 의 긍정적 인 부분을 나타냅니다 zzz. 이 분포를 대상으로하는 Metropolis-Hastings 샘플러는 생각할 수 있지만 알고리즘의 순서를 줄이기 위해 많은 수의 제로 확률을 사용하는 효율적인 직접 샘플러가 있는지 궁금합니다.O(2d)O(2d)O(2^d) 에 O(d)O(d)O(d).

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.