প্রায় কোনও দ্বিখণ্ডিত কপুলা কিছু ননজারো পারস্পরিক সম্পর্কের সাথে একসাথে সাধারণ র্যান্ডম তারতম্যের উত্পাদন করে (কিছু শূন্য দেবে তবে তারা বিশেষ ক্ষেত্রে)। তাদের মধ্যে বেশিরভাগ (প্রায় সবগুলি) একটি অ-সাধারণ যোগফল উত্পাদন করবে।
কিছু কপুলার পরিবারে কোনও পছন্দসই (জনসংখ্যা) স্পিয়ারম্যান পারস্পরিক সম্পর্ক তৈরি করা যেতে পারে; অসুবিধাটি কেবল সাধারণ মার্জিনের জন্য পিয়ারসন পারস্পরিক সম্পর্ক সন্ধানে; এটি নীতিগতভাবে কার্যকর, তবে বীজগণিত সাধারণভাবে বেশ জটিল হতে পারে। [তবে, আপনার যদি জনসংখ্যার স্পিয়ারম্যান পারস্পরিক সম্পর্ক থাকে তবে পিয়ারসন পারস্পরিক সম্পর্ক - কমপক্ষে গাউসের মতো হালকা লেজযুক্ত মার্জিনের জন্য - এটি অনেক ক্ষেত্রে খুব বেশি দূরে নাও থাকতে পারে]]
কার্ডিনালের প্লটের প্রথম দুটি উদাহরণ বাদে সমস্তকে অ-স্বাভাবিক অঙ্ক দেওয়া উচিত।
কিছু উদাহরণ - প্রথম দুটি হ'ল একই কপুলার পরিবার থেকে কার্ডিনালের উদাহরণ বিভারিয়েট বিতরণের পঞ্চম হিসাবে, তৃতীয়টি হ'ল অবক্ষয়।
উদাহরণ 1:
ক্লেটন কোপুলা ( )θ=−0.7
এখানে যোগফলটি খুব স্পষ্টভাবে পিক এবং মোটামুটি দৃ strongly়ভাবে ডান স্কিউ w
উদাহরণ 2:
ক্লেটন কোপুলা ( )θ=2
এখানে যোগফলটি হালকাভাবে বাম স্কু হয়। এটি সবার কাছে একেবারেই সুস্পষ্ট না হলে, আমি এখানে বিতরণটি উল্টিয়েছি (অর্থাত আমাদের ফ্যাকাশে বেগুনিতে একটি হিস্টোগ্রাম রয়েছে ) এবং এটি সুপারিপোজড করেছি যাতে আমরা অসম্পূর্ণতা আরও পরিষ্কারভাবে দেখতে পারি:−(x+y)
আমরা সহজেই যোগফলের স্কিউনেসের দিকটি বিনিময় করতে পারি যাতে নেতিবাচক সম্পর্কটি বাম স্কিউ এবং ডান স্কিউয়ের সাথে ধনাত্মক পারস্পরিক সম্পর্কের সাথে চলে যায় (উদাহরণস্বরূপ, এবং প্রতিটির মধ্যে নিয়ে উপরের কেসগুলি - নতুন ভেরিয়েবলের পারস্পরিক সম্পর্ক পূর্বের মতই হবে তবে যোগফলের বিভাজনটি 0-এর আশপাশে উল্টে থাকবে, স্কিউনেসকে বিপরীত করে তুলবে)।X∗=−XY∗=−Y
অন্যদিকে, আমরা যদি কেবল তাদের একটিটিকে অবহেলা করি তবে আমরা পারস্পরিক সম্পর্কের চিহ্ন (তবে এটির দিক নয়) দিয়ে স্কিউনেসের শক্তির মধ্যে সংযোগ পরিবর্তন করব change
বাইভারিয়েট বিতরণ এবং সাধারণ মার্জিনের সাথে কী ঘটতে পারে তার অনুভূতি পেতে কয়েকটি পৃথক কপুলা নিয়ে ঘুরে বেড়ানোর পক্ষে এটি মূল্যবান।
টি-কোপুলার সাথে গাউসিয়ান মার্জিনগুলি কপুলাসের বিশদ সম্পর্কে খুব বেশি চিন্তা না করে পরীক্ষা-নিরীক্ষা করা যেতে পারে (পারস্পরিকভাবে বিভাজনীয় টি থেকে উত্পন্ন, যা সহজ, তারপরে সম্ভাব্য ইন্টিগ্রাল ট্রান্সফর্মের মাধ্যমে ইউনিফর্ম মার্জিনে রূপান্তরিত হবে, তারপরে ইউনিফর্ম মার্জিনগুলি গাউসিয়ানের মাধ্যমে রূপান্তর করবে) বিপরীত সাধারণ সিডিএফ)। এটিতে একটি নন-নরমাল-তবে-সিমমেট্রিক যোগফল থাকবে। সুতরাং আপনার কাছে কপুলা-প্যাকেজগুলি সুন্দর না থাকলেও আপনি কিছু কিছু নিখরচায় করতে পারেন (উদাহরণস্বরূপ, যদি আমি এক্সেলের মধ্যে চূড়ান্তভাবে একটি উদাহরণ দেখানোর চেষ্টা করছিলাম, আমি সম্ভবত টি-কোপুলা দিয়ে শুরু করব)।
-
উদাহরণ 3 : (এটি আমার শুরুতে শুরু করা উচিত এটির মতো)
স্ট্যান্ডার্ড ইউনিফর্ম উপর ভিত্তি করে একটি কপুলা বিবেচনা করুন এবং এবং জন্য জন্য প্রদান করুন । ফলাফলটি এবং জন্য সমান প্রান্তিক রয়েছে , তবে দ্বিপরিচালিত বিতরণ অধঃপতন হয়। উভয় মার্জিনকে সাধারণ to এ রূপান্তর করা , আমরা জন্য এমন একটি বিতরণ পেয়েছি যা দেখতে এরকম দেখাচ্ছে:UV=U0≤U<12V=32−U12≤U≤1UVX=Φ−1(U),Y=Φ−1(V)X+Y
এই ক্ষেত্রে তাদের মধ্যে পারস্পরিক সম্পর্ক 0.66 এর কাছাকাছি।
সুতরাং আবারও এবং সাথে একটি স্বাভাবিক সম্পর্ক রয়েছে (এই ক্ষেত্রে, স্বতন্ত্রভাবে) অ-সাধারণ যোগফল - কারণ তারা বিভাজনযুক্ত সাধারণ নয়।XY
[ONE কেন্দ্রে আলোকসম্পাতের দ্বারা সম্পর্কযুক্তরূপে একটি সীমার উৎপন্ন পারে (ইন , জন্য মধ্যে ), প্রাপ্ত করতে । এগুলির স্পাইকটি 0 এ থাকবে এবং তারপরে একটি ফাঁক দিয়ে স্বাভাবিক লেজ থাকবে]Uসি[0,1(12−c,12+c)cভি[0,12]V
কিছু কোড:
library("copula")
par(mfrow=c(2,2))
# Example 1
U <- rCopula(100000, claytonCopula(-.7))
x <- qnorm(U[,1])
y <- qnorm(U[,2])
cor(x,y)
hist(x,n=100)
hist(y,n=100)
xysum <- rowSums(qnorm(U))
hist(xysum,n=100,main="Histogram of x+y")
plot(x,y,cex=.6,
col=rgb(0,100,0,70,maxColorValue=255),
main="Bivariate distribution")
text(-3,-1.2,"cor = -0.68")
text(-2.5,-2.8,expression(paste("Clayton: ",theta," = -0.7")))
দ্বিতীয় উদাহরণ:
#--
# Example 2:
U <- rCopula(100000, claytonCopula(2))
x <- qnorm(U[,1])
y <- qnorm(U[,2])
cor(x,y)
hist(x,n=100)
hist(y,n=100)
xysum <- rowSums(qnorm(U))
hist(xysum,n=100,main="Histogram of x+y")
plot(x,y,cex=.6,
col=rgb(0,100,0,70,maxColorValue=255),
main="Bivariate distribution")
text(3,-2.5,"cor = 0.68")
text(2.5,-3.6,expression(paste("Clayton: ",theta," = 2")))
#
par(mfrow=c(1,1))
তৃতীয় উদাহরণের জন্য কোড:
#--
# Example 3:
u <- runif(10000)
v <- ifelse(u<.5,u,1.5-u)
x <- qnorm(u)
y <- qnorm(v)
hist(x+y,n=100)