现在做模拟时遇到很大的问题，我的模型是时间序列模型，参数的估计值需要用条件最小二乘得到，所有R中关于lasso的包我用不了，所以只能自己编写程序，我的编程思路是用范剑青2001年SCAD的那篇论文里提到的用局部二次逼近的方法，近似出惩罚函数，然后再用迭代法求得估计值，这个编程我已经卡了很久了，有哪位对这种编程熟悉的，请指点一下，有偿指点也可以！！！