有的时候在处理数据时,在处理某个函数的极小值的时候,用optimize函数比较好用,
用法是optimize(函数公式,interval=c()包含极小值区间,tol=精度)
详见例子:
> x<-1:100 #首先产生x的取值
> x
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
[19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
[37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
[55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
[73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
[91] 91 92 93 94 95 96 97 98 99 100
> f<-function(p) sum(log(1+(x-p)^2)) #然后建立包含P的公式
> f
function(p) sum(log(1+(x-p)^2))
> out<-optimize(f,c(0,20)) #使用optimize求在区间0,20间的极小值
> out
$minimum
[1] 19.99993
$objective
[1] 628.5662
> out<-optimize(f,c(20,40)) #求20,40间的极小值
> out
$minimum
[1] 39.99993
$objective
[1] 593.0818