Statystyki i duże zbiory danych r

1

Jak interpretować współczynniki z regresji beta?

Mam pewne dane, które są ograniczone od 0 do 1. Użyłem betaregpakietu w R, aby dopasować model regresji z ograniczonymi danymi jako zmienną zależną. Moje pytanie brzmi: jak interpretować współczynniki z regresji?

15 r regression interpretation beta-distribution regression-coefficients

5

Generuj normalnie rozmieszczone liczby losowe z nieokreśloną dodatnio macierzą kowariancji

Oszacowałem próbkę macierzy kowariancji próbki i otrzymałem macierz symetryczną. Z C , to proszę utworzyć n -variate normalnego rozproszonego rn a zatem potrzebny jest rozkład Cholesky'iego z C . Co powinienem zrobić, jeśli C nie jest pozytywnie określony?CCCCCCnnnCCCCCC

15 r random-generation covariance-matrix multivariate-normal cholesky

1

Prognozy szeregów czasowych z danymi dziennymi: ARIMA z regresorem

Korzystam z dziennych szeregów czasowych danych o sprzedaży, które zawierają około 2 lata codziennych punktów danych. Na podstawie niektórych samouczków / przykładów online próbowałem zidentyfikować sezonowość danych. Wydaje się, że istnieje cotygodniowa, miesięczna i prawdopodobnie roczna okresowość / sezonowość. Na przykład są dni wypłaty, szczególnie w przypadku efektu pierwszego dnia …

15 r time-series forecasting arima multiple-seasonalities

1

Jak ggplot oblicza przedziały ufności dla regresji?

Pakiet kreślarski R ggplot2 ma niesamowitą funkcję o nazwie stat_smooth do kreślenia linii regresji (lub krzywej) z powiązanym pasmem ufności. Trudno mi jednak dokładnie ustalić, w jaki sposób generowany jest ten przedział ufności, dla każdej linii regresji (lub „metody”). Jak mogę znaleźć te informacje?

15 r regression confidence-interval ggplot2

1

Jak oszacować proces Poissona za pomocą R? (Lub: jak korzystać z pakietu NHPoisson?)

Mam bazę danych zdarzeń (tj. Zmienną dat) i powiązanych zmiennych towarzyszących. Zdarzenia są generowane przez niestacjonarny proces Poissona z parametrem będącym nieznaną (ale być może liniową) funkcją niektórych zmiennych towarzyszących. Myślę, że pakiet NHPoisson istnieje właśnie w tym celu; ale po 15 godzinach nieudanych badań wciąż nie jestem w pobliżu, …

15 r poisson-distribution poisson-process

1

Znaczenie terminów wyjściowych w pakiecie GBM?

Korzystam z pakietu GBM do klasyfikacji. Zgodnie z oczekiwaniami wyniki są dobre. Ale staram się zrozumieć wyniki klasyfikatora. W produkcji występuje pięć terminów. `Iter TrainDeviance ValidDeviance StepSize Improve` Czy ktoś mógłby wyjaśnić znaczenie każdego terminu, zwłaszcza znaczenie Ulepszenia .

15 r machine-learning

1

Wizualizacja wyników modeli mieszanych

Jednym z problemów, które zawsze miałem z modelami mieszanymi, jest wymyślanie wizualizacji danych - takich, które mogłyby skończyć się na papierze lub plakacie - gdy tylko uzyska się wyniki. Obecnie pracuję nad modelem efektów mieszanych Poissona z formułą, która wygląda mniej więcej tak: a <- glmer(counts ~ X + Y …

15 r data-visualization mixed-model

5

Biblioteka Java typu open source do statystyk na poziomie oferowanym przez magisterski kurs statystyczny

Biorę udział w kursie dla absolwentów statystyki stosowanej, który wykorzystuje następujący podręcznik (aby poczuć poziom badanego materiału): koncepcje i metody statystyczne , prowadzone przez GK Bhattacharyya i RA Johnsona. Profesor wymaga od nas korzystania z SAS do odrabiania lekcji. Moje pytanie brzmi: czy istnieje biblioteka (y) Java, której można używać …

15 r sas java

3

Jakie są opcje w modelu proporcjonalnej regresji hazardu, gdy resztki Schoenfelda nie są dobre?

Robię regresję proporcjonalną hazardu Coxa przy użyciu R coxph, która obejmuje wiele zmiennych. Resztki Martingale wyglądają świetnie, a reszty Schoenfelda są świetne dla PRAWIE wszystkich zmiennych. Istnieją trzy zmienne, których reszty Schoenfelda nie są płaskie, a charakter zmiennych jest taki, że ma sens, że mogą się zmieniać w czasie. Są …

15 r model cox-model hazard schoenfeld-residuals

1

Dlaczego szacunki współczynnika regresji rlm () różnią się od lm () w R?

Korzystam z rlm w pakiecie R MASS do regresji wielowymiarowego modelu liniowego. Działa dobrze dla wielu próbek, ale otrzymuję quasi-zerowe współczynniki dla konkretnego modelu: Call: rlm(formula = Y ~ X1 + X2 + X3 + X4, data = mymodel, maxit = 50, na.action = na.omit) Residuals: Min 1Q Median 3Q …

15 r multiple-regression robust

2

Wygeneruj trzy skorelowane równomiernie rozmieszczone zmienne losowe

Załóżmy, że mamy X 2 ∼ unif ( n , 0 , 1 ) ,X1∼unif(n,0,1),X1∼unif(n,0,1),X_1 \sim \textrm{unif}(n,0,1), X2∼unif(n,0,1),X2∼unif(n,0,1),X_2 \sim \textrm{unif}(n,0,1), gdzie jest jednorodną losową próbką o wielkości n, iunif(n,0,1)unif(n,0,1)\textrm{unif}(n,0,1) Y=X1,Y=X1,Y=X_1, Z=0.4X1+1−0.4−−−−−−√X2.Z=0.4X1+1−0.4X2.Z = 0.4 X_1 + \sqrt{1 - 0.4}X_2. Zatem korelacja między i Z wynosi 0,4 .YYYZZZ0.40.40.4 Jak mogę to rozszerzyć …

15 r correlation random-generation uniform

2

Jak przetestować efekt zmiennej grupującej za pomocą modelu nieliniowego?

Mam pytanie dotyczące zastosowania zmiennej grupującej w modelu nieliniowym. Ponieważ funkcja nls () nie zezwala na zmienne czynnikowe, starałam się dowiedzieć, czy można przetestować wpływ czynnika na dopasowanie modelu. Poniżej zamieściłem przykład, w którym chcę dopasować „sezonowy sezon von Bertalanffy'ego” do różnych metod leczenia wzrostu (najczęściej stosowanych do wzrostu ryb). …

15 r mixed-model nls

3

Najlepszy sposób na wizualizację ścierania za pomocą R?

Na tej stronie niedawno odkryłem diagramy Sankeya, świetny sposób na wizualizację tego, co dzieje się w tradycyjnym schemacie blokowym. Oto dobry przykład diagramu Sankeya autorstwa George'a M. Whitesidesa i George'a W. Crabtree , Source; Nie zapomnij o długoterminowych podstawowych badaniach w dziedzinie energii , nauki 9 lutego 2007 r .: …

15 r data-visualization sankey-diagram

3

Jak rozwinąć ramkę danych w R.

Zablokowana . To pytanie i odpowiedzi są zablokowane, ponieważ pytanie jest nie na temat, ale ma znaczenie historyczne. Obecnie nie akceptuje nowych odpowiedzi ani interakcji. Mam problem z analizą R. Mam ramkę danych taką jak ta: Name | Group | Count Person 1 | A | 3 Person 2 | …

15 r

1

Jak uzyskać R-kwadrat dla lepszego dopasowania?

Jak obliczyć statystyki R-kwadrat ( ) w R dla i / lub wyniku funkcji? Na przykład dla tych danych:r2r2r^2loesspredict cars.lo <- loess(dist ~ speed, cars) cars.lp <- predict(cars.lo, data.frame(speed = seq(5, 30, 1)), se = TRUE) cars.lpma dwie tablice fitdla modelu i se.fitstandardowego błędu.

15 r r-squared loess

Pytania otagowane jako r