Как вариант
require(MultinomialCI)
m = round(multinomialCI(c(50,100,120,88,44), 0.05),3) #собственно, расчеты, вместо вектора вставить ваши значения
#вывод значений в более-менее понятной форме
print(paste("0 баллов: [", m[1,1], m[1,2], "]"));
print(paste("1 балл: [", m[2,1], m[2,2], "]"));
print(paste("2 балла: [", m[3,1], m[3,2], "]"));
print(paste("3 балла: [", m[3,1], m[3,2], "]"));
print(paste("4 балла: [", m[3,1], m[3,2], "]"));