इस फ़ाइल को क्रिएटिव कॉमन्स श्रेय-समानसांझा 2.5 साधारण लाइसेंस के अंतर्गत लाइसेंस किया गया है।

आप खुलकर:

• बाँट सकते हैं – रचना की प्रतिलिपि बना सकते हैं, बाँँट सकते हैं और संचारित कर सकते हैं
• रीमिक्स कर सकते हैं – कार्य को अनुकूलित कर सकते हैं

निम्नलिखित शर्तों के अंतर्गत:

• श्रेय – यह अनिवार्य है कि आप यथोचित श्रेय प्रदान करें, लाइसेंस की कड़ी प्रदान करें, और अगर कोई बदलाव हुए हों तो उन्हें इंगित करें। आप ऐसा किसी भी उचित तरीके से कर सकते हैं, लेकिन किसी भी तरह उससे यह नहीं संकेत नहीं किया जाना चाहिए कि लाइसेंसधारी द्वारा आपको अथवा आपके इस प्रयोग का समर्थन किया जा रहा हो।
• समानसांझा – अगर आप इस रचना में कोई बदलाव करते हैं या इसपर आधारित कुछ रचित करते हैं तो आप अपने योगदान को सिर्फ इसी या इसके सामान किसी लाइसेंस के अंतर्गत बाँट सकते हैं।

https://creativecommons.org/licenses/by-sa/2.5CC BY-SA 2.5 Creative Commons Attribution-Share Alike 2.5 truetrue

R source:

library(tseries)
library(RSvgDevice)

quote<-function(inst, nDs) {
  if(!inherits(try(open(url("http://quote.yahoo.com")), silent = TRUE), "try-error")) {
    start <- strftime(as.POSIXlt(Sys.time() - nDs*24*3600), format="%Y-%m-%d") 
    end <- strftime(as.POSIXlt(Sys.time()), format = "%Y-%m-%d") 
    x <- get.hist.quote(instrument=inst, start=start, end=end, quote = c("AdjClose"), compression="m",  quiet=T, retclass="zoo")
  }
  return(x)
}

yield=2.9
nDays=1042
comp<-c("^DJI", "AA", "AIG", "AXP", "BA", "C", "CAT", "DD", "DIS", "GE", "GM", "HD", "HON", "HPQ", "IBM", "INTC", "JNJ", "JPM", "KO", "MCD", "MMM", "MO", "MRK", "MSFT", "PFE", "PG", "T", "UTX", "VZ", "WMT", "XOM")

comp.quote<-quote(comp[1],nDays)
from= as.Date(index(comp.quote[1,1]),format="%Y-%m-%d")
to=   as.Date(index(comp.quote[length(comp.quote[,1]),1]),format="%Y-%m-%d")

for (c in 2:length(comp)) {
  comp.quote<-merge(comp.quote,quote(comp[c],nDays))
}

rdt<-12*diff(comp.quote)/comp.quote*100

points=array(NA,dim=c(length(comp),2))
for (p in 1:length(comp)) {
  points[p,1]= cov(rdt[,1],rdt[,p])/var(rdt[,1])
  points[p,2]=mean(rdt[,p])
}

devSVG(file="CAPM-SML.svg", width=10, height=8, bg="transparent", fg="black", onefile=TRUE, xmlHeader=TRUE)

par(lwd=2)
plot(points[2:length(comp),],type="p",col="blue",xlab="beta",ylab="return (in % p.a.)",main="CAPM: Security Market Line",xlim=range(0,points[,1]))
lines(points[1,1],points[1,2],col="red",type="p")
abline(a=yield,b=points[1,2]-yield,col="purple")
text(points[2:length(comp),1],points[2:length(comp),2],comp,pos=3,col="blue")
text(points[1,1],             points[1,2],           comp[1],pos=3,col="red")

lm1<-lm(points[,2]~points[,1])
abline(lm1,col="darkgrey",lwd=1)
 
text(x=par("usr")[2], y=par("usr")[3]+0.5, labels=paste("Dow Jones monthly data from", format.Date(from,"%d %b %Y"), "to", format.Date(to,"%d %b %Y"),""), col="grey",adj=c(1,0) )
grid(lwd=1)

dev.off()

 

इस W3C-अनिर्दिष्ट वेक्टर चित्र को R की मदद से बनाया गया था.