Dopasowałem model z kilkoma zmiennymi niezależnymi, z których jedną jest opóźnienie zmiennej zależnej, używając pakietu dynlm. Zakładając, że mam prognozy o krok do przodu dla moich zmiennych niezależnych, w jaki sposób mogę uzyskać prognozy o krok do przodu dla moich zmiennych zależnych? Oto przykład: library(dynlm) y<-arima.sim(model=list(ar=c(.9)),n=10) #Create AR(1) dependant variable …
Oto, co chcę zrobić, ale wydaje się, że nie ma predictmetody dla mlogit. Jakieś pomysły? library(mlogit) data("Fishing", package = "mlogit") Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode") Fish_fit<-Fish[-1,] Fish_test<-Fish[1,] m <- mlogit(mode ~price+ catch | income, data = Fish_fit) predict(m,newdata=Fish_test)
Przeprowadziłem test Levene'a i Bartletta na grupach danych z jednego z moich eksperymentów, aby potwierdzić, że nie naruszam założenia ANOVA o jednorodności wariancji. Chciałbym sprawdzić z wami, że nie robię żadnych niewłaściwych założeń, jeśli nie macie nic przeciwko: D Wartość p zwrócona przez oba te testy to prawdopodobieństwo, że moje …
Powiedzmy, że mam dane: x1 <- rnorm(100,2,10) x2 <- rnorm(100,2,10) y <- x1+x2+x1*x2+rnorm(100,1,2) dat <- data.frame(y=y,x1=x1,x2=x2) res <- lm(y~x1*x2,data=dat) summary(res) Chcę wykreślić ciągłe przez ciągłą interakcję tak, że x1 jest na osi X, a x2 jest reprezentowany przez 3 linie, jedna reprezentująca x2 przy wyniku Z wynoszącym 0, jedna przy …
Poszukuję, jak (wizualnie) wyjaśnić studentom pierwszego roku prostą korelację liniową. Klasycznym sposobem wizualizacji byłoby stworzenie wykresu rozproszenia Y ~ X z prostą linią regresji. Ostatnio wpadłem na pomysł rozszerzenia tego typu grafiki, dodając do wykresu 3 kolejne obrazy, pozostawiając mi: wykres rozproszenia y ~ 1, następnie y ~ x, res …
Każdy, kto podąży za baseballem, prawdopodobnie słyszał o nieoczekiwanym występie Jose Bautisty w Toronto typu MVP. W ciągu czterech poprzednich lat osiągnął około 15 przebiegów u siebie w sezonie. W zeszłym roku osiągnął 54 lata, a liczba ta przekroczyła zaledwie 12 graczy w historii baseballu. W 2010 roku wypłacono mu …
Natknąłem się na zdjęcie prototypu aplikacji, który znajduje znaczące zmiany („trendy” - nie wzrosty / wartości odstające) w danych o ruchu: Chcę napisać program (Java, opcjonalnie R), który jest w stanie zrobić to samo - ale ponieważ moje umiejętności statystyczne są nieco zardzewiałe, muszę ponownie zagłębić się w ten temat. …
Jestem początkującym do statystyki (tylko jeden kurs college'u), ale mam doświadczenie w programowaniu. Właśnie zacząłem bawić się biblioteką klasyfikatora Bayesa dla Ruby i szukam pomysłów na rzeczy do analizy. W tej chwili mam do czynienia z kategoryzacją Tweet, ale czy masz jakieś pomysły? Co ważniejsze, w jaki sposób mogę dowiedzieć …
Zakładając, że istnieje sens w testowaniu założenia normalności dla anova (patrz 1 i 2 ) Jak można to przetestować w R? Spodziewałbym się zrobić coś takiego: ## From Venables and Ripley (2002) p.165. utils::data(npk, package="MASS") npk.aovE <- aov(yield ~ N*P*K + Error(block), npk) residuals(npk.aovE) qqnorm(residuals(npk.aov)) Co nie działa, ponieważ „reszty” …
Zmienne zależne w MANOVA nie powinny być „zbyt silnie skorelowane”. Ale jak silna korelacja jest zbyt silna? Interesujące byłoby uzyskanie opinii ludzi na ten temat. Na przykład, czy kontynuowałbyś MANOVA w następujących sytuacjach? Y1 i Y2 są skorelowane z ip < 0,005r = 0,3r=0.3r=0.3p < 0,005p<0.005p<0.005 Y1 i Y2 są …
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 w zeszłym roku . Jakie implementacje typu open source - w jakimkolwiek języku - istnieją, które mogą obliczać ścieżki regulacyjne lasso dla regresji liniowej poprzez …
Dotyczy to tylko mojego projektu dla licealistów, więc nie musi być idealny. Realizuję projekt dotyczący globalnego ocieplenia i chcę badać ludzi pod kątem ich opinii. Wiem, że jeśli użyję wygodnej próbki moich kolegów z klasy, będę miał dużo uprzedzeń. Zastanawiałem się, czy istnieje strona internetowa, w której mogłabym opublikować ankietę …
Chcę zaimplementować (w R) następujący bardzo prosty dynamiczny model liniowy, dla którego mam 2 nieznane parametry zmieniające się w czasie (wariancja błędu obserwacji i wariancja błędu stanu ). ϵ 2 tϵ1tϵt1\epsilon^1_tϵ2)tϵt2\epsilon^2_t Ytθt + 1==θt+ ϵ1tθt+ ϵ2)tYt=θt+ϵt1θt+1=θt+ϵt2 \begin{matrix} Y_t & = & \theta_t + \epsilon^1_t\\ \theta_{t+1} & = & \theta_{t}+\epsilon^2_t \end{matrix} …
Jeśli zdecydujesz się przeanalizować projekt kontrolny przed leczeniem z ciągłą zmienną zależną przy użyciu mieszanej ANOVA, istnieją różne sposoby oceny wpływu przebywania w grupie leczenia. Efekt interakcji jest jedną z głównych opcji. Ogólnie, szczególnie podoba mi się miara typu Cohena d (tj. ). Nie podoba mi się wyjaśnienie miary wariancji, …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.