if you use functions like lm and glm inside a function, you may run
into scoping problems. Use the assign function to deal with this.
For example,
>myfun_function(rowvals){
assign(yvals,rowvals[3:52],0)
mod_glm(yvalssoi.std+nao.std+np.std+tvals[1:50],family=quasi(link=log,var=mu))
return(mod$coef)
}