- R-br - DINF UFPR

Adicionar restrições a parâmetros em regressão não
by Wenceslau Teixeira 09 Ago '17

09 Ago '17

> > > Para o codigo que o Walmes disponibilizou para as curvas bimodais de > Durner e preciso instalar o package > minpack.lm install.packages("minpack.lm") > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 8 Aug 2017 12:02:29 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: Delcio Rudinei Bortolanza <delcio.bortolanza(a)outlook.com>, a > lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Adicionar restrições a parâmetros em regressão > não linear no R > Message-ID: > <CAFU=EkYt-YJiO-3gcyRYfW_4uJ+9qyZFr4O5zO+DbAoVhhcjbQ@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Delcio, > > Basta fazer assim: > > durner <- as.formula(theta ~ > w0 * (thetar + ((thetas - thetar)/(1 + (alpha0 * > psi)^n0)^(1 - 1/n0))) + > (1 - w0) * (thetar + ((thetas - thetar)/(1 + > (alpha1 * psi)^n1)^(1 - 1/n1)))) > > dur <- nlsLM(durner, > start = list(w0 = 0.6, > thetas = 0.5, > thetar = 0.25, > alpha0 = 0.1, > n0 = 2, > alpha1 = 0.1, > n1 = 1), > data = dados, > algorithm = "LM", > trace = "F", > control = list(maxiter = 500)) > summary(dur) > > À disposição. > Walmes. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170808/737511e5/attachment-0001.html> > > ------------------------------ > > Message: 2 > Date: Tue, 8 Aug 2017 14:35:29 -0300 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: Cleber Oliveira <oliveiramacedocleber(a)gmail.com>, a lista > Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Análise de dados antes e depois > Message-ID: > <CAKrF98=P-xCSjHHV8_qiCVtAprJryOfQs8kNL+ > HoutQdX0cwiA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Cleber, > > Sua pergunta extrapola o objetivo deste grupo (não CMR, etc.). > > Ademais ela enseja muitas questões antes que possa ser iniciada uma > discussão com assertividade. > > Minha recomendação é que você a refaça nos foruns mais voltados para > questões de Estatística (p.ex.: https://stats.stackexchange.com/) e no > seu > caso Projeto e Análise de Experimentos. > > HTH > -- > Cesar Rabak > > > 2017-08-08 9:50 GMT-03:00 Cleber Oliveira via R-br < > r-br(a)listas.c3sl.ufpr.br > >: > > > Bom dia listeiros (@) > > > > Estou realizando um experimento e de uma única amostra inicial dividi em > > várias subamostras e em cada uma apliquei um tratamento e ao longo dos > dias > > fiz a análise dos parâmetros. Gostaria de verificar o efeito dos > > tratamentos sob as alterações dos parâmetros em relação ao inicial e em > > relação aos tratamentos e tempo. Estou com muitas dúvidas de como > iniciar. > > Alguém que já tenha realizado este tipo de análise poderia me ajudar. > > Obrigado. > > -- > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170808/00668a09/attachment-0001.html> > > ------------------------------ > > Message: 3 > Date: Tue, 08 Aug 2017 19:17:55 -0300 > From: sznelwar(a)uol.com.br > To: walmeszeviani(a)gmail.com, r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Adicionar restrições a parâmetros em regressão > não linear no R > Message-ID: <598a38938f354_38515c87c2813d89477e(a)a4-weasel17.mail> > Content-Type: text/plain; charset="utf-8" > > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170808/703b580c/attachment-0001.html> > > ------------------------------ > > Message: 4 > Date: Tue, 08 Aug 2017 19:21:33 -0300 > From: sznelwar(a)uol.com.br > To: walmeszeviani(a)gmail.com, r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Fwd: Adicionar restrições a parâmetros em > regressão não linear no R > Message-ID: <598a396da05ef_57cd155cbcdd73e448736(a)a4-weasel12.mail> > Content-Type: text/plain; charset="utf-8" > > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170808/a6bfb3b1/attachment-0001.html> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 80, assunto 5 > **************************************** > -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

1 0

Re: [R-br] Adicionar restrições a parâmetros em regressão não linear no R
by sznelwar＠uol.com.br 08 Ago '17

08 Ago '17

1 1

Análise de dados antes e depois
by Cleber Oliveira 08 Ago '17

08 Ago '17

Bom dia listeiros (@) Estou realizando um experimento e de uma única amostra inicial dividi em várias subamostras e em cada uma apliquei um tratamento e ao longo dos dias fiz a análise dos parâmetros. Gostaria de verificar o efeito dos tratamentos sob as alterações dos parâmetros em relação ao inicial e em relação aos tratamentos e tempo. Estou com muitas dúvidas de como iniciar. Alguém que já tenha realizado este tipo de análise poderia me ajudar. Obrigado. --

2 1

Adicionar restrições a parâmetros em regressão não linear no R
by Delcio Rudinei Bortolanza 08 Ago '17

08 Ago '17

Vejam o seguinte script: rm(list=ls(all=TRUE)) library(minpack.lm) dados<-read.table(header = TRUE, dec=",",text= "psi theta 0,01 0,5 0,5 0,48 1 0,42 3 0,4 10 0,39 30 0,35 100 0,3 300 0,28 1000 0,27 ") dados theta=dados$theta psi=dados$psi durner<-as.formula(theta~0.6*(thetar+((thetas-thetar)/(1+(alpha0*psi)^n0)^(1-1/n0)))+ 0.4*(thetar+((thetas-thetar)/(1+(alpha1*psi)^n1)^(1-1/n1)))) dur<-nlsLM(durner, start=list(thetas=0.5,thetar=0.25,alpha0=0.1,n0=2,alpha1=0.1,n1=1),data = dados, algorithm = "LM", trace = "F", control = list(maxiter = 500)) summary(dur) rsquared<-1-(sum(residuals(dur)^2)/sum((theta - mean(theta))^2)) rsquared thetas=coef(dur)[1] thetar=coef(dur)[2] alpha0=coef(dur)[3] n0=coef(dur)[4] alpha1=coef(dur)[5] n1=coef(dur)[6] par(mar=c(4,4,1,1)) par(mfrow = c(1, 1)) plot(psi,theta,xlim = c(min(psi),max(psi)),log = "x",las=1,xaxt="n",ylab=expression(theta~(m^{3}~m^{-3})),xlab=expression(abs(psi)~(kPa)),mgp=c(2.5, 1, 0)) labels <- 10^(-3:3) axis(1, at=labels, labels=labels,lwd.ticks = 2) #,tck=-0.05 axis(1, at=sort(sapply(1:9, function(x) x*10^(-3:3))), label=NA, lwd=0.5) curve(0.6*(thetar+((thetas-thetar)/(1+(alpha0*psi)^n0)^(1-1/n0)))+ 0.4*(thetar+((thetas-thetar)/(1+(alpha1*psi)^n1)^(1-1/n1))),add = T,xname = "psi") Gostaria de, ao invés de dar valores (0.6 e 0.4): durner<-as.formula(theta~0.6*(thetar+((thetas-thetar)/(1+(alpha0*psi)^n0)^(1-1/n0)))+ 0.4*(thetar+((thetas-thetar)/(1+(alpha1*psi)^n1)^(1-1/n1)))) adicionar dois parâmetros (w0 e w1): durner<-as.formula(theta~w0*(thetar+((thetas-thetar)/(1+(alpha0*psi)^n0)^(1-1/n0)))+ w1*(thetar+((thetas-thetar)/(1+(alpha1*psi)^n1)^(1-1/n1)))) mas de tal forma que: w0+w1==1. Alguém tem sabe como poderia fazer isso? Att., Delcio Bortolanza Doutorando em Agronomia-UPF

2 1

Re: [R-br] Digest R-br, volume 80, assunto 1
by Juan C. R. Soto Sotelo 01 Ago '17

01 Ago '17

Concordo que a chamada da query é muito grande, você viu como ficaram os campos na tabela tempo_permanencia_models depois do insert : tipo e tamanho? Juan Error in odbcQuery(channel, query, rows_at_time) : 'Calloc' could not allocate memory (214748364800 of 1 bytes) <https://mailtrack.io/> Enviado com Mailtrack <https://mailtrack.io/install?source=signature&lang=pt&referral=sotojcr@gmai…> Em 1 de agosto de 2017 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > Tópicos de Hoje: > > 1. Re: Erro de memória importando objeto pequeno de SQL > (Cesar Rabak) > 2. Re: Erro de memória importando objeto pequeno de SQL > (Zhu Sha Zang) > > > ---------- Mensagem encaminhada ---------- > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com> > Cc: "a lista Brasileira oficial de discussão do programa R." < > r-br(a)listas.c3sl.ufpr.br> > Bcc: > Date: Mon, 31 Jul 2017 14:53:52 -0300 > Subject: Re: [R-br] Erro de memória importando objeto pequeno de SQL > Correto Pedro! > > Contudo, não entendo porque. . . > > Minha forma de atacar esse problema seria mandar os dados para o BD por > outra maneira (que ñ via R) e tentar conectar apenas para puxar os dados: > aí a gente vê se é na criação do BD via R que ele deixa em aberto um > tamanho default muito grande ou se é a chamada que por alguma coincidência > cria um default exigindo tanta memória (embora nos OSs modernos isso não > deveria gerar problema porque existem os gerenciadores de memória virtual). > > HTH > -- > Cesar Rabak > > > 2017-07-31 10:32 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com>: > >> Cesar, >> >> A função RODBC::sqlQuery não tem exemplo, mas já usei diversas vezes com >> outros scripts e funciona. Esse problema está acontecendo somente no >> exemplo dessa conexão em específico. Então estou entendendo a partir de >> agora que alguma configuração da conexão está pedindo pra reservar memória >> mais do que o a máquina possui. >> >> Pedro Brasil >> >> Em 29 de julho de 2017 20:45, Cesar Rabak <cesar.rabak(a)gmail.com> >> escreveu: >> >>> Pedro, >>> >>> Calloc é uma chamada "primitiva" para o R, uma função da biblioteca C e >>> é uma função feita para "alocar" (reservar) memória e zerá-la. >>> >>> O estranho, para mim, é a tentativa do aplicativo de reservar *tanta* memória >>> assim na chamada da função do R. >>> >>> Essa função sqlQuery() pode ser testada fazendo: >>> >>> > example(sqlQuery) >>> >>> ? >>> >>> Se o exemplo do próprio pacote funcionar, aí teremos que fazer uma busca >>> para ver de onde vem esse comportamento. >>> >>> HTH >>> -- >>> Cesar Rabak >>> >>> >>> >>> >>> 2017-07-28 17:12 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil >>> <emmanuel.brasil(a)gmail.com>: >>> >>>> Ei Cesar, >>>> >>>> Eu tenho pouco familiaridade com SQL, e o amigo da TI que tem e faz >>>> manutenção do servidor diz que o erro é na minha estação com R. Eu acho que >>>> não, mas não tenho certeza porque eu não sei ao certo o que quer dizer esse >>>> erro. O que seria essa "Calloc"? É uma coisa natural do R? >>>> >>>> Abraço, >>>> >>>> Pedro Brasil >>>> >>>> Em 28 de julho de 2017 16:59, Cesar Rabak <cesar.rabak(a)gmail.com> >>>> escreveu: >>>> >>>>> Pedro, >>>>> >>>>> Mesmo sem poder fazer uma tentativa, me parece que o erro é estranho >>>>> demais para não ser notado: a chamada à primitiva da ling. C "calloc" está >>>>> pedindo 200 Gbytes, que falha e me parece algo grande demais para um >>>>> "default" da chamada a sqlQuery(). >>>>> >>>>> HTH >>>>> >>>>> >>>>> 2017-07-28 9:42 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil >>>>> via R-br <r-br(a)listas.c3sl.ufpr.br>: >>>>> >>>>>> Amigos de R, >>>>>> >>>>>> Estou envolvido numa tarefa que importo um banco SQL de um servidor, >>>>>> faço uma modelagem, salvo o modelo no SQL através de uma seriliazação, e >>>>>> depois carrego esse modelo do SQL para fazer previsões em novos dados. >>>>>> >>>>>> Estou seguindo um exemplo disponível nessa página: >>>>>> https://docs.microsoft.com/pt-br/sql/advanced-analytics/tuto >>>>>> rials/walkthrough-deploy-and-use-the-model >>>>>> >>>>>> Sei que esse exemplo não é reproduzível porque depende da conexão com >>>>>> o servidor e credenciais que não posso divulgar. Mas se alguem já passou >>>>>> por esse problema poderia pelo menos indicar aonde está o problema que não >>>>>> estou sabendo identificar? Acho que muito provavelmente seria uma >>>>>> configuração do SQL no servidor... >>>>>> >>>>>> library(RODBC) >>>>>> >>>>>> # define database connection string >>>>>> dbconn <- 'Driver={SQL Server};Server=BLABLABLA;Datab >>>>>> ase=BLABLABLA;Trusted_Connection=yes;' >>>>>> >>>>>> # define a connection >>>>>> conn <- odbcDriverConnect(connection = dbconn) >>>>>> >>>>>> # Estou fazendo o exemplo com o objeto LETTERS só pra ter certeza que >>>>>> é um objeto pequeno e que não deveria ser um problema de memória do R. >>>>>> > LETTERS >>>>>> [1] "A" "B" "C" "D" "E" "F" "G" "H" "I" "J" "K" "L" "M" "N" "O" "P" >>>>>> "Q" "R" "S" "T" "U" "V" "W" "X" >>>>>> [25] "Y" "Z" >>>>>> > teste_pequeno <- serialize(LETTERS, connection = NULL) >>>>>> > teste_pequeno <- paste(teste_pequeno, collapse = "") >>>>>> >>>>>> # o EXEC PersistModel faz >>>>>> # insert into [tempo_permanencia_models] (model, [description]) >>>>>> values (convert(varbinary(max),@m,2), @description) >>>>>> # insere o modelo e uma descricao >>>>>> >>>>>> > p <- paste("EXEC PersistModel @m='", teste_pequeno,"', @description >>>>>> = 'Teste de serialização objeto pequeno'", sep = "") >>>>>> > sqlQuery(conn, p) >>>>>> character(0) >>>>>> > rm(teste_pequeno) >>>>>> > p <- "SELECT model FROM [tempo_permanencia_models] WHERE id_model = >>>>>> 1022" >>>>>> > mod.teste.peq <- sqlQuery(conn, p) >>>>>> Error in odbcQuery(channel, query, rows_at_time) : >>>>>> 'Calloc' could not allocate memory (214748364800 of 1 bytes) >>>>>> > p <- paste0("SELECT id_model, description FROM >>>>>> [tempo_permanencia_models] WHERE id_model = 1022") >>>>>> > sqlQuery(conn, p) # >>>>>> id_model description >>>>>> 1 1022 Teste de serialização objeto pequeno >>>>>> >>>>>> Reparem que o objeto está salvo no servidor. Foi salvo a partir do R, >>>>>> mas não retorna para o R por conta do erro de memoria. Eu não estou sabendo >>>>>> nem por onde começar. >>>>>> >>>>>> Que a força esteja sempre com vocês e abraço forte. >>>>>> >>>>>> Pedro Brasil >>>>>> >>>>>> _______________________________________________ >>>>>> R-br mailing list >>>>>> R-br(a)listas.c3sl.ufpr.br >>>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>>> código mínimo reproduzível. >>>>>> >>>>> >>>>> >>>> >>> >> > > > ---------- Mensagem encaminhada ---------- > From: Zhu Sha Zang <zhushazang(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Cc: > Bcc: > Date: Mon, 31 Jul 2017 18:52:05 -0400 > Subject: Re: [R-br] Erro de memória importando objeto pequeno de SQL > > Parece que RODBC está "deprecated", tente usar o pacote RSQLServer. > > > Atenciosamente > > On 07/31/2017 01:53 PM, Cesar Rabak via R-br wrote: > > Correto Pedro! > > Contudo, não entendo porque. . . > > Minha forma de atacar esse problema seria mandar os dados para o BD por > outra maneira (que ñ via R) e tentar conectar apenas para puxar os dados: > aí a gente vê se é na criação do BD via R que ele deixa em aberto um > tamanho default muito grande ou se é a chamada que por alguma coincidência > cria um default exigindo tanta memória (embora nos OSs modernos isso não > deveria gerar problema porque existem os gerenciadores de memória virtual). > > HTH > -- > Cesar Rabak > > > 2017-07-31 10:32 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com>: > >> Cesar, >> >> A função RODBC::sqlQuery não tem exemplo, mas já usei diversas vezes com >> outros scripts e funciona. Esse problema está acontecendo somente no >> exemplo dessa conexão em específico. Então estou entendendo a partir de >> agora que alguma configuração da conexão está pedindo pra reservar memória >> mais do que o a máquina possui. >> >> Pedro Brasil >> >> Em 29 de julho de 2017 20:45, Cesar Rabak <cesar.rabak(a)gmail.com> >> escreveu: >> >>> Pedro, >>> >>> Calloc é uma chamada "primitiva" para o R, uma função da biblioteca C e >>> é uma função feita para "alocar" (reservar) memória e zerá-la. >>> >>> O estranho, para mim, é a tentativa do aplicativo de reservar *tanta* memória >>> assim na chamada da função do R. >>> >>> Essa função sqlQuery() pode ser testada fazendo: >>> >>> > example(sqlQuery) >>> >>> ? >>> >>> Se o exemplo do próprio pacote funcionar, aí teremos que fazer uma busca >>> para ver de onde vem esse comportamento. >>> >>> HTH >>> -- >>> Cesar Rabak >>> >>> >>> >>> >>> 2017-07-28 17:12 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil >>> <emmanuel.brasil(a)gmail.com>: >>> >>>> Ei Cesar, >>>> >>>> Eu tenho pouco familiaridade com SQL, e o amigo da TI que tem e faz >>>> manutenção do servidor diz que o erro é na minha estação com R. Eu acho que >>>> não, mas não tenho certeza porque eu não sei ao certo o que quer dizer esse >>>> erro. O que seria essa "Calloc"? É uma coisa natural do R? >>>> >>>> Abraço, >>>> >>>> Pedro Brasil >>>> >>>> Em 28 de julho de 2017 16:59, Cesar Rabak <cesar.rabak(a)gmail.com> >>>> escreveu: >>>> >>>>> Pedro, >>>>> >>>>> Mesmo sem poder fazer uma tentativa, me parece que o erro é estranho >>>>> demais para não ser notado: a chamada à primitiva da ling. C "calloc" está >>>>> pedindo 200 Gbytes, que falha e me parece algo grande demais para um >>>>> "default" da chamada a sqlQuery(). >>>>> >>>>> HTH >>>>> >>>>> >>>>> 2017-07-28 9:42 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil >>>>> via R-br <r-br(a)listas.c3sl.ufpr.br>: >>>>> >>>>>> Amigos de R, >>>>>> >>>>>> Estou envolvido numa tarefa que importo um banco SQL de um servidor, >>>>>> faço uma modelagem, salvo o modelo no SQL através de uma seriliazação, e >>>>>> depois carrego esse modelo do SQL para fazer previsões em novos dados. >>>>>> >>>>>> Estou seguindo um exemplo disponível nessa página: >>>>>> https://docs.microsoft.com/pt-br/sql/advanced-analytics/tuto >>>>>> rials/walkthrough-deploy-and-use-the-model >>>>>> >>>>>> Sei que esse exemplo não é reproduzível porque depende da conexão com >>>>>> o servidor e credenciais que não posso divulgar. Mas se alguem já passou >>>>>> por esse problema poderia pelo menos indicar aonde está o problema que não >>>>>> estou sabendo identificar? Acho que muito provavelmente seria uma >>>>>> configuração do SQL no servidor... >>>>>> >>>>>> library(RODBC) >>>>>> >>>>>> # define database connection string >>>>>> dbconn <- 'Driver={SQL Server};Server=BLABLABLA;Datab >>>>>> ase=BLABLABLA;Trusted_Connection=yes;' >>>>>> >>>>>> # define a connection >>>>>> conn <- odbcDriverConnect(connection = dbconn) >>>>>> >>>>>> # Estou fazendo o exemplo com o objeto LETTERS só pra ter certeza que >>>>>> é um objeto pequeno e que não deveria ser um problema de memória do R. >>>>>> > LETTERS >>>>>> [1] "A" "B" "C" "D" "E" "F" "G" "H" "I" "J" "K" "L" "M" "N" "O" "P" >>>>>> "Q" "R" "S" "T" "U" "V" "W" "X" >>>>>> [25] "Y" "Z" >>>>>> > teste_pequeno <- serialize(LETTERS, connection = NULL) >>>>>> > teste_pequeno <- paste(teste_pequeno, collapse = "") >>>>>> >>>>>> # o EXEC PersistModel faz >>>>>> # insert into [tempo_permanencia_models] (model, [description]) >>>>>> values (convert(varbinary(max),@m,2), @description) >>>>>> # insere o modelo e uma descricao >>>>>> >>>>>> > p <- paste("EXEC PersistModel @m='", teste_pequeno,"', @description >>>>>> = 'Teste de serialização objeto pequeno'", sep = "") >>>>>> > sqlQuery(conn, p) >>>>>> character(0) >>>>>> > rm(teste_pequeno) >>>>>> > p <- "SELECT model FROM [tempo_permanencia_models] WHERE id_model = >>>>>> 1022" >>>>>> > mod.teste.peq <- sqlQuery(conn, p) >>>>>> Error in odbcQuery(channel, query, rows_at_time) : >>>>>> 'Calloc' could not allocate memory (214748364800 of 1 bytes) >>>>>> > p <- paste0("SELECT id_model, description FROM >>>>>> [tempo_permanencia_models] WHERE id_model = 1022") >>>>>> > sqlQuery(conn, p) # >>>>>> id_model description >>>>>> 1 1022 Teste de serialização objeto pequeno >>>>>> >>>>>> Reparem que o objeto está salvo no servidor. Foi salvo a partir do R, >>>>>> mas não retorna para o R por conta do erro de memoria. Eu não estou sabendo >>>>>> nem por onde começar. >>>>>> >>>>>> Que a força esteja sempre com vocês e abraço forte. >>>>>> >>>>>> Pedro Brasil >>>>>> >>>>>> _______________________________________________ >>>>>> R-br mailing list >>>>>> R-br(a)listas.c3sl.ufpr.br >>>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>>> código mínimo reproduzível. >>>>>> >>>>> >>>>> >>>> >>> >> > > > _______________________________________________ > R-br mailing listR-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >

1 0

Aggregate com função de dois argumentos
by Marcelo Santana Camacho 27 Jul '17

27 Jul '17

Bom dia, Alguém poderia me dar uma dica com esse problema? Estou tentando agregar um conjunto de dados, calculando o gini (pacote reldist), que recebe a variável principal e uma variável de ponderação - gini(x, weights=rep(1,length=length(x))). Minha dificuldade é passar os dois argumentos para a função. library(reldist) ineq_df<-read.csv("amostra_input.csv") aggregate(c(qtd,renda)~cod_setor,data=ineq_df,gini) -- *__________________________________Marcelo Santana Camacho*UNIFESSPA

2 2

Alterações no Plot 'extRemes'
by Bruce Kelly 25 Jul '17

25 Jul '17

Bom dia listeiros(@) Estou utilizando o pacote 'extRemes', minha duvida é como posso melhorar o Plot, já que ele disponibiliza 4 figuras em uma só. Como posso alterar legendas, inserir valores medios, etc ou trabalhar com cada gráfico separado. codigo AE2=fevd(EST1, dad1, type=c('GEV'), method=c('Bayesian'), time.units='months', period.basis='year') postmode(AE2) plot(AE2) plot(AE2, 'trace') ci(AE2) ci(AE2,type='parameter') [cid:9141c91f-0af9-4dec-8f9d-c83d492e9598] Bruce Kelly da N. Silva Dr. em Ciências Climáticas - UFRN Pesquisador - PPGCC/UFRN Lattes: http://lattes.cnpq.br/4323502426271751

1 0

tratamento de pontos discrepantes (pesos de bovinos)
by Adriele Giaretta Biase 24 Jul '17

24 Jul '17

Boa tarde pessoal, estou trabalhando com um banco de dados que possuem várias coletas de pesos ao longo do tempo de vários bois. O meu objetivo é deixar um ajuste simples (pode até ser um modelo não- linear ao invés de uma regressão simples), porém, preciso remover os possíveis pontos discrepantes que atrapalham o ajuste, de forma recursiva (automática). O algoritmo será executado em vários momentos, enquanto acompanha o crescimento do animal. O mais importante seria a remoção dos *outiliers* do banco de dados. Eu criei um algoritmo para fazer isso. Gostaria de saber se alguém possui uma ideia melhor ou sugestão. Segue a função criada. *OBS:* Muitos pontos discrepantes ainda ficam no banco, isso devido a um erro operacional da balança que pesa os animais. Precisaria de algo mais robusto para remover os pontos destoantes de forma recursiva para todos os animais. DADOS é um data.frame com as seguintes colunas: BRINCO (identificação ou numero do animal), DATA (dia da pesagem do peso do animal) e PESO (peso do animal correspondente ao dia) fit_linear <- function(DADOS){ fit_aux = c() BRINCO = as.numeric(as.vector(levels( as.factor(DADOS[,"BRINCO"])))) pdf(file = "ajuste%03d.pdf") # salvar os gráficos em pdf for (i in BRINCO){ dados = subset(DADOS, BRINCO == i, selec = c( BRINCO, DATA, PESO)) j=1 while (j<=30){ dias = as.vector(as.numeric(as.factor(dados[,"DATA"]))) pesos = as.numeric(as.vector(dados$PESO)) dad = cbind(dias, pesos) mod <- lm(pesos~dias) INFLUENTES = influence.measures(mod) dados_analise = cbind(INFLUENTES$is.inf, dad) dados_influentes = as.data.frame(dados_analise) colnames(dados_influentes) = c("x1", "x2", "x3", "x4", "x5", "x6", "DATA", "PESO") dados = subset(dados_influentes, (x1==FALSE & x2==FALSE & x3==FALSE & x4==FALSE & x5==FALSE & x6==FALSE), selec = c(x1, x2, x3, x4, x5, x6, DATA, PESO)) j = j+1 } res <- residuals(mod) h <- hatvalues(mod) P = length(mod$coefficient); N=length(peso); P; N; hc<-3*P/N;hc limt <- list(c(DFB=2/sqrt(N),FDFits=2*sqrt(P/N),cov.r=3*P/N, Cook=qf(0.5,2,8, lower.tail = TRUE, log.p = FALSE),hat=3*P/N)) N=length(peso); P=length(mod$coefficient) rs <- rstudent(mod) h <- lm.influence(mod)$hat; lc <- 3*P/N minrs=min(min(rs),-3) maxrs=max(max(rs),3) ymin=minrs-.1 ymax=maxrs+.1 maxh=max(max(h),lc) minh=min(h) xmin=minh-.1 xmax=maxh+.1 par(mfrow=c(1,1)) plot(c(xmin,xmax),c(ymin,ymax), type="n", xlab="h - leverage", ylab="RStudent", main= i) abline(h=-2.5, col="red") abline(h=2.5,col="red"); abline(v=lc, col="blue") points(h,rs) plot(dad, main= i) abline(mod,lty=2) } dev.off() } Agradeço, -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

3 2

Re: [R-br] tratamento de pontos discrepantes (pesos de bovinos)
by sznelwar＠uol.com.br 23 Jul '17

23 Jul '17

1 0

geom_boxplot
by fabianidenise＠gmail.com 21 Jul '17

21 Jul '17

Bom dia listeiros, estou com uma dúvida, e talvez alguém possa me indicar alguma maneira de resolver. Estou com uma série histórica de dados diários de chuva, e preciso fazer um boxplot dos valores agregados mensalmente. DATA CHUVA 01/01/1980 5 02/01/1980 0 03/01/1980 0 04/01/1980 0 05/01/1980 0 06/01/1980 0 07/01/1980 4,9 08/01/1980 34 09/01/1980 5,6 10/01/1980 8 11/01/1980 0 . . . Até aí td bem, consegui executar pelas linhas de comando abaixo: P1 <- ggplot(dados,aes(format(as.Date(DATA), "%m"), CHUVA)) + geom_boxplot() P1 No entanto, como para a variável chuva esses valores diários não fazem muito sentido para a minha análise, preciso fazer um boxplot aggregando esses valores pela soma mensal. Alguém teria alguma dica de como posso fazer isso? Muito obrigada, Fabiani _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

2 1