哈哈，我是刚学R，能不能给个详细的代码啊？那个a[,i]是不是要先定义啊？

这code写的。
为什么upper和lower不一样长阿？为什么比start长呢? 你的数据到底叫什么阿？这个bb1啥的都是列名字吧？
再改改吧

# define likelihood function(normal distribution)
       log.likelihood<-function(para,x){
            y<-sum(-log((1/sqrt(2*pi*para[2]))*exp(-(x-para[1])^2/sqrt(2*pi*para[2]))))
      }
vector<-rnorm(1000,2,0.5);
a<-matrix(vector,100,10);
b<-numeric(20);
b<-matrix(b,2,10);
p<-c(1.8,0.6);
for(i in 1:10){
            
            out<-nlm(log.likelihood,p,a[,i])
            parameter<-out$estimate;
            b[1,i]<-parameter[1];b[2,i]=parameter[2];
            }

看看这个希望对你有帮助，另外
你写的没有定义优化的目标函数，故无法运行

还有一个简单的方法，用apply函数
例如
apply(bb,2,FUNCTION,...)
具体用法参见help