Estou tentando gerar uma função geral, que rode vários modelos de regressão. O problema é que, quando um modelo dá erro, a função principal para. Eu gostaria de usar uma condicional que caso tenha erro, pule para a próxima etapa:
Exemplo: ao tentar rodar um modelo nls()
, apresenta essa mensagem descrita abaixo, como essa função é uma das que estou testando dentro de uma função geral, ele trava a função neste ponto.
um exemplo, quando eu quero carregar algum pacote, eu normalmente faço assim:
if(!require(nome do pacote)){install.packages("nome do pacote")}
no meu intendimento é assim: se(não existe esse pacote){ instale o pacote} por isso eu aqui informo que "!" representaria o "erro".
supondo que eu crie uma função:
melhor.ajuste<-function(df,x,y){
n1<-nls(y~x+c/e*x^d,start=....,data=df)
n2<-nls(y~x-b*e^d/log(x),start=....,data=df)
n3<-nls(y~x+a*b/c^x,start=....,data=df)
list(modelo1=n1,modelo2=n2,modelo3=n3)}
caso por exemplo a função n2 apresente uma erro qualquer, como exemplo:
stop(simpleError("Error in nlsModel(formula, mf, start, wts) :\n matriz gradiente singular com estimativas de parâmetros iniciais"))
minha função não me apresentar valor nenhum, pois houve um erro!
o que eu desejo era algo do tipo:
if(!nls(y~x+c/e*x^d,start=....,data=df)){n1=0}else{n1<-nls(y~x+c/e*x^d,start=....,data=df)}
if(!nls(y~x+a*b/c^x,start=....,data=df)){n2=0}else{n1<-nls(y~x+a*b/c^x,start=....,data=df)}
if(!nls(y~x+a*b/c^x,start=....,data=df)){n3=0}else{n1<-nls(y~x+a*b/c^x,start=....,data=df)}
alguém sabe de alguma forma de eu fazer isso! o excel faz, pela função SEERRO(), onde seerro(caso tenha erro; faça isso), no meu caso seria, se(seerro(nls(y~x+a*b/c^x,start=....,data=df);0)==0;n1=0;n1=nls(y~x+a*b/c^x,start=....,data=df))
, ou seja seria uma condicional que identifica erro, dentro de uma condicional.
Deixando claro, que em nenhum momento, meu principal objetivo é descobrir qual é o erro e tentar concertar, pois eu irei aplicar a função em mais de 50 modelos, sabendo que alguns não se ajustam mesmo, por ser matriz singular, por estourar o minFactor, e ou coisa parecida. meu objetivo, é encontrar uma forma de que caso a função apresente erro ele concatene 0 ao objeto, e passe para o proximo modelo.
e ao final me forneça a lista de todos os modelos.
try
?