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(yvals~
soi.std+nao.std+np.std+tvals[1:50],family=quasi(link="log",var="mu"))
return(mod$coef)
}