Agosto 2017 - R-br

Modelo de regressão envolvendo rasters e vetores numéricos
by Jefferson Ferreira-Ferreira 09 Ago '17

09 Ago '17

Pessoal, Estou tentando ajustar a Estimativa de Máxima Máxima Penalizada de Firth para a regressão logística implementada no pacote logistf para um conjunto de rasters (imagens), mas estou com problemas para implementá-lo usando rasters e vetores numéricos. Para ilustrar o problema: library(raster) library(logistf) # Criando 2 rasters para reproduzir o que estou enfrentando r <- raster(nrow=10, ncol=10) # Variável-raster resposta (binária) com 9 bandas/layers s1 <- lapply(1:9, function(i) setValues(r, sample(0:1,ncell(r),replace = T))) s1 <- stack(s1) # um raster-variável preditora val <- sample(0:60,ncell(r),replace = T) s2 <- raster(nrow=10, ncol=10,vals=val) plot(s1) plot(s2) # Uma segunda variável preditora armazenada como vetor numérico: 9 valores assim como existem 9 bandas/layres em s1 (não é coincidência) exp_2 <- c(27.00,30.02,31.07,32.72,33.73,35.12,35.65,36.06,38.32) Agora, eu quero ajustar o modelo usando o logistf (eu tenho motivos estatísticos para isso, não reproduzido aqui) usando a função calc do pacote raster. É aí que o mistério mora. O raciocínio na formulação do problema é que cada célula/pixel seguiria a segunite fórmula glm: s1/layer1 ~ 27.00 + pixel correspondente em s2 + 27.00:pixel correspondente em s2 s1/layer2 ~ 30.02 + pixel correspondente em s2 + 30.02:pixel correspondente em s2 s1/layer3 ~ 31.07 + pixel correspondente em s2 + 31.07:pixel correspondente em s2 ... e assim por diante para todas as 9 bandas/layers do meu raster-variável resposta s1, que são pareadas com valores de exp_2. Eu tentei algo assim: fun <- function (x) { logistf (x ~ exp_2 + s2 + exp_2:s2)$coeficients } coefs <- calc(s1, fun) Mas estava claro que não funcionaria. A parte complicada é infomar ao script que eu quero que cada valor de exp_2 seja usado com cada camada raster de s1 para este modelo. Qualquer idéia será muito útil. Idéias? Obrigado e abraços. *Jefferson Ferreira-Ferreira, **PhD (abd)* *Geographer* *Ecosystem Dynamics Observatory <http://tscanada.wix.com/ecodyn> - EcoDyn/UNESP* *Department of **Geography * *Institute of Geosciences and Exact Sciences** (IGCE)* *São Paulo State University (UNESP)* *Rio Claro, SP - Brazil*

1 1

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