Statystyki i duże zbiory danych r

2

Mam zestaw danych miesięcznych dotyczących temperatury powierzchni morza (SST) i chcę zastosować metodologię klastrową do wykrywania regionów o podobnych wzorcach SST. Mam zestaw miesięcznych plików danych od 1985 do 2009 roku i chcę zastosować klastrowanie do każdego miesiąca jako pierwszy krok. Każdy plik zawiera dane siatkowe dla 358416 punktów, z …

12 r clustering spatial

1

Testuj równoważność modeli zagnieżdżonych

Powiedzmy, że jest funkcją liniową i manekina d . Moja hipoteza jest taka, że d sama jest jak hedonistyczny indeks wektora innych zmiennych, Z . Mam poparcie dla tego podejścia w MANOVA z Z (tj z_1 , z_2 , ..., z_n ) na d . Czy istnieje sposób przetestowania równoważności …

12 r hypothesis-testing regression model-selection

3

Właściwy sposób radzenia sobie z 3-poziomową tabelą awaryjną

Mam trzypoziomową tabelę zdarzeń awaryjnych, z danymi zliczeń dla kilku gatunków, rośliny żywicielskiej, z której zostały zebrane, oraz tego, czy zbiór ten miał miejsce w deszczowy dzień (to naprawdę ma znaczenie!). Przy użyciu R fałszywe dane mogą wyglądać mniej więcej tak: count <- rpois(8, 10) species <- rep(c("a", "b"), 4) …

12 r categorical-data log-linear

7

Co to są pakiety do eksploracji tekstu dla R i czy istnieją inne programy do eksploracji tekstu typu open source?

Czy możesz polecić pakiet do eksploracji tekstu w języku R, który może być używany w przypadku dużych ilości danych? Po drugie, czy jest dostępny GUI dla któregokolwiek z pakietów eksploracji tekstu w R? Po trzecie, czy istnieje inny program do eksploracji tekstu typu open source, który jest łatwy i intuicyjny …

12 r text-mining

4

Analiza danych wiatru za pomocą R.

Cześć, analizuję dane dotyczące wiatru w celu oszacowania energii z turbiny wiatrowej. Wziąłem 10 lat danych wiatrowych i wykreśliłem histogram; moim drugim etapem było dopasowanie rozkładu Weibulla do danych. Użyłem R z pakietem lmomdo obliczenia kształtu i skali Weibula. Użyłem tego kodu: >library(lmom) wind.moments<-samlmu(as.numeric(pp$WS)) moments<-pelwei(wind.moments) x.wei<-rweibull(n=length(pp$WS), shape=moments["delta"], scale=moments["beta"]) hist(as.numeric(pp$WS), freq=FALSE) …

12 r distributions

4

Konfiguracja Sweave, R, Latex, Eclipse StatET [zamknięte]

Zamknięte. To pytanie jest nie na temat . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby było tematem dotyczącym weryfikacji krzyżowej. Zamknięte 3 lata temu . Kilka dni temu widziałem post, jak skonfigurować SweaveR, który pozwoliłby użytkownikowi na bezpośrednie eksportowanie rzeczy takich jak tabele, wykresy itp. Do …

12 r

2

ANOVA Split-Plot: testy porównawcze modeli w R.

Jak mogę testować efekty w analizie ANOVA typu split-plot, używając odpowiednich porównań modeli do użycia z argumentami Xi w R? Znam i Dalgaard (2007) [1]. Niestety szczotkuje tylko projekty Split-Plot. Robi to w całkowicie losowy sposób z dwoma czynnikami wewnątrz badanych:Manova.mlm()?anova.mlm N <- 20 # 20 subjects total P <- …

12 r anova multivariate-analysis repeated-measures split-plot

1

Filtrowanie ramki danych

Wciąż uczysz się podstawowych funkcji w R, funkcja podzbioru wydaje się filtrować tylko warunek oparty na pojedynczej kolumnie z wieloma warunkami lub bez? Jak mogę łatwo filtrować dane z ramki danych? gdy otrzymasz wiele warunków Kiedy warunek należy zastosować w dostępnych kolumnach. Przykład: biorąc pod uwagę ramkę danych zawierającą name …

12 r

2

Jak wykonać test hipotezy MCMC na modelu regresji z efektem mieszanym z losowymi nachyleniami?

Biblioteka languageR zapewnia metodę (pvals.fnc) do testowania istotności MCMC ustalonych efektów w modelu dopasowania regresji mieszanej przy użyciu lmera. Jednak pvals.fnc podaje błąd, gdy model Lmer zawiera losowe zbocza. Czy istnieje sposób przeprowadzenia testu hipotetycznego MCMC dla takich modeli? Jeśli tak to jak? (Aby zostać zaakceptowanym, odpowiedź powinna zawierać sprawdzony …

12 r mixed-model statistical-significance monte-carlo

5

Lepsza klasyfikacja domyślna w regresji logistycznej

Pełne ujawnienie: To zadanie domowe. Zamieściłem link do zestawu danych ( http://www.bertelsen.ca/R/logistic-regression.sav ) Moim celem jest zmaksymalizowanie prognozy osób spłacających zaległości kredytowe w tym zbiorze danych. Każdy model, który do tej pory wymyśliłem, przewiduje> 90% domyślnych, ale <40% domyślnych, co daje ogólną skuteczność klasyfikacji ~ 80%. Zastanawiam się więc, czy …

12 r logistic spss self-study

2

Analiza przedmiotu dla początkującego R.

Próbuję ocenić 20-elementowy test wielokrotnego wyboru. Chcę przeprowadzić analizę przedmiotu, taką jak w tym przykładzie . Tak więc dla każdego pytania chcę wartość P i korelację z sumą oraz rozkład wybranych opcji. Nie wiem nic o różnych pakietach oprogramowania statystycznego, ale chciałbym używać R, ponieważ czuję się swobodnie w programowaniu, …

12 r correlation psychometrics scales

3

Czy istnieje poważny problem z pomijaniem obserwacji z brakującymi wartościami podczas obliczania macierzy korelacji?

Mam ten ogromny zestaw danych z około 2500 zmiennymi i podobnymi 142 obserwacjami. Chcę uruchomić korelację między zmienną X a resztą zmiennych. Ale w wielu kolumnach brakuje wpisów. Próbowałem to zrobić w R za pomocą argumentu „pairwise-complete” ( use=pairwise.complete.obs) i uzyskałem wiązkę korelacji. Ale potem ktoś na StackOverflow opublikował link …

12 r correlation missing-data correlation-matrix

1

Dlaczego lm i biglm w R dają różne wartości p dla tych samych danych?

Oto mały przykład: MyDf<-data.frame(x=c(1,2,3,4), y=c(1.2, .7, -.5, -3)) Teraz z base::lm: > lm(y~x, data=MyDf) %>% summary Call: lm(formula = y ~ x, data = MyDf) Residuals: 1 2 3 4 -0.47 0.41 0.59 -0.53 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 3.0500 0.8738 3.491 0.0732 . x -1.3800 0.3191 …

12 r regression p-value linear-model

3

Jak zaimplementować model mieszany za pomocą funkcji betareg w R?

Mam zestaw danych składający się z proporcji, które mierzą „poziom aktywności” poszczególnych kijanek, dzięki czemu wartości są powiązane od 0 do 1. Dane te zostały zebrane przez zliczenie liczby ruchów danej osoby w określonym przedziale czasu (1 dla ruchu, 0 za brak ruchu), a następnie uśrednia się, aby utworzyć jedną …

12 r mixed-model random-effects-model beta-distribution beta-regression

1

Gradient dla funkcji straty logistycznej

Zadałbym pytanie związane z tym . Znalazłem przykład pisania niestandardowej funkcji utraty dla xgboost tutaj : loglossobj <- function(preds, dtrain) { # dtrain is the internal format of the training data # We extract the labels from the training data labels <- getinfo(dtrain, "label") # We compute the 1st and …

12 r machine-learning gradient-descent boosting loss-functions

Pytania otagowane jako r