r - How to determine whether a points lies in an ellipse -

-

i had posted similar question earlier. trying determine whether point lies within ellipse. generate bivariate normal data , create ellipse. heres code use

 library(mass)  set.seed(1234)  x1<-null  x2<-null  k<-1  sigma2 <- matrix(c(.72,.57,.57,.46),2,2)  sigma2  rho <- sigma2[1,2]/sqrt(sigma2[1,1]*sigma2[2,2])   eta<-replicate(300,mvrnorm(k, mu=c(-2.503,-1.632), sigma2))    p1<-exp(eta)/(1+exp(eta))  n<-60  x1<-replicate(300,rbinom(k,n,p1[,1]))  x2<-replicate(300,rbinom(k,n,p1[,2]))   rate1<-x1/60  rate2<-x2/60   library(car)  dataellipse(rate1,rate2,levels=c(0.05, 0.95))

i need find out whether pair (p1[,1],p1[,2]) lies within area of ellipse above.

dataellipse returns ellipses polygons, use point.in.polygon function sp library check whether points inside ellipse:

ell = dataellipse(rate1, rate2, levels=c(0.05, 0.95))  point.in.polygon(rate1, rate2, ell$`0.95`[,1], ell$`0.95`[,2])

when run following code...

library(mass) set.seed(1234) x1<-null x2<-null k<-1 sigma2 <- matrix(c(.72,.57,.57,.46),2,2) sigma2 rho <- sigma2[1,2]/sqrt(sigma2[1,1]*sigma2[2,2]) eta<-replicate(300,mvrnorm(k, mu=c(-2.503,-1.632), sigma2)) p1<-exp(eta)/(1+exp(eta)) n<-60 x1<-replicate(300,rbinom(k,n,p1[,1])) x2<-replicate(300,rbinom(k,n,p1[,2])) rate1<-x1/60 rate2<-x2/60 library(car) ell = dataellipse(rate1, rate2, levels=c(0.05, 0.95)) library(sp) point.in.polygon(rate1, rate2, ell$`0.95`[,1], ell$`0.95`[,2])

... following output

  [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1  [56] 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [111] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [166] 1 1 1 1 1 1 1 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 [221] 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [276] 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Comments

Post a Comment

Popular posts from this blog

css - Which browser returns the correct result for getBoundingClientRect of an SVG element? -

-
this svg contains rect overflows svg element: <svg id='svg' width='10' height='10'> <rect x='-10' y='-10' width='30' height='30'/> </svg> chrome 28 , opera 12 return getboundingclientrect() svg element width , height of 10. firefox 23 reports width , height of 30. correct? jsfiddle the relevant spec cssom , delegates svg spec if svg element not "have associated css layout box". haven't found definition of "having associated css layout box", correct result seem hinge on definition, getbbox returns 30x30 rect in browsers. this firefox bug fixed now, fix in firefox nightlies , should make through firefox 33 released on 14 october 2014. see bug 530985 details.
Read more

gcc - Calling fftR4() in c from assembly -

-
i trying asm fft implementation www.embeddedsignals.com work. running make gives error: undefined reference fftr4(short*, short*, int) at top of asm have .syntax unified .thumb .global fftr4 .type fftr4, %function and there #define's associate c functions arguments registers followed macros fftr4 function uses. after that, fftr4 called .thumb_func fftr4: stmfd sp!, {r4-r11, lr} mov tmp0, #0 // bit reversed counter mov tmp1, #0 .word 0xf3a2fa92 //rbit r3,r2 //rbit r,n lsl r,#3 and c code follows: void fftr4(short *y, short *x, int n); short x[512]; short y[512]; int i; void loop() { (i=0;i<512;i++) x[i]=0; (i=0;i<512;i=i+8) { x[i+0]=16384; x[i+2]=16384; x[i+4]=-16384; x[i+6]=-16384;} fftr4(y, x, 256); serialusb.println(y[128]); serialusb.println(y[129]); serialusb.println(y[384]); serialusb.println(y[385]); // expected results: //y[128]is 8191; y[129] -8190 y[...
Read more

.htaccess - Matching full URL in RewriteCond -

-
i'm trying block post access site (and return 403), except 2 ip's , 1 url ( http://mydomain.com/about-us/contact-us/ ). here's have: rewritecond %{request_method} post rewritecond %{remote_addr} !=1.2.3.4 rewritecond %{remote_addr} !=4.3.2.1 rewritecond %{request_uri} !^http://mydomain.com/about-us/contact-us/$ rewriterule ^ / [f] these rules allowing access ip properly, contact page still not accessible. i'm not sure i'm using request_uri correctly? i've tried variations relative urls rewritecond %{request_uri} !^/about-us/contact-us$ , rewritecond %{request_uri} !^/about-us/contact-us/$ without luck. what doing wrong please? it looks me rule on file name specific. if relax bit, should work. made few other esthetic tweaks. give try: rewritecond %{request_method} post rewritecond %{remote_addr} !^1.2.3.4 rewritecond %{remote_addr} !^4.3.2.1 rewritecond %{request_uri} !about-us/contact-us rewriterule ^ - [f]
Read more