方差比率计算

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

##加载数据处理需要的程序包pacman::p_load(readxl,foreign,plyr,plm,stringr,dplyr,tidyr,knitr,  officer,flextable,stargazer,ggthemes,ggplot2,shiny,DT,reticulate)##读入收益率数据load("F:/科研/开题准备/实证/analsisold/dayreturns.RData")##Lo and MacKinlay 方法计算函数vR1<-function(x,k){  T<-length(x)  m=k*(T-k+1)*(1-k/T)  mu=mean(x)  len<-k:T  if(k==1){    sim=cumsum(x)  }else{    sim=tail(cumsum(x),-(k-1))  }  sigmak<-sum((sim-len*mu)^2)/m  sigma1<-var(x)  VRk<-sigmak/(sigma1)  return(VRk)}##Wright方法计算函数vR2<-function(x,k){  T<-length(x)  m=k*T  mu=mean(x)  len<-k:T  if(k==1){    sim=cumsum(x)  }else{    sim=tail(cumsum(x),-(k-1))  }  sigmak<-sum((sim-len*mu)^2)/m  sigma1<-var(x)  VRk<-sigmak/(sigma1)  return(VRk)}##计算每月统计量及噪声交易指标VR<-stockdata%>%mutate(yms= format(as.Date(trade_date),"%Y-%m"))%>%  group_by(ts_code,yms)%>%summarise(VR51=vR1(dayreturn,5),VR52=vR2(dayreturn,5))%%  ungroup()%>%mutate(NoiseTrad1=abs(VR51-1),NoiseTrad2=abs(VR52-1))

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：方差比 Mackinlay Analsis Foreign Officer

要多写一些介绍就好了。

非常感谢您的分享

cheetahfly 发表于 2020-9-27 19:15
要多写一些介绍就好了。

详细介绍请参考https://yjx.netlify.app/research ... %E8%AE%A1%E7%AE%97/

您好，方便贴出来您的RData吗？是用R软件跑的数据吗？

lilachan 发表于 2021-10-26 10:23
您好，方便贴出来您的RData吗？是用R软件跑的数据吗？

是的，上面那个连接有数据