【发布时间】:2020-10-01 20:51:57
【问题描述】:
解析错误(文本 = 脚本):文本参数解析错误:符号 必须使用反引号(第 26 行,字符 47 到 67)编写
当我尝试使用以下代码时出现上述错误。奇怪的是我的代码中没有任何 ``.. 对我所缺少的有什么建议吗?第 26 行以 if(sum(gdata..
data = data.frame( Well, Prod, ProdDate )
data$ProdDate = as.POSIXct(data$ProdDate, origin = "1970-01-01", tz="UTC")
data$ProdDate <- interpNA(data$ProdDate, method = c("linear")
minDate = aggregate(x=data$ProdDate, by=list(data$Well), FUN=min)
colnames(minDate) <- c("Well","MinDate")
data = merge(data, minDate, by="Well", all.x=TRUE)
# Calculate Days On
data$MinDate = as.POSIXct(data$MinDate, origin="1970-01-01", tz="UTC")
data$DaysOn = as.numeric(data$ProdDate - data$MinDate, units="days")
data$Predict = rep(NA, nrow(data))
data$Qi = rep(NA, nrow(data))
data$A = rep(NA, nrow(data))
data$B = rep(NA, nrow(data))
gdata <- data;
min.RSS <- function(data, par) {
with(data, sum(((par[1]/((1+par[2]*par[3] * DaysOn)^(1/par[3]))) - Oil)^2))
}
wells = unique( na.omit(gdata$Well) )
for(well in wells) {
idx = gdata$Well == well
if(sum(gdata$Prod[idx], na.rm=TRUE) > 0) {
startQi = max(data$Prod[idx], na.rm=TRUE)
fit <- optim(par = c(startQi, .2), min.RSS, x = gdata[idx,])
gdata$Predict[idx] = fit$par[1]/((1+fit$par[2]*1*gdata$DaysOn[idx])^(1/1))
gdata$Qi[idx] = fit$par[1]
gdata$A[idx] = fit$par[2]
gdata$B[idx] = 1
}
}
gdata$ProdDateCheck <- gdata$ProdDate
Predict <- gdata #gdata[, c("MinDate", "DaysOn", "Predict", "Qi", "A", "B")]
【问题讨论】: