Мне нужно построить круговую плотность, определенную в функции (f2), в следующем коде R
library(circular)
mu=c(1,1.166) # true values of the mean vector of bivariate normal variables
var.1 =1.096# sigma.1.squared for first normal variable
var.2=1#sigma.2.squared of second normal variable
cov.1.2=-0.31 # covariance between the two normal variables
var.cov<-matrix(c(var.1,cov.1.2,cov.1.2,var.2),nrow=2,ncol=2,byrow=T) #variance covariance matrix of the bivariate normal random variables
f2=function(x){
u=cbind(cos(x),sin(x))
a=drop((u)%*%solve(var.cov)%*%t(u))
b=drop((u)%*%solve(var.cov)%*%(mu))
c=drop(-0.5*t(mu)%*%solve(var.cov)%*%mu)
d=drop(b/sqrt(a))
aa=(cos(x))*(1/(2*pi*a*sqrt(det(var.cov))))*exp(c)*(1+((d*pnorm(d))/dnorm(d)))
return(aa)
}
curve.circular(f2,from=0,to=2*pi,n=101)
проблема в том, что я получаю сообщения об ошибках:
Ошибка в lines.circular (x, y, ...): объект 'll' не найден
Дополнительно: предупреждающее сообщение:
В sqrt (a): произведено NaN
Я не знаю, как справиться с этой ошибкой