firstfun<-function(data1=x1,data2=x2,key=c)
{
result<-merge(data1,data2,by=key)
return(result)
}
report1<-read.table(file="ReportCard1.txt",header = TRUE)
report2<-read.table(file="ReportCard2.txt",header = TRUE)
newdata<-firstfun(report1,report2,"xh")
#应用一:数据基本处理(车险数据)
Des.Fun<-function(x,...){
Av<-mean(x,na.rm = TRUE)
Sd<-sd(x,na.rm = TRUE)
N<-length(x[!is.na(x)])
Sk<-sum((x[!is.na(x)]-Av)^3/Sd^3)/N
Ku<-sum((x[!is.na(x)]-Av)^4/Sd^4)/N-3
result<-list(avg=Av,sd=Sd,skew=Sk,kurt=Ku)
return(result)
}
ClaimData<-read.table(file="车险数据.txt",header=TRUE)#分隔符不加
with(ClaimData,{
Des1<-sapply(ClaimData[,c(1,5)],FUN = Des.Fun)
print(Des1)
Des2<-tapply(nclaims, INDEX=vehiclegroup, FUN=mean,na.rm
{
result<-merge(data1,data2,by=key)
return(result)
}
report1<-read.table(file="ReportCard1.txt",header = TRUE)
report2<-read.table(file="ReportCard2.txt",header = TRUE)
newdata<-firstfun(report1,report2,"xh")
#应用一:数据基本处理(车险数据)
Des.Fun<-function(x,...){
Av<-mean(x,na.rm = TRUE)
Sd<-sd(x,na.rm = TRUE)
N<-length(x[!is.na(x)])
Sk<-sum((x[!is.na(x)]-Av)^3/Sd^3)/N
Ku<-sum((x[!is.na(x)]-Av)^4/Sd^4)/N-3
result<-list(avg=Av,sd=Sd,skew=Sk,kurt=Ku)
return(result)
}
ClaimData<-read.table(file="车险数据.txt",header=TRUE)#分隔符不加
with(ClaimData,{
Des1<-sapply(ClaimData[,c(1,5)],FUN = Des.Fun)
print(Des1)
Des2<-tapply(nclaims, INDEX=vehiclegroup, FUN=mean,na.rm