- 1、本文档共21页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
R语言入门_十三章、十四章课件
第十三章
调试;案例:寻找连续为1的游程
findruns1-function(x,k){
n-length(x)
runs-vector(length=n)
count-0
for(i in 1:(n-k+1)){
if(all(x[i:i+k-1]==1)){
count-count+1
runs[count]-i
}
}
if(count0){
runs-runs[1:count]
}
else runs-NULL
return(runs)
}
findruns1(c(1,0,0,1,1,0,0,1,1,1),2)
[1] 3 4 7 8 9
显然,出错了;调试方法一:手动,在程序中添加多个print(),把程序的内部细节输出来查看是哪里出的问题。
function(x,k){
print(x) #输出x
print(k) #输出k
n-length(x)
print(n) #输出n
runs-vector(length=n)
count-0
for(i in 1:(n-k+1)){
if(all(x[i:i+k-1]==1)){
count-count+1
runs[count]-i
print(runs) #输出循环的每一步
}
}
if(count0){
runs-runs[1:count]
print(runs) #输出结果
}
else runs-NULL
return(runs)
};findruns1(c(1,0,0,1,1,0,0,1,1,1),2)
[1] 1 0 0 1 1 0 0 1 1 1 #正确
[1] 2 #正确
[1] 10 #正确
[1] 3 0 0 0 0 0 0 0 0 0 #出错!!
[1] 3 4 0 0 0 0 0 0 0 0
[1] 3 4 7 0 0 0 0 0 0 0
[1] 3 4 7 8 0 0 0 0 0 0
[1] 3 4 7 8 9 0 0 0 0 0
[1] 3 4 7 8 9
[1] 3 4 7 8 9
可以看出,从把实参带入形参,到建立长度等于x长度d的向量,都是正确的。错误开始于循环的第一步。因此我们应该去检查循环语句,不难发现,问题出在all(x[i:i+k-1]==1)中,我们没有把i+k-1扩起来。;调试方法二:使用R的调试工具
在代码中插入print()显然比较繁琐,特别是程序较长的时候,使用R的调试工具会更便捷。
R的核心调试工具有浏览器构成,它可以让你逐行运行代码,并在运行过程中检查,可以用dubug(f)或browser()打开这个浏览器。
debug(f)函数可以把函数f()设置成调试状态,意味着每次调用f()都会进入这一状态,取消这一状态需调用undebug(f)。在2.10版本之后,可以用debugonce()代替(检查漏洞一次)。; debugonce(findruns1)
findruns1(c(1,0,0,1,1,0,0,1,1,1),2)
然后开始调试….
Browse[2] x
[1] 1 0 0 1 1 0 0 1 1 1
Browse[2]
debug at findruns1.R#2: n - length(x)
Browse[2]
debug at findruns1.R#3: runs - vector(length = n)
Browse[2] print(n)
[1] 10
Browse[2]
debug at findruns1.R#4: count - 0
Browse[2]
debug at findruns1.R#5: for (i in 1:(n - k + 1)) {
if (all(x[i:i + k - 1] == 1)) {
count - count + 1
runs[count] - i
}
};Browse[2]
debug at findruns1.R#5: i
Browse[2]
debug at findruns1.R#6: if (all(x[i:i + k - 1] == 1)) {
c
您可能关注的文档
- Rare allele of OsPPKL1 associated with grain length causes课件.ppt
- Rapamycin Induces MAP Kinase Phosphatase (MKP)2-课件.ppt
- Ramtron FRAM 介绍课件.ppt
- RAS单服务器实施配置20070816课件.ppt
- RAS远程访问服务器配置课件.ppt
- Razvojno-podporni center za odprto kodo课件.ppt
- RAS远程访问服务课件.ppt
- RBEs and MPCs in MSC.Nastran2课件.ppt
- RC-135侦察机课件.ppt
- RBS2206V2介绍课件.ppt
文档评论(0)