Uhm, it does just not sound right that LsqFit.jl is a solution here, but just a workaround.
Not an expert on numerical analysis but there should be a best practice for data normalization (assuming this is the problem) when the input variables are so different, before throwing them into GLM (ex: subtract mean and divide by standard deviation, or something like that).