m0_67028469 2022-10-19 21:22
浏览 12
已结题

R语言错误,请问如何改正

原始代码:
grid_id<-read.table("C:\Users\Admin\Desktop\VIC\WZ_VICmodel\soil\arcgis_Output_soil.txt",header =T,sep=",")[,3]
texture_id<-read.table("C:\Users\Admin\Desktop\VIC\WZ_VICmodel\soil\arcgis_Output_soil.txt",header =T,sep=",")[,34]
a<-factor(texture_id)
fdata<-data.frame(grid_id,texture_id)
fdata<-fdata[order(fdata$grid_id),]
texture_type=rep(0,13)
grid_code_T = c()
r=1
t=0
vn=0
l1 = nrow(fdata)
for(i in 1:l1){
grid=fdata$grid_id[i]
texture=fdata$texture_id[i]
vn=vn+1
if(texture>=1&&texture<=13)
texture_type[texture]=texture_type[texture]+1
if(i==l1||grid!=fdata$grid_id[i+1]){
max = texture_type[1]
for(i in 1:13){
if(texture_type[i]>max){
max=texture_type[i]
t=i
}
}
grid_code_T[r] = t
r=r+1
vn=0
texture_type=rep(0,13)
}
}
grid_code_T
write.table(grid_code_T,"C:\Users\Admin\Desktop\VIC\WZ_VICmodel\soil\output_finally_soil.txt",col.names=F,row.names=F,quote=F)
出现的问题
Error in if (texture >= 1 && texture <= 13) texture_type[texture] = texture_type[texture] + :
missing value where TRUE/FALSE needed
请问一下怎么改变

  • 写回答

0条回答 默认 最新

    报告相同问题?

    问题事件

    • 系统已结题 10月27日
    • 创建了问题 10月19日

    悬赏问题

    • ¥15 三分类机器学习模型可视化分析
    • ¥15 本地测试网站127.0.0.1 已拒绝连接,如何解决?(标签-ubuntu)
    • ¥50 Qt在release捕获异常并跟踪堆栈(有Demo,跑一下环境再回答)
    • ¥30 python,LLM 文本提炼
    • ¥15 关于将inet引入的相关问题
    • ¥15 关于一个倒计时的操作和显示设计
    • ¥15 提问STK的问题,哪位航天领域的同学会啊
    • ¥15 苹果系统的mac m1芯片的笔记本使用ce修改器使用不了
    • ¥15 单相逆变的电压电流双闭环中进行低通滤波PID算法改进
    • ¥15 关于#java#的问题,请各位专家解答!