Проблемы с построением экспоненциальной формулы - PullRequest
Новая
       59

Проблемы с построением экспоненциальной формулы

0 голосов
/ 21 июня 2020

Я пытаюсь построить экспоненциальный график переменной. Коэффициент переменной очень высок (350 миллионов) по результатам GLM. С другими переменными, у которых были более низкие коэффициенты, я смог легко построить их график без проблем. Я пытался установить интервал последовательности все меньше и меньше, но он продолжает давать сбой r, когда я пытаюсь построить его. Какие-либо предложения? Я уже безуспешно пытался разбить данные. Мои векторы тоже очень большие (18Мб). 

chlautcnod<-seq.int(0, 2.45259, 0.000001)

chlautcnodline<- glmnodosaALL$coefficients[1] + 
glmnodosaALL$coefficients[2]*mean(bornodosaAP$Chl_spring) + 
glmnodosaALL$coefficients[3]*chlautcnod + glmnodosaALL$coefficients[4]*mean(bornodosaAP$Dist_coast) +
glmnodosaALL$coefficients[5]*mean(bornodosaAP$Chl_winter)+ glmnodosaALL$coefficients[6]*mean(bornodosaAP$Depth) +
glmnodosaALL$coefficients[7]*mean(bornodosaAP$Chl_yr_avg)+ glmnodosaALL$coefficients[8]*mean(bornodosaAP$Dist_complete_river) +
glmnodosaALL$coefficients[9]*mean(bornodosaAP$Temp_yr_min)+ glmnodosaALL$coefficients[10]*mean(bornodosaAP$Chl_summer)+
glmnodosaALL$coefficients[11]*mean(bornodosaAP$Chl_yr_max)+ glmnodosaALL$coefficients[12]*mean(bornodosaAP$SWH_summer)+
glmnodosaALL$coefficients[13]*mean(bornodosaAP$SWH_yr_min)+ glmnodosaALL$coefficients[14]*mean(bornodosaAP$SWH_spring)

gc(plot(exp(1)^chlautcnodline~chlautcnod, xlab = (expression(paste("Chlorophyll-α Autumn (mg/m"^"3"~")"))), ylab= "Probability of C. nodosa occurance",ylim=c(0,0.05),xlim=c(0.15,0.17), type="l", bty="l")

1

...