2017-01-19 4 views
0

Я пытаюсь сохранить n значение из pwr.t.test в моем FWD_SummaryStats df, но каждый раз, когда я запускаю код, я сбой R. Точная ошибка код «Это приложение попросило Runtime прекратить его необычным способом ...«Сохранение значения pwr.t.test в итоговой статистике dplyr

Я не уверен, почему это происходит, так как работа каждой статистики дает правильный результат. Любая помощь будет очень признательна!

FWD_SummaryStats<-summarize(FWD1, count = n(), min = min(FWD_1_Fuel_Load, na.rm = T), 
        mean = mean(FWD_1_Fuel_Load, na.rm = T), median = median(FWD_1_Fuel_Load, na.rm = T), 
        max = max(FWD_1_Fuel_Load, na.rm = T), varience = var(FWD_1_Fuel_Load, na.rm = T), 
        SD = sd(FWD_1_Fuel_Load, na.rm = T), 
        n = pwr.t.test(d=((1.15*mean)-mean)/SD, sig.level=0.05,power = 0.95,type= "one.sample",alternative = "two.sided")$n 
) 

Спасибо, Кев

+0

Вы можете указать данные образца? –

+0

Это не похоже на отчет о сбое R. Используете ли вы что-то иное, чем CRAN, поставляемый вариант R? –

+0

@ 42 Нет, это официальная версия. Весьма странно, что функция суммирования может упасть, когда все эти показатели работают нормально при вычислении индивидуально. – AlphaKevy

ответ

0

Для кого-то, кто имеет эту проблему, проблема заключается в этой строке:

n = pwr.t.test(d=((1.15*mean)-mean)/SD, sig.level=0.05, power = 0.95,type="one.sample", alternative = "two.sided")$'n' 

$ N должно быть в окружении цитат, чтобы сформировать $ н ', чтобы сохранить значение.

SummaryStatistics<- 
      summarize(FWD1, count = n(), min = min(FWD_1_Fuel_Load, na.rm = T), 
        mean = mean(FWD_1_Fuel_Load, na.rm = T), median = median(FWD_1_Fuel_Load, na.rm = T), 
        max = max(FWD_1_Fuel_Load, na.rm = T), varience = var(FWD_1_Fuel_Load, na.rm = T), 
        SD = sd(FWD_1_Fuel_Load, na.rm = T), 
        n = pwr.t.test(d=((1.15*mean)-mean)/SD, sig.level=0.05, power = 0.95,type="one.sample", alternative = "two.sided")$'n' 
) 

Я надеюсь, что это поможет тем, кто застрял в одной и той же проблеме.

Kev

+0

Если вы действительно можете сбой CRAN R, вам следует отправить отчет об ошибке, отправив по электронной почте автора пакета и список рассылки r-devel. –