Jeśli masz macierz z n rzędami im kolumnami, możesz użyć SVD lub innych metod, aby obliczyć przybliżenie niskiej macierzy dla podanej macierzy. Jednak przybliżenie niskiego rzędu nadal będzie miało n wierszy i m kolumn. W jaki sposób przybliżenia niskiego poziomu mogą być przydatne w uczeniu...
Załóżmy że są ciągłymi zmiennymi losowymi ze skończonymi sekundami. wersję współczynnika korelacji rang Spearmana można zdefiniować jako współczynnik iloczynu iloczynu Pearsona ρ całek prawdopodobieństwa przekształca F_X (X) i F_Y (Y) , gdzie F_X, F_Y są cdf dla X i Y , tj.ρ s F X ( X ) F Y ( Y ) F...
Powiedzmy, że mam zestaw danych „cewnika nerkowego”. Próbuję modelować krzywą przeżycia za pomocą modelu Coxa. Jeśli wezmę pod uwagę model Coxa: potrzebuję oszacowania podstawowego zagrożenia. Korzystając z wbudowanej funkcji pakietu R , mogę łatwo to zrobić w następujący sposób:h ( t , Z) =...
Używam R i analizowałem moje dane za pomocą GLM z łączem dwumianowym. Chcę wiedzieć, jakie jest znaczenie przecięcia w tabeli wyjściowej. Punkt przecięcia dla jednego z moich modeli jest znacząco inny, jednak zmienna tak nie jest. Co to znaczy? Co to jest przechwycenie. Nie wiem, czy po prostu...
Wikipedia i winieta pakietu warstwowego R dostarczają dobrych informacji na temat założeń wspierających standardowe błędy współczynnika OLS oraz matematycznego tła estymatorów warstwowych. Nadal nie jestem pewien, w jaki sposób rozwiązany jest problem heteroscedastyczności resztek, prawdopodobnie...
Czy ktoś może udowodnić następujący związek między wskaźnikiem informacji Fishera a względną entropią (lub dywergencją KL) w czysto matematyczny, rygorystyczny sposób? D(p(⋅,a+da)∥p(⋅,a))=12gi,jdaidaj+(O(∥da∥3)D(p(⋅,a+da)∥p(⋅,a))=12gi,jdaidaj+(O(‖da‖3)D( p(\cdot , a+da) \parallel p(\cdot,a) )...
Chcę, aby zaimplementować algorytm EM ręcznie, a następnie porównać je do wyników działań normalmixEMz mixtoolsopakowania. Oczywiście byłbym szczęśliwy, gdyby oba doprowadziły do tych samych rezultatów. Głównym odniesieniem jest Geoffrey McLachlan (2000), Finite Mixture Models . Mam gęstość...
Czy caretpakiet R sprawdza poprawność krzyżową dla obu modeli alphai lambdadla glmnetmodelu? Uruchamianie tego kodu eGrid <- expand.grid(.alpha = (1:10) * 0.1, .lambda = (1:10) * 0.1) Control <- trainControl(method = "repeatedcv",repeats = 3,verboseIter =TRUE) netFit <- train(x...
Jak mogę próbkować z rozkładu mieszaniny, a w szczególności mieszaniny rozkładów normalnych w R? Na przykład, jeśli chcę próbkować z: 0.3×N(0,1)+0.5×N(10,1)+0.2×N(3,.1)0.3×N(0,1)+0.5×N(10,1)+0.2×N(3,.1) 0.3\!\times\mathcal{N}(0,1)\; + \;0.5\!\times\mathcal{N}(10,1)\; +...
Uczę się R i eksperymentowałem z analizą wariancji. Prowadziłem oba kruskal.test(depVar ~ indepVar, data=df) i anova(lm(depVar ~ indepVar, data=dF)) Czy istnieje praktyczna różnica między tymi dwoma testami? Rozumiem, że oboje oceniają hipotezę zerową, że populacje mają ten sam środek....
Cześć, próbuję znaleźć nieparametryczny odpowiednik dwukierunkowej ANOVA (konstrukcja 3x4), która może zawierać interakcje. Z mojego czytania z Zar 1984 „Analiza biostatystyczna” jest to możliwe przy użyciu metody przedstawionej w Scheirer, Ray i Hare (1976), jednak zgodnie z innymi postami online...
W tym modelu regresji liniowej jednowymiarowej yja= β0+ β1xja+ ϵjayja=β0+β1xja+ϵjay_i = \beta_0 + \beta_1x_i+\epsilon_i podane zestaw danych D = { ( x1, y1) , . . . , ( xn, yn) }re={(x1,y1),...,(xn,yn)}D=\{(x_1,y_1),...,(x_n,y_n)\} , oszacowania współczynników β 1 = Σ i x i y i - n ˉ x ˉ Y β 0=...
Rozważ obserwacje dobrze ocenzurowane z wydarzeniami czasami . Liczba podatnych osobników w czasie wynosi , a liczba zdarzeń w czasie wynosi .t1, t2), …t1,t2,…t_1, t_2, \dotsn I I d Ijaiinjanin_ijaiirejadid_i Kaplan-Meier lub estymator produktu powstaje naturalnie jako MLE, gdy funkcja przeżycia...
Usiłuję zrozumieć wyprowadzenie oczekiwanego błędu prognozy na niższy poziom (ESL), szczególnie na podstawie wyprowadzenia 2.11 i 2.12 (warunkowanie, krok w kierunku minimum punktowego). Wszelkie wskazówki lub linki są mile widziane. Poniżej raportuję fragment z ESL str. 18. Pierwsze dwa równania...
Mam zestaw danych składa się z 5 funkcji: A, B, C, D, E. Wszystkie są wartości liczbowe. Zamiast tworzyć klastrowanie oparte na gęstości, chcę skupić dane w sposób podobny do drzewa decyzyjnego. Mam na myśli takie podejście: Algorytm może dzielić dane na X początkowych klastrów w oparciu o cechę...
Niech model regresji liniowej uzyskany przez funkcję R lm chciałby wiedzieć, czy można to uzyskać za pomocą polecenia Mean Squared Error. Miałem następujący wynik w przykładzie > lm <- lm(MuscleMAss~Age,data) > sm<-summary(lm) > sm Call: lm(formula = MuscleMAss ~ Age, data =...
W R, kiedy mam (uogólniony) model liniowy ( lm, glm, gls, glmm, ...), jak mogę sprawdzić współczynnik (nachylenie regresji) przed jakąkolwiek inną wartość niż 0? W podsumowaniu modelu wyniki testu t współczynnika są automatycznie raportowane, ale tylko dla porównania z 0. Chcę porównać go z inną...
Interesuje mnie lepsze zrozumienie metody delta do aproksymacji standardowych błędów średnich efektów krańcowych modelu regresji, który obejmuje termin interakcji. Patrzyłem na powiązane pytania metodą delta, ale żadne nie dostarczyło tego, czego szukam. Rozważ następujące przykładowe dane jako...
Podczas gdy te dwa wszechobecne terminy są często używane jako synonimy, czasem wydaje się, że istnieje rozróżnienie. Czy rzeczywiście jest różnica, czy są one dokładnie
Wiem, że przeprowadzanie strojenia hiperparametrów poza walidacją krzyżową może prowadzić do stronniczo wysokich oszacowań zewnętrznej ważności, ponieważ zestaw danych używany do mierzenia wydajności jest taki sam, jak użyty do dostrojenia funkcji. Zastanawiam się, jak poważny jest to problem ....