мне кажется, это что-то странное в circular

> tmp <- data.table(x = rep(1:5, each = 2), y = c('A','B'))
> tmp[, list(A = mean(x), A1 = (mean(x)), B = mean(x), B1 = (mean(x))), by = y]
   y A A1 B B1
1: A 3  3 3  3
2: B 3  3 3  3

Дорогие, а никто не знает, можно ли в rmarkdown оргаинизовать partial upload страниц таблицы DT::datatable()? Я хочу чтобы вот такие сайты не грузились столетиями.

> tmp <- data.table(x = rep(1:5, each = 2), y = c('A','B'))
> tmp[, list(A = mean(x), A1 = (mean(x)), B = mean(x), B1 = (mean(x))), by = y]
   y A A1 B B1
1: A 3  3 3  3
2: B 3  3 3  3

сорри, rvm из CircStats
проблема в том, что среднее в скобках отличается от среднего
я тоже думаю, что это прикол Circular, но я не понимаю, как взятие результата вызова функции в скобки может оказывать влияние на ее выполнение

по идее, никак
вон, у меня в примере все одинаково

вооот, я и говорю, колдунство

rvm кстати не важен, вот так тоже разные значения
data.table(x = rnorm(500, pi, 0.5), y = c('A','B'))[,.(A = (mean(circular(x))), B = mean(circular(x))), by = y]

судя по всему, вызов mean без скобок не приводит к вызову mean.circular

dt[,.(A = (mean(circular(x))), B = mean(circular(x)), C = circular::mean.circular(circular(x))), by = y]
   y         A        B         C
1: A 0.6033983 1.949538 0.6033983
2: B 0.3794051 2.191128 0.3794051

но почему, убей не понимаю

dt[,.(A = (mean(circular(x))), B = mean(circular(x)), C = circular::mean.circular(circular(x))), by = y]
   y         A        B         C
1: A 0.6033983 1.949538 0.6033983
2: B 0.3794051 2.191128 0.3794051

Гуглил меньше минуты.
https://github.com/Rdatatable/data.table/issues/3533

ааа, спасибо огромное
я просто не понимал, что гуглить

Коллеги, добрый день. Я тут почти доделал пакет для AmoCRM, кто-нибудь хочет помочь потестить? Было бы круто, так как у меня только один аккаунт и чем раньше отловлю баги, тем раньше выкачу рабочую версию. В общем, желающие — напишите в личку, как допилю доку до MVP, скину вам

Мой опыт в академии, показывает, что лучше как можно скорее выкатить, что есть, а потом ждать issues на гитхабе, писем на почте и прочих пожеланий. Сразу идеально не получиться. Да и до идеала можно доводить столетиями

Ну, если никто не отзовется, то так и поступлю, да. Подумал, что проще грубые баги поправить при наличии помощи, а изи баги уже версиями фиксить

а я даже не понимаю, про что речь %(

В чате аналитиков не так давно спрашивали

Я и спрашивал)

Не, наоборот спрашивали есть ли готовый пакет

а, не видел

Привет всем, подскажите, пожалуйста, что считает эа функция?

Phi <- function(a, b) 
{
  r<-cov(a,b)/(var(a)*var(b))
  Sr<-(1-r^2)/sqrt(length(unique(a))-1)
  return(abs(r/Sr))
}

Думаю, что гуглить, но не могу же я формулу загуглить

Понимаю только, что r - корреляция

Phi <- function(a, b) 
{
  r<-cov(a,b)/(var(a)*var(b))
  Sr<-(1-r^2)/sqrt(length(unique(a))-1)
  return(abs(r/Sr))
}

Чем-то похоже на отношение коэффициента корреляции к его стандартной ошибке