Novembro 2014 - R-br

Resolução de grafico word - pdf
by Leonardo Monteiro 26 Nov '14

26 Nov '14

Pessoal, não é uma pergunta relacionada ao R, mas eu estou em fase de termino da tese, e quando passo o work para pdf, os graficos e mapas perdem mto a qualidade (Fica com aparencia de borrado). Eu estou gerando ambos em .tiff com 300 dpi. Por favor, se alguem tiver uma sugestão, eu agradeço. Leonardo

6 5

Facebook
by Vinicius Brito Rocha 25 Nov '14

25 Nov '14

Pessoal, ninguém nunca extraiu dados do facebook com o R? Perguntei mas ninguém respondeu abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário * *M.Sc. Engenharia de Produção/PO* *Ph.D. Estudante Sistemas Computacionais*

6 7

Modelo misto - Coeficiente parcial de determinação
by Fernando Souza 25 Nov '14

25 Nov '14

caros amigos, Estou tentando calcular o coeficiente parcial de determinação, para um efeito de modelo misto, conforme recomendação do livro do Kutner et al 2004 (Applied linear statistical models, 5º edição; para edições anteriores a 4º o autor principal era o Neter). Nesse livro ele calcula ess R2parcial com: R2parcial = (SSE(X1)-SSE(X1,X2)) / SSE(X1) ; onde SSE(X1) é a soma de quadrado do erro do modelo reduzido (modelo sem a variável de interesse) e SSE(X1,X2) a soma de quadrado mo modelo completo (contendo todas as variáveis) Esse R2 mede a porcentagem de redução da variação da variável resposta quando a variável de interesse é adicionada ao modelo. Ele propõe para modelo linear generalizado. Desconfio que para modelo misto não seja realizado da mesma forma devido ao modelo aleatório. Na função abaixo utilizei a função resid() para obtenção dos resíduos e realizar as soma de quadrado, no entanto esse resíduo é resultado do modelo de efeito fixo e aleatório. Há uma forma de remover o resíduo do modelo fixo dos modelos por meio do comando : modelo.lme$residuals[,1]. minha dúvida é como eu devo me proceder para obter esta estatística. R2parcial <- function(modelocompleto,modeloreduzido){ SSR.compl<-sum(resid(modelocompleto)^2) SSR.reduz<- sum(resid(modeloreduzido)^2) R2P <- (SSR.reduz-SSR.compl)/SSR.reduz return(paste("R2 parcial:",round(R2P,3))) } #um pequeno exemplo. somente para demonstração, não considerei os ajustes necessários a esses dados. library(nlme) #o banco de dados Gasoline, está presente no pacote nlme modelocompleto <- lme(yield~endpoint+API+vapor+ASTM,random=~1|Sample,data=Gasoline) modeloreduzido <-lme(yield~API+vapor+ASTM,random=~1|Sample,data=Gasoline) # reduçao gerada pela variavel endpoint R2parcial(modelocompleto,modeloreduzido)

1 0

Re: [R-br] Diferença entre métodos
by Ari Clecius 25 Nov '14

25 Nov '14

Boa noite, alguém poderia me ajudar o porque da diferença? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 Em Sábado, 22 de Novembro de 2014 10:06, Ari Clecius <ari072000(a)yahoo.com.br> escreveu: Bom dia, estou rodando um experimento do tipo Box Behnken e usei a função reg e comparei com o obtido pela função rsm, o problema é que o que os dois consideram significativo difere, aonde posso estar errando? Testei também a função AIC para o modelo que ele fornece e dá diferente pela função AIC(). Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CMR library(DoE.wrapper)Design.1_2<- bbd.design(nfactors= 3 , ncenter= 4 ,randomize= FALSE,seed=7480, factor.names=list( A=c(0,400),B=c(600,2200),C=c(1,10) ) ) res=c(73.450,75.725,63.980,67.670,76.070,77.345,66.005,67.335,84.740,81.595,74.830,70.325,72.460,73.170,73.925,73.800)length(res)#Medidas resumoDesign.1_2<- add.response(Design.1_2,res, replace=FALSE)plancoded<- code.design(Design.1_2)plancoded LinearModel.1 <- aov(res ~(x1+x2+x3)^2+ I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)sreg = step(LinearModel.1,direction="both") #Melhor modelo#Step: AIC=14.07#res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2) # Df Sum of Sq RSS AIC#<none> 16.074 14.074#+ x1:x2 1 0.501 15.574 15.568#+ x2:x3 1 0.462 15.612 15.607#+ x1:x3 1 0.001 16.073 16.073#- x1 1 9.181 25.255 19.303#- I(x2^2) 1 9.310 25.384 19.385#- I(x3^2) 1 36.195 52.269 30.941#- x2 1 79.223 95.297 40.551#- I(x1^2) 1 86.793 102.867 41.774#- x3 1 212.747 228.821 54.566 summary(sreg) # Df Sum Sq Mean Sq F value Pr(>F) #x1 1 9.18 9.18 5.140 0.04959 * #x2 1 79.22 79.22 44.357 9.27e-05 ***#x3 1 212.75 212.75 119.118 1.72e-06 ***#I(x1^2) 1 86.79 86.79 48.596 6.53e-05 ***#I(x2^2) 1 9.31 9.31 5.213 0.04831 * #I(x3^2) 1 36.20 36.20 20.266 0.00149 ** #Residuals 9 16.07 1.79 #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Design.1_2.coded<-plancodedrsmModel.2 <- rsm(res ~ FO(x1, x2, x3) + TWI(x1, x2, x3) + PQ(x1, x2, x3), data=Design.1_2.coded)summary(rsmModel.2) # Estimate Std. Error t value Pr(>|t|) #(Intercept) 73.33875 0.79347 92.4275 1.081e-10 ***#x1 1.07125 0.56107 1.9093 0.104810 #x2 -3.14688 0.56107 -5.6087 0.001370 ** #x3 -5.15688 0.56107 -9.1911 9.350e-05 ***#x1:x2 0.35375 0.79347 0.4458 0.671360 #x1:x3 0.01375 0.79347 0.0173 0.986736 #x2:x3 -0.34000 0.79347 -0.4285 0.683246 #x1^2 -4.65813 0.79347 -5.8705 0.001081 ** #x2^2 1.52563 0.79347 1.9227 0.102878 #x3^2 3.00813 0.79347 3.7911 0.009062 ** #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Multiple R-squared: 0.9664, Adjusted R-squared: 0.916 #F-statistic: 19.17 on 9 and 6 DF, p-value: 0.0009315 LinearModel.2 <- aov(res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)summary(LinearModel.2 )AIC(LinearModel.2)#[1] 61.47997 #valor obtido anteriormente Step: AIC=14.07

1 0

Ler arquivo .xls - RDOBC e gdata
by Andre Oliveira 24 Nov '14

24 Nov '14

Bom dia,sempre usei o gdata no linux e o RDOBC e gdata no win para ler arquivos .xls, no entanto venho tendo problemas para colocar o gdata e o xlsx mesmo no linux, e quando passo algum trabalho para alunos sempre dá problemas com rjava etc etc. Diante tudo isto recorro a lista e pergunto existe alguma forma mais simples de fazer o serviço de ler um .xls em linux e windows? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

grafico de efeitos espaciais - est. espacial
by Edson Ferreira 22 Nov '14

22 Nov '14

Paz e Bem !! estou com o seguinte codigo, para modelos CAR com dados de contagem (Poisson) e gostaria de plotar no grid criado um grafico em tom de cinzacom os efeitos espaciais e um outro com a média de cada regiãonão estou conseguindo.neste modelo qual a saída q mostra o valor estimado pro rho.?desde ja obrigado ############################ CARBayes Leroux vs CAR ####################################### ##################################################################### rm(list=ls()) require(spatstat) require(lattice) require(MASS) require(coda) require (MCMCpack) require(plotrix) require(CARBayes) ########## simulação x.easting= 1:10 x.northing <- 1:10 Grid <- expand.grid(x.easting, x.northing) n <- nrow(Grid) distance <-array(0, c(n,n)) W <-array(0, c(n,n)) for(i in 1:n) { for(j in 1:n) { temp <- (Grid[i,1] - Grid[j,1])^2 + (Grid[i,2] - Grid[j,2])^2 distance[i,j] <- sqrt(temp) if(temp==1) W[i,j] <- 1 } } x1= rnorm(n) # variavel 1 x2= rnorm(n) # variavel 2 E= rep(50,n) # valor esperado para cada regiao theta= rnorm(n, sd=0.05) # phi <- mvrnorm(n=1, mu=rep(0,n), Sigma=0.4 * exp(-0.1 * distance)) risk <- exp(-0.2 + 0.1 * x1 + 0.1*x2 + theta + phi) fitted <- E * risk Y <- rpois(n=n, lambda=fitted) formula <- Y ~ x1 + x2 + offset(log(E)) model <- poisson.lerouxCAR(formula=formula, W=W, burnin=50, n.sample=1000) model$formula model$samples model$fitted.values model$random.effects model$residuals model$W.summary model$DIC model$p.d model$summary.results model$model model$accept AttEdson

1 0

Diferença entre métodos
by Ari Clecius 22 Nov '14

22 Nov '14

Bom dia, estou rodando um experimento do tipo Box Behnken e usei a função reg e comparei com o obtido pela função rsm, o problema é que o que os dois consideram significativo difere, aonde posso estar errando? Testei também a função AIC para o modelo que ele fornece e dá diferente pela função AIC(). Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CMR library(DoE.wrapper)Design.1_2<- bbd.design(nfactors= 3 , ncenter= 4 ,randomize= FALSE,seed=7480, factor.names=list( A=c(0,400),B=c(600,2200),C=c(1,10) ) ) res=c(73.450,75.725,63.980,67.670,76.070,77.345,66.005,67.335,84.740,81.595,74.830,70.325,72.460,73.170,73.925,73.800)length(res)#Medidas resumoDesign.1_2<- add.response(Design.1_2,res, replace=FALSE)plancoded<- code.design(Design.1_2)plancoded LinearModel.1 <- aov(res ~(x1+x2+x3)^2+ I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)sreg = step(LinearModel.1,direction="both") #Melhor modelo#Step: AIC=14.07#res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2) # Df Sum of Sq RSS AIC#<none> 16.074 14.074#+ x1:x2 1 0.501 15.574 15.568#+ x2:x3 1 0.462 15.612 15.607#+ x1:x3 1 0.001 16.073 16.073#- x1 1 9.181 25.255 19.303#- I(x2^2) 1 9.310 25.384 19.385#- I(x3^2) 1 36.195 52.269 30.941#- x2 1 79.223 95.297 40.551#- I(x1^2) 1 86.793 102.867 41.774#- x3 1 212.747 228.821 54.566 summary(sreg) # Df Sum Sq Mean Sq F value Pr(>F) #x1 1 9.18 9.18 5.140 0.04959 * #x2 1 79.22 79.22 44.357 9.27e-05 ***#x3 1 212.75 212.75 119.118 1.72e-06 ***#I(x1^2) 1 86.79 86.79 48.596 6.53e-05 ***#I(x2^2) 1 9.31 9.31 5.213 0.04831 * #I(x3^2) 1 36.20 36.20 20.266 0.00149 ** #Residuals 9 16.07 1.79 #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Design.1_2.coded<-plancodedrsmModel.2 <- rsm(res ~ FO(x1, x2, x3) + TWI(x1, x2, x3) + PQ(x1, x2, x3), data=Design.1_2.coded)summary(rsmModel.2) # Estimate Std. Error t value Pr(>|t|) #(Intercept) 73.33875 0.79347 92.4275 1.081e-10 ***#x1 1.07125 0.56107 1.9093 0.104810 #x2 -3.14688 0.56107 -5.6087 0.001370 ** #x3 -5.15688 0.56107 -9.1911 9.350e-05 ***#x1:x2 0.35375 0.79347 0.4458 0.671360 #x1:x3 0.01375 0.79347 0.0173 0.986736 #x2:x3 -0.34000 0.79347 -0.4285 0.683246 #x1^2 -4.65813 0.79347 -5.8705 0.001081 ** #x2^2 1.52563 0.79347 1.9227 0.102878 #x3^2 3.00813 0.79347 3.7911 0.009062 ** #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Multiple R-squared: 0.9664, Adjusted R-squared: 0.916 #F-statistic: 19.17 on 9 and 6 DF, p-value: 0.0009315 LinearModel.2 <- aov(res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)summary(LinearModel.2 )AIC(LinearModel.2)#[1] 61.47997 #valor obtido anteriormente Step: AIC=14.07

1 0

Erro na instalação do R-2.15.0 no Ubuntu 14.04
by ASANTOS 21 Nov '14

21 Nov '14

Boa tarde Pessoal, Por motivos de pacotes que os autores não passaram para R3.0, no ubuntu 14.04 utilizo o R-2.15.0, no entanto, precisei formatar o computador e tentando instalar novamente o R-2.15.0, abaixado no com o arquivo R-2.15.0.tar.gz, não funciona, sendo que após a descompactação do arquivo, na pasta criada que fiz os passos: santosa@santosa-K43E:~$ cd Documentos santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install E recebi o erro no checking: configure: error: No F77 compiler found e quando tentei fazer o make install :santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install make: *** Sem regra para processar o alvo `install'. Pare. Alguém sabe se existe alguma incompatibilidade entre o R-2.15.0 e o ubuntu 14.04 ou estou fazendo algo errado? Obrigado, segue uma lista do resultado completo com vários no que recebi no processo de checking usando ./configure: santosa@santosa-K43E:~$ cd Documentos santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure checking build system type... i686-pc-linux-gnu checking host system type... i686-pc-linux-gnu loading site script './config.site' loading build-specific script './config.site' checking for pwd... /bin/pwd checking whether builddir is srcdir... yes checking for working aclocal... missing checking for working autoconf... missing checking for working automake... missing checking for working autoheader... missing checking for gawk... no checking for mawk... mawk checking whether ln -s works... yes checking for bison... no checking for byacc... no checking for ar... ar checking for a BSD-compatible install... /usr/bin/install -c checking for sed... /bin/sed checking for which... /usr/bin/which checking for less... /usr/bin/less checking for gtar... no checking for gnutar... no checking for tar... /bin/tar checking for tex... no checking for pdftex... no configure: WARNING: you cannot build PDF versions of the R manuals checking for pdflatex... no configure: WARNING: you cannot build PDF versions of vignettes and help pages checking for makeindex... no checking for makeinfo... no configure: WARNING: you cannot build info or HTML versions of the R manuals checking for texi2dvi... no checking for kpsewhich... no checking for unzip... /usr/bin/unzip checking for zip... /usr/bin/zip checking for gzip... /bin/gzip checking for bzip2... /bin/bzip2 checking for firefox... /usr/bin/firefox using default browser ... /usr/bin/firefox checking for acroread... no checking for acroread4... no checking for xdg-open... /usr/bin/xdg-open checking for notangle... false checking for pkg-config... /usr/bin/pkg-config checking for gcc... gcc checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -E checking for grep that handles long lines and -e... /bin/grep checking for egrep... /bin/grep -E checking whether gcc needs -traditional... no checking how to run the C preprocessor... gcc -E checking for gfortran... no checking for f95... no checking for fort... no checking for xlf95... no checking for ifort... no checking for ifc... no checking for efc... no checking for pgf95... no checking for lf95... no checking for gfortran... no checking for ftn... no checking for g95... no checking for f90... no checking for xlf90... no checking for pgf90... no checking for pghpf... no checking for epcf90... no checking for g77... no checking for f77... no checking for xlf... no checking for frt... no checking for pgf77... no checking for cf77... no checking for fort77... no checking for fl32... no checking for af77... no checking for fc... no configure: error: No F77 compiler found santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install make: *** Sem regra para processar o alvo `install'. Pare. -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

4 6

Re: [R-br] Erro na instalação do R-2.15.0 no Ubuntu 14.04
by Alisson Lucrécio 20 Nov '14

20 Nov '14

Voce baixou o arquivo .deb ou .taz? Segue umas sugestões que acho pertinente. É por isso, que a partir de agora salvo tmb minha library, assim mesmo que tenha o R-3.1.2, se tiver minha biblioteca salva posso carregar os pacotes dela. Vc já tenteu instalar com biocLite? source("http://bioconductor.org/biocLite.R") biocLite("biology") Att. On 20 November 2014 15:22, Alisson Lucrécio < alisson.lucrecio(a)ifgoiano.edu.br> wrote: > É por isso, que a partir de agora salvo tmb minha library, assim mesmo que > tenha o R-3.1.2, se tiver minha biblioteca salva posso carregar os pacotes > dela. > > Vc já tenteu instalar com biocLite? > > > On 20 November 2014 15:05, ASANTOS [via R-br] < > ml-node+s2285057n4663561h27(a)n4.nabble.com> wrote: > >> Boa tarde Pessoal, >> >> Por motivos de pacotes que os autores não passaram para >> R3.0, no ubuntu 14.04 utilizo o R-2.15.0, no entanto, precisei formatar >> o computador e tentando instalar novamente o R-2.15.0, abaixado no com o >> arquivo R-2.15.0.tar.gz, não funciona, sendo que após a descompactação >> do arquivo, na pasta criada que fiz os passos: >> >> santosa@santosa-K43E:~$ cd Documentos >> santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> >> >> E recebi o erro no checking: configure: error: No F77 >> compiler found e quando tentei fazer o make install >> :santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> make: *** Sem regra para processar o alvo `install'. Pare. >> >> Alguém sabe se existe alguma incompatibilidade entre >> o R-2.15.0 e o ubuntu 14.04 ou estou fazendo algo errado? >> Obrigado, segue uma lista do resultado completo com vários no que recebi >> no processo de checking usando ./configure: >> >> santosa@santosa-K43E:~$ cd Documentos >> santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure >> checking build system type... i686-pc-linux-gnu >> checking host system type... i686-pc-linux-gnu >> loading site script './config.site' >> loading build-specific script './config.site' >> checking for pwd... /bin/pwd >> checking whether builddir is srcdir... yes >> checking for working aclocal... missing >> checking for working autoconf... missing >> checking for working automake... missing >> checking for working autoheader... missing >> checking for gawk... no >> checking for mawk... mawk >> checking whether ln -s works... yes >> checking for bison... no >> checking for byacc... no >> checking for ar... ar >> checking for a BSD-compatible install... /usr/bin/install -c >> checking for sed... /bin/sed >> checking for which... /usr/bin/which >> checking for less... /usr/bin/less >> checking for gtar... no >> checking for gnutar... no >> checking for tar... /bin/tar >> checking for tex... no >> checking for pdftex... no >> configure: WARNING: you cannot build PDF versions of the R manuals >> checking for pdflatex... no >> configure: WARNING: you cannot build PDF versions of vignettes and help >> pages >> checking for makeindex... no >> checking for makeinfo... no >> configure: WARNING: you cannot build info or HTML versions of the R >> manuals >> checking for texi2dvi... no >> checking for kpsewhich... no >> checking for unzip... /usr/bin/unzip >> checking for zip... /usr/bin/zip >> checking for gzip... /bin/gzip >> checking for bzip2... /bin/bzip2 >> checking for firefox... /usr/bin/firefox >> using default browser ... /usr/bin/firefox >> checking for acroread... no >> checking for acroread4... no >> checking for xdg-open... /usr/bin/xdg-open >> checking for notangle... false >> checking for pkg-config... /usr/bin/pkg-config >> checking for gcc... gcc >> checking whether the C compiler works... yes >> checking for C compiler default output file name... a.out >> checking for suffix of executables... >> checking whether we are cross compiling... no >> checking for suffix of object files... o >> checking whether we are using the GNU C compiler... yes >> checking whether gcc accepts -g... yes >> checking for gcc option to accept ISO C89... none needed >> checking how to run the C preprocessor... gcc -E >> checking for grep that handles long lines and -e... /bin/grep >> checking for egrep... /bin/grep -E >> checking whether gcc needs -traditional... no >> checking how to run the C preprocessor... gcc -E >> checking for gfortran... no >> checking for f95... no >> checking for fort... no >> checking for xlf95... no >> checking for ifort... no >> checking for ifc... no >> checking for efc... no >> checking for pgf95... no >> checking for lf95... no >> checking for gfortran... no >> checking for ftn... no >> checking for g95... no >> checking for f90... no >> checking for xlf90... no >> checking for pgf90... no >> checking for pghpf... no >> checking for epcf90... no >> checking for g77... no >> checking for f77... no >> checking for xlf... no >> checking for frt... no >> checking for pgf77... no >> checking for cf77... no >> checking for fort77... no >> checking for fl32... no >> checking for af77... no >> checking for fc... no >> configure: error: No F77 compiler found >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> make: *** Sem regra para processar o alvo `install'. Pare. >> >> -- >> ====================================================================== >> Alexandre dos Santos >> Proteção Florestal >> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso >> Campus Cáceres >> Caixa Postal 244 >> Avenida dos Ramires, s/n >> Bairro: Distrito Industrial >> Cáceres - MT CEP: 78.200-000 >> Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) >> e-mails:[hidden email] >> <http://user/SendEmail.jtp?type=node&node=4663561&i=0> >> [hidden email] >> <http://user/SendEmail.jtp?type=node&node=4663561&i=1> >> Lattes: http://lattes.cnpq.br/1360403201088680 >> ====================================================================== >> >> >> --- >> Este email foi escaneado pelo Avast antivírus. >> http://www.avast.com >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4663561&i=2> >> 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. >> >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Erro-na-instalacao-do-R-2-15-0-no-Ub… >> To unsubscribe from R-br, click here >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

1 1

Dificuldade na utilização de um loop duplo para a criação de várias amostras
by Thiago Touya 20 Nov '14

20 Nov '14

Boa tarde! Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'. Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'. Segue o código que estou usando para a geração de uma amostra: ## Tamanho da Amostra n=100 ## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5) ## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) } ## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] } ## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) } O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra. Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

4 5