经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
R语言论坛
R分段函数编写
楼主
慧(会)幸福
9280
5
收藏
2014-02-16
我的代码如下,为什么运行完之后画出来的图不是我想要的,理论上讲画出来的图应该是两边是水平线,中间是光滑曲线,
但是最终画出来是条直线,求高手帮忙检查!多谢了!
a=-0.618034
b=1.618034
l=function(x){
s=2*log(1+x^2)
t=2*log(1+(x-1)^2)
s-t
}
de=function(x){
4*x/(1+x^2)-4*(x-1)/(1+(x-1)^2)
}
lstar=function(x){
if ((x > a)&(x < b)){
result=l(x)
}
if (x<=a){
result=de(a)*(x-a)+l(a)
} else {
result=de(b)*(x-b)+l(b)
}
result
}
x=seq(from=-4,to=4,by=0.01)
y=lstar(x)
plot(x,y,type='l',ylab= "log-likelihood ratio")
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
全部回复
沙发
yywan0913
2014-2-17 17:02:00
复制代码
因为x是向量形式,所以你的function不对,想办法将function中得x定义为向量形式,然后再比
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
藤椅
yywan0913
2014-2-17 17:04:01
此处的ifelse就相当于那么多判断了
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
板凳
ntsean
2014-2-18 08:44:50
你的lstar函数可以简化,并且vectorize
lstar <- function(x) {
x0 <- pmin(pmax(a, x), b)
de(x0)*(x-x0) + l(x0)
}
这样就可以了
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
报纸
ntsean
2014-2-18 08:47:04
或者你可以直接用 y = sapply(x, lstar) 也是可以的
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
地板
慧(会)幸福
2014-2-20 20:44:04
ntsean 发表于 2014-2-18 08:47
或者你可以直接用 y = sapply(x, lstar) 也是可以的
非常感谢
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
相关推荐
想不通的分段函数; 没搞懂的核光滑>///<
请教:多个报价曲线加总为一条报价曲线问题
求助:分段函数的拟合问题?
求助匹配问题,急啊
【求助】一点程序问题。。
请问R可以做分段函数吗?
R语言关于循环跳出的问题
R中画简单的分段函数图形
R语言求分段函数
分段函数的计量建模。
栏目导航
R语言论坛
商学院
教师之家与经管教育
爱问频道
经管文库(原现金交易版)
文献求助专区
热门文章
你的SSCI发表焦虑,AI真的能懂吗?——一篇 ...
CDA数据分析脱产就业班于2025年08月02日开班 ...
CDA考试模拟题库:新增章节练习题(更新于1 ...
【AI Agent可靠性】 智能体Agent记忆系统: ...
全球数字经贸规则年度观察报告(2025年)
《以日为鉴》阅读分享
2025骑行配件出海研究报告
Modern Computer Algebra
河南话破圈:从“土味方言”到“潮流符号” ...
Matrix Algebra-James E·Gentle
推荐文章
AI狂潮席卷学术圈,不会编程也能打造专属智 ...
10月重磅来袭|《打造Coze/Dify专属学术智能 ...
最快1年拿证,学费不足5W!热门美国人工智能 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群