经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
R语言论坛
R分段函数编写
楼主
慧(会)幸福
9381
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语言论坛
卫生经济学
真实世界经济学(含财经时事)
经管高考
新手入门区
EViews专版
热门文章
我该如何记住你?智能体记忆系统的演化之路
CDA数据分析脱产就业班于2026年3月7日开班! ...
无上高明的“无为”“无住”哲学在传统中国
新宏观丨中美经济总量差距拉大的根源
艾瑞咨询 - 2025年中国早教行业白皮书
第一太平戴维斯 - 2026年中国房地产市场展望 ...
2025中国居民退休准备指数调研报告-清华大学 ...
Measure Theory for Analysis and Probabil ...
现代数学基础19 偏微分方程 孔德兴
科研时间70%耗在“下载-复制-粘贴”?零代码 ...
推荐文章
2026JG学术冬训营:从Stata初高到Python机器 ...
【必看】【本版版规,欢迎发悬赏贴求助】
26年寒假天津站|Gemini论文写作&数据分析 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群