Styjun

Why are grass strips more dangerous than tarmac?

Is /home directory in root partition mapped to /home partition

What is the best option to connect old computer to modern TV

Imperfective Aspect in German "not since" constructions

The term for the person/group a political party aligns themselves with to appear concerned about the general public

Explain Ant-Man's "not it" scene from Avengers: Endgame

Can you please explain this joke: "I'm going bananas is what I tell my bananas before I leave the house"?

What does War Machine's "Canopy! Canopy!" line mean in "Avengers: Endgame"?

Why was it possible to cause an Apple //e to shut down with SHIFT and paddle button 2?

How do you translate “is all” used at the end of a sentence?

Why is Colorado so different politically from nearby states?

Beginner's snake game using PyGame

What is the correct expression of 10/20, 20/30, 30/40 etc?

What is the right way to float a home lab?

Concise way to draw this pyramid

Why were the Night's Watch required to be celibate?

Is it OK to bring delicacies from hometown as tokens of gratitude for an out-of-town interview?

Was the 1959 Tibetan Uprising really an uprising?

Unconventional Opposites

How to provide realism without making readers think grimdark

Can an old DSLR be upgraded to match modern smartphone image quality

How is it possible for Mordenkainen to be alive during the Curse of Strahd adventure?

Why does my electric oven present the option of 40A and 50A breakers?

Please help me identify this plane

plotting circular density using curve.circular

Side-by-side plots with ggplot2Plot two graphs in same plot in RShading a kernel density plot between two points.How to create a density plot in matplotlib?How to overlay density plots in R?R: Using MLE to do a CFA by handHow to compute the inverse of a close to singular matrix in R?Circular density plot using ggplot2Probability over a circular regionConducting MLE for multivariate case (bivariate normal) in R

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty height:90px;width:728px;box-sizing:border-box;

0

I need to plot the circular density defined in function (f2) in the following R code

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)

the problem is that I get error messages:

Error in lines.circular(x, y, ...) : object 'll' not found
In addition: Warning message:
In sqrt(a) : NaNs produced

I do not know how to deal with that error

r

share|improve this question

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I'm not familiar with library(circular). But if we add aa <- aa[,1:2] before return(aa), your function will work until aa <- aa[,1:50] after that it will fail. Try to use browser() inside your function to check if it behaves as you expected.
  
  – A. Suliman
  Mar 24 at 12:42
  
  

  
  
  
  

add a comment | 

0

I need to plot the circular density defined in function (f2) in the following R code

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)

the problem is that I get error messages:

Error in lines.circular(x, y, ...) : object 'll' not found
In addition: Warning message:
In sqrt(a) : NaNs produced

I do not know how to deal with that error

r

share|improve this question

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I'm not familiar with library(circular). But if we add aa <- aa[,1:2] before return(aa), your function will work until aa <- aa[,1:50] after that it will fail. Try to use browser() inside your function to check if it behaves as you expected.
  
  – A. Suliman
  Mar 24 at 12:42
  
  

  
  
  
  

add a comment | 

I need to plot the circular density defined in function (f2) in the following R code

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)

the problem is that I get error messages:

Error in lines.circular(x, y, ...) : object 'll' not found
In addition: Warning message:
In sqrt(a) : NaNs produced

I do not know how to deal with that error

r

share|improve this question

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

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)

share|improve this question

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

share|improve this question

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

edited Mar 24 at 12:42

A. Suliman

6,62341225

edited Mar 24 at 12:42

A. Suliman

6,62341225

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13

asked Mar 24 at 11:42

Amani Abu BakrAmani Abu Bakr

13