OK
German
Это будет неверное но достаточное определение
Но зачем? Если человеку это знание реально нужно, он и правильное определение осилит, а если не осилит, значит, для его задач это нерелевантно
Size: a a a
2021 March 15
АС
А самое главное — большинству биологов для их задач она и не нужна
Тем, которые работают в готовых программках - да. Тем, кто сами пишут - желательно бы представлять, что есть такое понятие, как сложность алгоритма и иногда код может выполняться несколько месяцев и простой оптимизацией уменьшить время исполнения до нескольких часов
G
У матана есть и в программировании, и в статистике, кроме прочего, есть функция - выработать привычку нее приходя в сознание спрашивать - "а какие допущения/плюсы/минусу у моей модели/структуры данных?"
В той же статистике - простой гуглеж на проверку соответствия распределения нормальному выдаст в числе прочего тест Колмогорова-Смирнова. Можно его посчитать, а можно задаться вопросом, а можно ли. И еще через какое-то время понять, что все не так просто и как минимум стандартная функция в R не подходит (о чем у нее хэлпе и сказано:) )
В той же статистике - простой гуглеж на проверку соответствия распределения нормальному выдаст в числе прочего тест Колмогорова-Смирнова. Можно его посчитать, а можно задаться вопросом, а можно ли. И еще через какое-то время понять, что все не так просто и как минимум стандартная функция в R не подходит (о чем у нее хэлпе и сказано:) )
А почему колмогоров смирнов не подходит?
OK
Тем, которые работают в готовых программках - да. Тем, кто сами пишут - желательно бы представлять, что есть такое понятие, как сложность алгоритма и иногда код может выполняться несколько месяцев и простой оптимизацией уменьшить время исполнения до нескольких часов
Я за разделение труда. В таком случае надо идти к биоинформатикам и не заниматься самодеятельностью
DP
German
А почему колмогоров смирнов не подходит?
если параметры оценивать из выборки, то не подходит. У него распределение p-value начинает жестко жаться в сторону 1. Надо поправлять. В R стандартном это не сделано
G
Но зачем? Если человеку это знание реально нужно, он и правильное определение осилит, а если не осилит, значит, для его задач это нерелевантно
Потому что в программировании дофига собственно программирования, паттерны проектирования, особенности языков, знание работы ОС - зачем загружать человека тем, что он будет использовать в простейшей версии
OK
German
Потому что в программировании дофига собственно программирования, паттерны проектирования, особенности языков, знание работы ОС - зачем загружать человека тем, что он будет использовать в простейшей версии
Вот как раз этим людям матан не пустой звук и определение они осилят. Мы тут все-таки про скриптинг для биологов
G
если параметры оценивать из выборки, то не подходит. У него распределение p-value начинает жестко жаться в сторону 1. Надо поправлять. В R стандартном это не сделано
Он же непараметрический, не?
G
Вот как раз этим людям матан не пустой звук и определение они осилят. Мы тут все-таки про скриптинг для биологов
Ну, в российских вузах их еще и интегрировать по двойному контуру учат - ну а что, осилят же :)
DP
он проверяет простую гипотезу вида выборка x пришла из распределения с заданными параметрами. Если вместо заданных параметров подставить оцененныее по выборке x, то вы подгоняете свою гипотезу под данные и распределение p-value смещается в сторону 1.
G
он проверяет простую гипотезу вида выборка x пришла из распределения с заданными параметрами. Если вместо заданных параметров подставить оцененныее по выборке x, то вы подгоняете свою гипотезу под данные и распределение p-value смещается в сторону 1.
Коррекция Бесселя на станд отклонение не решает проблему? Если я просто стандартизирую
ИП
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)
И действительно, спасибо, не знал
G
Окей, верю, а чтобы не перекашивало что надо сделать
ИП
German
Коррекция Бесселя на станд отклонение не решает проблему? Если я просто стандартизирую
не это про оценку сд по оценке среднего, немного не про то, но идея похожая, мне кажется
DP
German
Коррекция Бесселя на станд отклонение не решает проблему? Если я просто стандартизирую
нет
DP
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})VK
Ребята, а что такое Tg' ?
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
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})