Przed przesłaniem mojej metaanalizy chcę utworzyć wykres lejka w celu przetestowania niejednorodności i stronniczości publikacji. Mam sumaryczny rozmiar efektu i rozmiary efektu z każdego badania, które przyjmują wartości od -1 do +1. Mam wielkości próbek n1, n2 dla pacjentów i kontroli z każdego...
Jak powinienem złożyć składnię rmafunkcji z pakietu metafor , aby uzyskać wyniki w następującym prawdziwym przykładzie małej metaanalizy? (statystyczny efekt SMD z efektem losowym) study, mean1, sd1, n1, mean2, sd2, n2 Foo2000, 0.78, 0.05, 20, 0.82, 0.07, 25 Sun2003, 0.74, 0.08, 30, 0.72, 0.05,...
Chciałbym dowiedzieć się o wartościach (x, y)używanych podczas kreślenia plot(b, seWithMean=TRUE)w pakiecie mgcv . Czy ktoś wie, jak mogę wyodrębnić lub obliczyć te wartości? Oto przykład: library(mgcv) set.seed(0) dat <- gamSim(1, n=400, dist="normal", scale=2) b <- gam(y~s(x0),...
Pytanie dla początkujących tutaj. Obecnie wykonuję regresję nieparametryczną przy użyciu pakietu np w R. Mam 7 funkcji i stosując metodę brutalnej siły zidentyfikowałem najlepszy 3. Ale wkrótce będę miał więcej niż 7 funkcji! Moje pytanie brzmi: jakie są obecnie najlepsze metody wyboru cech dla...
Muszę narysować złożoną grafikę do wizualnej analizy danych. Mam 2 zmienne i dużą liczbę przypadków (> 1000). Na przykład (liczba wynosi 100, jeśli dyspersja jest mniej „normalna”): x <- rnorm(100,mean=95,sd=50) y <- rnorm(100,mean=35,sd=20) d <- data.frame(x=x,y=y) 1) Muszę...
Powiedzmy, że regresuję Y na X1 i X2, gdzie X1 jest zmienną numeryczną, a X2 jest czynnikiem o czterech poziomach (A: D). Czy jest jakiś sposób na zapisanie funkcji regresji liniowej lm(Y ~ X1 + as.factor(X2)), abym mógł wybrać konkretny poziom X2 - powiedzmy B - jako linię...
Kończę analizę dużej grupy danych. Chciałbym wziąć model liniowy zastosowany w pierwszej części pracy i ponownie go dopasować za pomocą liniowego modelu mieszanego (LME). LME byłby bardzo podobny, z tym wyjątkiem, że jedna ze zmiennych zastosowanych w modelu byłaby zastosowana jako efekt losowy. Te...
AKTUALIZACJA: daszek używa teraz foreachwewnętrznie, więc to pytanie nie jest już tak naprawdę istotne. Jeśli możesz zarejestrować działający backend równoległy foreach, Caret go użyje. Mam pakiet karetki dla R i jestem ciekawy w użyciu trainfunkcji do krzyżowej weryfikacji moich modeli. Chcę...
Chciałbym uzyskać graficzną reprezentację korelacji w artykułach, które zebrałem do tej pory, aby łatwo zbadać relacje między zmiennymi. Kiedyś rysowałem (niechlujny) wykres, ale teraz mam za dużo danych. Zasadniczo mam stół z: [0]: nazwa zmiennej 1 [1]: nazwa zmiennej 2 [2]: wartość...
Mam pewne dane w R, zapisane na liście. Myśleć d <- c(1,2,3,4) chociaż to nie są moje dane. Jeśli następnie wprowadzę polecenie plot(density(d, kernel="gaussian", width=1)) następnie otrzymuję oszacowanie prawdopodobieństwa gęstości jądra, gdzie jądro jest normalnie normalne. Jeśli...
tło Mam dane z badań terenowych, w których istnieją cztery poziomy leczenia i sześć powtórzeń w każdym z dwóch bloków. (4x6x2 = 48 obserwacji) Bloki są oddalone od siebie o około 1 milę, aw obrębie bloków znajduje się siatka o powierzchni 42, 2m x 4m i chodnik o szerokości 1m; moje badanie...
Mam 5 serii całkowitych zwrotów z rynków wschodzących, dla których prognozuję przyszłe zwroty z jednego okresu (1 rok). Chciałbym zbudować portfel zoptymalizowany pod kątem średnich wariancji Markowitza serii 5, wykorzystując historyczne wariancje i kowariancje (1) oraz własne prognozy oczekiwanych...
Muszę przeprowadzić analizę czynnikową dla zestawu danych złożonego ze zmiennych dychotomicznych (0 = tak, 1 = nie) i nie wiem, czy jestem na dobrej drodze. Za pomocą tetrachoric()tworzę macierz korelacji, na której działam fa(data,factors=1). Wynik jest bardzo zbliżony do wyników, które otrzymuję...
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 R działającą na Amazon EC2, używając zmodyfikowanej wersji bioprzewodnika AMI . Obecnie...
Mam pytanie dotyczące tego, czy należy użyć przesunięcia. Załóż bardzo prosty model, w którym chcesz opisać (ogólną) liczbę bramek w hokeju. Masz więc bramki, liczbę rozegranych gier i zmienny manekin „napastnik”, który jest równy 1, jeśli gracz jest napastnikiem, a 0 w przeciwnym razie. Który z...
W poprzednim pytaniu zapytałem o dopasowanie rozkładów do niektórych niegaussowskich danych empirycznych. Zasugerowano mi offline, że mogę spróbować założyć, że dane są gaussowskie i najpierw dopasować filtr Kalmana. Następnie, w zależności od błędów, zdecyduj, czy warto opracować coś bardziej...
Wiem, że ad.test () może być używany do testowania normalności. Czy można uzyskać ad.test w celu porównania dystrybucji z dwóch próbek danych? x <- rnorm(1000) y <- rgev(2000) ad.test(x,y) Jak mogę wykonać test Andersona-Darlinga na 2
Zawsze uważałem, że czas nie powinien być wykorzystywany jako predyktor w regresjach (w tym gam), ponieważ wtedy po prostu „opisałby” sam trend. Jeśli celem badań jest znalezienie parametrów środowiskowych, takich jak temperatura itp., Które wyjaśniają wariancję, powiedzmy, aktywności zwierzęcia,...
Używam pakietu randomForest w R do opracowania losowego modelu lasu w celu wyjaśnienia ciągłego wyniku w „szerokim” zestawie danych z większą liczbą predyktorów niż próbek. W szczególności dopasowuję jeden model RF, umożliwiając procedurze wybranie zestawu ~ 75 zmiennych predykcyjnych, które moim...
Chciałbym połączyć prognozę i prognozę wsteczną (mianowicie prognozowane wartości przeszłe) zestawu danych szeregów czasowych w jeden szereg czasowy, minimalizując średni błąd przewidywania kwadratu. Powiedzmy, że mam szeregi czasowe z lat 2001–2010 z luką dla roku 2007. Byłem w stanie prognozować...