请教程序

请教一个R程序问题
给定初始P值；
While(|b1-a1|+…+|b6-a6|>0.001)
{
fit=function f(x, p);
fit1=[a1, a2,a3,a4,a5,a6]
b1=a1,b2=a2,b3=a3,b4=a4,b5=a5,b6=a6; #保存计算的六个参数
fit2=function g(x, a1,a2,a3,a4,a5,a6).###利用fit1计算的六个参数来计算p值
p=fit2
}
完成上面的循环后，保留最后一次的结果：即p, a1, a2,a3,a4,a5,a6
输出p, a1, a2,a3,a4,a5,a6

1. fitt <- function(x,p){
   
2. 
   
3.     return(p)
   
4. }
   
5. 
   
6. gfit <- function(x,p){
   
7.     return(p)
   
8. }
   
9. 
   
10. dista <- function(pa,pb){
    
11.     dista <- sum(abs(pa-pb))
    
12.     return(dista)   
    
13. }
    
14. output <- numeric(7)
    
15. pa <- numeric(6)
    
16. repeat{
    
17.     p <- fitt(x,pa)
    
18.     pb <- gfit(x,pa)
    
19.     if(dista(pa,pb)>0.001) {
    
20.         output <<-c(p,pa)
    
21.         break
    
22.     }
    
23.     pa <- pb
    
24. }

复制代码

nuomin 发表于 2015-4-9 10:45

十分感谢