W hierarchicznym modelu danych którym wydaje się typowe w praktyce, aby wybierać wartości ( takie, że średnia i wariancja rozkładu gamma w przybliżeniu odpowiadają średniej i wariancji danych (np. Clayton i Kaldor, 1987 „Empirical Bayes Estimates of Standaryzated Age Relative Risks for Disease Mapping”, Biometrics ). Najwyraźniej jest to jednak …
Chcę powiązać wady kodu z miernikami złożoności kodu, takimi jak bliskość. Jednym z powszechnych modeli jest postrzeganie tego jako procesu Poissona, w którym czas trwania to ilość czasu poświęcanego na kodowanie, a gęstość jest funkcją złożoności kodu. Jestem w stanie zrobić regresję i uzyskać wartości istotności itp. Jednak trudno mi …
Mam zestaw danych zawierający liczbę działań wykonanych przez osoby w ciągu 7 dni. Konkretne działanie nie powinno być istotne dla tego pytania. Oto kilka statystyk opisowych dla zestawu danych: RangeMeanVarianceNumber of observations0−77218.22791696Range0−772Mean18.2Variance2791Number of observations696 \begin{array}{|c|c|} \hline \text{Range} & 0 - 772 \\ \hline \text{Mean} & 18.2 \\ \hline \text{Variance} & …
Do analizy liczby ptaków z zerowym napełnieniem chciałbym zastosować modele zliczania z zerowym napełnieniem przy użyciu pakietu R pscl . Jednak patrząc na przykład podany w dokumentacji dla jednej z głównych funkcji ( ? Zeroinfl ), zaczynam wątpić, jaka jest prawdziwa zaleta tych modeli. Zgodnie z podanym tam kodem próbki …
Po pierwsze, mam pytanie, czy rozkład Poissona jest „stabilny”, czy nie. Bardzo naiwnie (i nie jestem zbyt pewny co do „stabilnych” rozkładów), opracowałem rozkład liniowej kombinacji rozproszonych RV Poissona, używając iloczynu MGF. Wygląda na to, że dostaję kolejnego Poissona z parametrem równym liniowej kombinacji parametrów poszczególnych RV. Stwierdzam więc, że …
Obecnie pracuję nad projektem z udziałem GLM (i ewentualnie GAM) niektórych danych zliczających w czasie. Normalnie zrobiłbym to w SAS, ale próbuję przejść do R i mam ... problemy. Kiedy dopasowuję GLM do zliczania danych, używając następujących elementów: cdi_model <- glm(counts ~ exposure + covariate + month, data=test, family = …
Czy ktokolwiek może pokazać, w jaki sposób oczekiwana wartość i wariancja zerowego nadciśnionego Poissona, z funkcją masy prawdopodobieństwa f(y)={π+(1−π)e−λ,(1−π)λye−λy!,if y=0if y=1,2....f(y)={π+(1−π)e−λ,if y=0(1−π)λye−λy!,if y=1,2.... f(y) = \begin{cases} \pi+(1-\pi)e^{-\lambda}, & \text{if }y=0 \\ (1-\pi)\frac{\lambda^{y}e^{-\lambda}}{y!}, & \text{if }y=1,2.... \end{cases} gdzie jest prawdopodobieństwem, że obserwacja wynosi zero w procesie dwumianowym, a jest średnią Poissona?ππ\piλλ\lambda …
Czy można zastosować rozkład Poissona do analizy danych ciągłych, a także danych dyskretnych? Mam kilka zestawów danych, w których zmienne odpowiedzi są ciągłe, ale bardziej przypominają rozkład Poissona niż rozkład normalny. Jednak rozkład Poissona jest rozkładem dyskretnym i zwykle dotyczy liczb lub zliczeń.
mgcvOpakowanie Rposiada dwie funkcje montowania interakcji produktów napinacz: te()i ti(). Rozumiem podstawowy podział pracy między nimi (dopasowanie interakcji nieliniowej vs. rozkładanie tej interakcji na główne efekty i interakcję). To, czego nie rozumiem, to dlaczego te(x1, x2)i ti(x1) + ti(x2) + ti(x1, x2)może powodować (nieznacznie) różne wyniki. MWE (dostosowany z ?ti): …
Powiedz, że mam następujący model: Poisson(λ)∼{λ1λ2if t<τif t≥τPoisson(λ)∼{λ1if t<τλ2if t≥τ\text{Poisson}(\lambda) \sim \begin{cases} \lambda_1 & \text{if } t \lt \tau \\ \lambda_2 & \text{if } t \geq \tau \end{cases} I wywnioskowałem, że tylne dla λ1λ1\lambda_1 i pokazano poniżej z moich danych. Czy istnieje bayesowski sposób stwierdzenia (lub określenia ilościowego), czy i …
Korzystając z wikipedii znalazłem sposób na obliczenie funkcji masy prawdopodobieństwa wynikającej z sumy dwóch zmiennych losowych Poissona. Myślę jednak, że moje podejście jest błędne. Niech będą dwiema niezależnymi losowymi zmiennymi Poissona ze średnimi i , gdzie i są stałymi, to funkcję generującą prawdopodobieństwo podaje Teraz, korzystając z faktu, że funkcją …
Jeśli mamy dwie niezależne zmienne losowe i X 2 ∼ P o i s ( λ ) , jaka jest funkcja masy prawdopodobieństwa X 1 + X 2 ?X1∼ B i n o m ( n , p )X1∼Binom(n,p)X_1 \sim \mathrm{Binom}(n,p)X2)∼ P o i s ( λ )X2∼Pois(λ)X_2 \sim \mathrm{Pois}(\lambda)X1+ …
Korzystam z modelu regresji Poissona do zliczania danych i zastanawiam się, czy istnieją powody, aby nie używać solidnego standardowego błędu do szacowania parametrów? Jestem szczególnie zaniepokojony, ponieważ niektóre z moich szacunków bez solidnej nie są znaczące (np. P = 0,13), ale z solidną są znaczące (p <0,01). W SAS jest …
Mam GLMM w postaci: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Kiedy używam drop1(model, test="Chi"), otrzymuję inne wyniki niż w przypadku korzystania Anova(model, type="III")z pakietu samochodowego lub summary(model). Te dwa ostatnie dają te same odpowiedzi. Korzystając z wielu sfabrykowanych danych, odkryłem, że te …
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.