VK
Size: a a a
2021 March 15
VK
Всем спасибо, вери хелпфул блин.
G
pvals <- replicate(1000, {x <- rnorm(100, mean = 0, sd=1);
mn <- mean(x);
s <- sd(x);
x <- (x - mn) / s;
ks.test(x, y='pnorm')$p.value})По ходу надо просто сравнивать не со станд нормальным а с нормальным с оцененными сигма и мю, и все в шоколаде
DP
Иван Поздняков
set.seed(42)
many_p <- replicate(100000, ks.test(rnorm(100), "pnorm")$p.value)
hist(many_p)
И действительно, спасибо, не знал
many_p <- replicate(100000, ks.test(rnorm(100), "pnorm")$p.value)
hist(many_p)
И действительно, спасибо, не знал
у вас как раз странно, ибо вы не передаете в pnorm доп параметров. То есть тестируете простую гипотезу, которая верна
ИП
стоп, вот это должно решать проблему (согласно , а результат получается тот же:
set.seed(42) many_p_withparams <- replicate(100000, ks.test(rnorm(100), "pnorm", mean = 0, sd = 1)$p.value) hist(many_p_withparams)
set.seed(42) many_p_withparams <- replicate(100000, ks.test(rnorm(100), "pnorm", mean = 0, sd = 1)$p.value) hist(many_p_withparams)
DP
у вас как раз странно, ибо вы не передаете в pnorm доп параметров. То есть тестируете простую гипотезу, которая верна
более-менее равномерное
G
И сравнивать с другим распределением :( тогда ок. Датасайнтисты захватили чат
IF
Кажись оно, да! Спасибо!
Хорошо. Тогда рекомендую гуглить, я сам был без понятия, но было любопытно, а гугл все знает :)
DP
German
По ходу надо просто сравнивать не со станд нормальным а с нормальным с оцененными сигма и мю, и все в шоколаде
pvals <- replicate(1000, {x <- rnorm(100, mean = 0, sd=1); mn <- mean(x); s <- sd(x); ks.test(x, y='pnorm', mn, s)$p.value})ИП
у вас как раз странно, ибо вы не передаете в pnorm доп параметров. То есть тестируете простую гипотезу, которая верна
честно говоря, я не особо тыкал
ks.test(), но вроде как если он оценивает параметры, то он их оценивает по выборке, а если передавать через элипсис, то они прокидываются в функцию распределения как заданные параметрыVK
Хорошо. Тогда рекомендую гуглить, я сам был без понятия, но было любопытно, а гугл все знает :)
Да, просто не было в выдаче ничего релевантного почему-то. Теперь изи!
DP
pvals <- replicate(1000, {x <- rnorm(100, mean = 0, sd=1); mn <- mean(x); s <- sd(x); ks.test(x, y='pnorm', mn, s)$p.value})
G
я уже прочитал, другой распределений нужно
IF
Да, просто не было в выдаче ничего релевантного почему-то. Теперь изи!
"tg solutions chemistry" был запрос, который дал ответ
VK
"tg solutions chemistry" был запрос, который дал ответ
Угу, я drug substance искала
IF
Умение гуглить - самое полезный навык в жизни
G
Умение гуглить - самое полезный навык в жизни
а вовсе не математика
DP
нагуглить можно любую хрень