R如何对变量自导

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

p1<-vector(length=9)
p2<-vector(length=9)
v1<-vector(length=9)
v2<-vector(length=9)
v1[10]=0
v2[10]=0
p2[9]=3
for(i in 9:1)
  v2=p2+v1[i+1]-p2*((1-p2)/7)
#在这里我想用v2对p2求导，然后令倒数等于0，怎么弄啊？？？求大神
  v1=0.8*(1-((1+p1)/10))*(p1+v1[i+1])+0.8*((1+p1)/10)*v2+0.2*v1[i+1]
#这里也是差不多的意思

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：如何

1. func <- function(p2) {
   
2.     eval(D(expression(p2 + v1 - p2*((1 - p2)/7)), "p2"))
   
3. }
   
4. uniroot(func, c(-1e2, 1e2))

复制代码