有了这个R代码,我现在有了每个Y的比例和形状参数:
library(QRM)
y <- read.csv(file="C:\\Users\\simon\\Desktop\\Durchschnit\\ESGHigh10%.csv", header= TRUE, sep=";", fileEncoding="UTF-8-BOM")
y1 <- unlist(y)
months <- 1:240
n <- 2557
time <- (rep(months, each=n))
z <- read.csv(file="C:\\Users\\simon\\Desktop\\Covariats\\Covariats-USA-Test1.csv", header= TRUE, sep=";", fileEncoding="UTF-8-BOM")
c1= (rep(z[,1], each=n))
c2= (rep(z[,2], each=n))
final <- data.frame(time, y1, c1, c2)
x <- final[complete.cases(final), ]
mod0 <- gamGPDfit(x, threshold=5.88, datvar="y1", xiFrhs=~c1+c2, nuFrhs=~c1+c2)
现在,我想知道哪些协变量会影响我的Y,您是否有任何想法或R代码可以帮助我。 预先谢谢你!