臨床統計もおもしろいですよ、その１

臨床統計もおもしろいですよ、その１at HOSP

臨床統計もおもしろいですよ、その１ - 暇つぶし2ch713:cos(pi/3),sin(pi/3),0)){ fn <- function(x,O,A,B,C){ AO=A-O BO=B-O CO=C-O HO=x[1]*AO+x[2]*BO+(1-x[1]-x[2])*CO AB=B-A AC=C-A c(HO%*%AB,HO%*%AC) } fn1 <- function(x) fn(x,O,A,B,C) x=nleqslv::nleqslv(c(1/3,1/3),fn1)$'x' AO=A-O BO=B-O CO=C-O HO=x[1]*AO+x[2]*BO+(1-x[1]-x[2])*CO h=sqrt(sum(HO^2)) a=sqrt(sum((B-C)^2)) b=sqrt(sum((C-A)^2)) c=sqrt(sum((A-B)^2)) s=(a+b+c)/2 base=sqrt(s*(s-a)*(s-b)*(s-c)) V=1/3*base*h return(V) } Tetra() sqrt(2)/12

次ページ

続きを表示

1を表示