- R-br - DINF UFPR

Criar função para tirar a média entre grupos sob condições dadas
by ASANTOS 31 Jan '16

31 Jan '16

Caros listeiros, Estou tentando criar uma função aggPestFE que para uma primeira variável se for "Formiga", vai tirar a media da coluna 5, mas se for "Besouros" ou "Barata" vai tira a média da coluna 4. Sendo que quero o resultado dado em forma de list, ficando o meu resultado: [[1]] [1] "Barata" [[1]][[3]] Projeto Talhao Injuria(%) 1 Abre Campo 15 25 2 Volta Redonda 73 0 3 Vitoria 74 50 [[2]] [1] "Besouros" [[1]][[3]] Projeto Talhao Injuria(%) 1 Abre Campo 15 0 2 Volta Redonda 73 25 3 Vitoria 74 50 [[2]] [1] "Formiga" [[1]][[3]] Projeto Talhao Area (m2) 1 Abre Campo 15 12.5 2 Volta Redonda 73 1.5 3 Vitoria 74 23.8 Para isso estou tentando: ##Dados artificiais Bichos<-c("Besouros","Barata","Formiga","Barata","Formiga","Besouros","Formiga", "Besouros","Barata","Formiga","Barata","Formiga","Besouros","Formiga","Formiga") Talhao<-c(73,15,74,75,15,15,15,73,15,15,73,15,73,74,74) Projeto<-c("Abre Campo", "Volta Redonda","Vitoria","Abre Campo", "Volta Redonda", "Vitoria","Abre Campo", "Volta Redonda","Vitoria","Abre Campo", "Volta Redonda", "Vitoria","Abre Campo", "Volta Redonda","Vitoria") Injuria<-c(25,100,0,25,0,100,0,50,25,0,25,0,25,0,0) Area<-c(0,0,12.5,0,7.5,0,1.5,0,0,23.8,0,5.3,0,2.0,11.3) tableFE<-data.frame(Bichos, Talhao, Projeto, Injuria, Area) ## Banco de dados criado ##Função aggPestFE<-function(x, db=tableFE){ mylist <- list() if (tableFE[,1] != "Formigas") { for (i in length(tableFE[,1])){ GP_FE2<-tableFE[tableFE[,1]==tableFE[i],] aggdata <-aggregate(GP_FE2[,4], list(GP_FE2[,3],GP_FE2[,2]), mean) colnames(aggdata)<-c("Regional","Projeto","Talhão","Injúria média (%)") tmp <- list(aggdata) mylist[NPRAGS[i]] <- tmp result<-mylist return(result) } else if (tableFE[,1] == "Formigas") { GP_FE2<-[tableFE[,1] aggdata <-aggregate(GP_FE2[,5], list(GP_FE2[,3],GP_FE2[,2]), mean) colnames(aggdata)<-c("Regional","Projeto","Talhão","Injúria média (%)") tmp <- list(aggdata) mylist[NPRAGS[i]] <- tmp result<-mylist return(result) } } RES<-list(result,aggdata3) result<-RES return(result) } #Teste aggPestFE(tableFE) # Sem sucesso, alguém poderia me ajudar? Obrigado e abraços, -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

5 23

Separar strings no R
by Wagner Tassinari 30 Jan '16

30 Jan '16

Saudações grupo, Tenho uma variável que representa o nome completo de uma pessoa, ex: Manuel dos Santos de Jesus. Gostaria de separar em duas variáveis o nome e o úlltimo sobrenome respectivamente, por exemplo: VarNome = Manuel VarSobrenome = Jesus. Agradeço muitíssimo a todos, -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

6 6

Erro função lme
by Victor Fioreze 29 Jan '16

29 Jan '16

Olá retorno aqui com a minha duvida anterior: Estou rodando um script com a função lme do pacote nlme e esta dando um erro.Após eu rodar o modelo usando lme aparece: Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Revisei os dados e o script e esta tudo ok.O que pode estar acontecendo? Att, pois não sei responder no fórum.... O arquivoa ser carregado esta no anexo Segue o meu comando para que alguém possa tentar> cruza <- read.table("C:/Users/Victor/Dropbox/projeto bezerras/Dados/cernelha.txt", header=T,)> cruza Bezerro Tratamento ACd Data DN DJN Sexo ACD Idade1 423 1 81.33 41937 41930 41948 FALSE 81.0 672 423 1 81.59 41944 41930 41948 FALSE 81.0 743 423 1 82.07 41951 41930 41948 FALSE 81.0 814 423 1 83.62 41958 41930 41948 FALSE 81.0 885 423 1 83.72 41965 41930 41948 FALSE 81.0 906 425 1 81.38 41937 41930 41962 FALSE 81.0 677 425 1 81.95 41944 41930 41962 FALSE 81.0 748 425 1 82.40 41951 41930 41962 FALSE 81.0 819 425 1 83.00 41958 41930 41962 FALSE 81.0 8810 425 1 83.07 41965 41930 41962 FALSE 81.0 9011 426 1 81.09 41937 41930 41968 FALSE 81.0 6712 426 1 81.69 41944 41930 41968 FALSE 81.0 7413 426 1 82.29 41951 41930 41968 FALSE 81.0 8114 426 1 83.20 41958 41930 41968 FALSE 81.0 8815 426 1 83.28 41965 41930 41968 FALSE 81.0 9016 428 1 81.69 41937 41930 41980 FALSE 81.0 6717 428 1 82.32 41944 41930 41980 FALSE 81.0 7418 428 1 83.34 41951 41930 41980 FALSE 81.0 8119 428 1 83.65 41958 41930 41980 FALSE 81.0 8820 428 1 83.74 41965 41930 41980 FALSE 81.0 9021 421 2 82.72 41937 41930 41931 FALSE 79.5 6722 421 2 84.07 41944 41930 41931 FALSE 79.5 7423 421 2 84.98 41951 41930 41931 FALSE 79.5 8124 421 2 85.77 41958 41930 41931 FALSE 79.5 8825 421 2 85.99 41965 41930 41931 FALSE 79.5 9026 429 2 81.29 41937 41930 41989 FALSE 81.0 6727 429 2 83.21 41944 41930 41989 FALSE 81.0 7428 429 2 84.47 41951 41930 41989 FALSE 81.0 8129 429 2 85.75 41958 41930 41989 FALSE 81.0 8830 429 2 85.91 41965 41930 41989 FALSE 81.0 9031 430 2 85.61 41937 41930 42007 FALSE 83.5 6732 430 2 87.20 41944 41930 42007 FALSE 83.5 7433 430 2 87.53 41951 41930 42007 FALSE 83.5 8134 430 2 87.87 41958 41930 42007 FALSE 83.5 8835 430 2 88.02 41965 41930 42007 FALSE 83.5 9036 438 2 90.66 41937 41930 42063 FALSE 89.0 6737 438 2 91.99 41944 41930 42063 FALSE 89.0 7438 438 2 92.80 41951 41930 42063 FALSE 89.0 8139 438 2 94.07 41958 41930 42063 FALSE 89.0 8840 438 2 94.25 41965 41930 42063 FALSE 89.0 9041 422 3 83.43 41937 41930 41933 FALSE 83.0 6742 422 3 84.17 41944 41930 41933 FALSE 83.0 7443 422 3 84.65 41951 41930 41933 FALSE 83.0 8144 422 3 84.91 41958 41930 41933 FALSE 83.0 8845 422 3 84.97 41965 41930 41933 FALSE 83.0 9046 427 3 79.06 41937 41930 41978 FALSE 78.4 6747 427 3 79.86 41944 41930 41978 FALSE 78.4 7448 427 3 80.65 41951 41930 41978 FALSE 78.4 8149 427 3 81.66 41958 41930 41978 FALSE 78.4 8850 427 3 81.78 41965 41930 41978 FALSE 78.4 9051 432 3 82.79 41937 41930 42018 FALSE 80.2 6752 432 3 83.80 41944 41930 42018 FALSE 80.2 7453 432 3 84.64 41951 41930 42018 FALSE 80.2 8154 432 3 85.60 41958 41930 42018 FALSE 80.2 8855 432 3 85.79 41965 41930 42018 FALSE 80.2 9056 435 3 89.95 41937 41930 42025 FALSE 88.5 6757 435 3 90.69 41944 41930 42025 FALSE 88.5 7458 435 3 91.34 41951 41930 42025 FALSE 88.5 8159 435 3 92.30 41958 41930 42025 FALSE 88.5 8860 435 3 92.43 41965 41930 42025 FALSE 88.5 9061 424 4 82.57 41937 41930 41954 FALSE 82.5 6762 424 4 83.02 41944 41930 41954 FALSE 82.5 7463 424 4 83.17 41951 41930 41954 FALSE 82.5 8164 424 4 83.33 41958 41930 41954 FALSE 82.5 8865 424 4 83.35 41965 41930 41954 FALSE 82.5 9066 431 4 82.99 41937 41930 42013 FALSE 81.0 6767 431 4 85.85 41944 41930 42013 FALSE 81.0 7468 431 4 87.54 41951 41930 42013 FALSE 81.0 8169 431 4 89.33 41958 41930 42013 FALSE 81.0 8870 431 4 89.62 41965 41930 42013 FALSE 81.0 9071 433 4 84.94 41937 41930 42020 FALSE 81.4 6772 433 4 86.87 41944 41930 42020 FALSE 81.4 7473 433 4 88.33 41951 41930 42020 FALSE 81.4 8174 433 4 88.99 41958 41930 42020 FALSE 81.4 8875 433 4 89.25 41965 41930 42020 FALSE 81.4 9076 436 4 87.64 41937 41930 42036 FALSE 87.3 6777 436 4 88.81 41944 41930 42036 FALSE 87.3 7478 436 4 89.17 41951 41930 42036 FALSE 87.3 8179 436 4 90.48 41958 41930 42036 FALSE 87.3 8880 436 4 90.58 41965 41930 42036 FALSE 87.3 90 str(cruza)'data.frame': 80 obs. of 9 variables: $ Bezerro : int 423 423 423 423 423 425 425 425 425 425 ... $ Tratamento: int 1 1 1 1 1 1 1 1 1 1 ... $ ACd : num 81.3 81.6 82.1 83.6 83.7 ... $ Data : int 41937 41944 41951 41958 41965 41937 41944 41951 41958 41965 ... $ DN : int 41930 41930 41930 41930 41930 41930 41930 41930 41930 41930 ... $ DJN : int 41948 41948 41948 41948 41948 41962 41962 41962 41962 41962 ... $ Sexo : logi FALSE FALSE FALSE FALSE FALSE FALSE ... $ ACD : num 81 81 81 81 81 81 81 81 81 81 ... $ Idade : int 67 74 81 88 90 67 74 81 88 90 ... library(nlme)> attach(cruza)> fTrat<-as.factor(Tratamento)> class(fTrat)[1] "factor"> class(ACd)[1] "numeric"> class(ACD)[1] "numeric"> class(Data)[1] "integer"> Data<-as.factor(Data)> simcomp.ran <- lme(ACd ~ DJN + ACD + fTrat + Data + fTrat:Data + ACD*fTrat + fTrat*DJN,+ random= ~1 | Bezerro/fTrat)Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Muito obrigado pela ajuda e desculpe não saber responder no fórum. Abraço a todos! Victor Ionatan FiorezeMédico Veterinário e ZootecnistaMestrando em Nutrição de Ruminantes, PPGZ-UFPelBolsista Capes Integrante GERUMEN (53) 81032887lattes: http://lattes.cnpq.br/5026587958713372

2 1

Aumentar a memória do R
by Wagner Tassinari 28 Jan '16

28 Jan '16

Olá pessoal, Como faço para liberar mais memoria para o R rodar ? Estou utilizando os comandos: > memory.size() [1] 3306.77 > memory.size(TRUE) [1] 3995 > memory.limit() [1] 32676 Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

3 2

Re: [R-br] Dúvida em looping e data.frame
by Paulo Henrique Pimenta 28 Jan '16

28 Jan '16

Usei a função simulacao <- function(proporcao, moedas){ for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) a <- print(moedas(n,p)) return(a) }b <- simulacao(proporcao, moedas)No final ele possui somente um vetor com os dados gerados com n = 30 e p = 0.5 e não todos os gerados pelas combinações. Como eu poderia ter uma vetor ou alguma outra forma de armazenar estes dados? Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435. > Message: 11 > Date: Wed, 27 Jan 2016 09:33:10 +0300 > From: Paulo Henrique Pimenta <paulopimenta6(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Dúvida em looping e data.frame > Message-ID: <BAY182-W241C2B8A79DFC8287003AFF4D90(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos > Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > ####Laboratorio 1 de R > ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL > moeda <- function(p){ > if(runif(1) < p){ > y <- 1 > }else{ > y <- 0 } return(y) } ####separando uma funcao da outra#### > moedas <- function(numeros,probabilidade){ > for (i in 1:numeros){ > vetor[i] <- moeda(probabilidade) > } > return(vetor) > } > > ####separando uma funcao da outra#### > proporcao <- function(moedas){ > numeros_1 <- sum(vetor_moedas) > numeros_totais <- length(vetor_moedas) > proporcao_final <- numeros_1/numeros_totais > return(proporcao_final) > } > ####separando uma funcao da outra#### > simulacao <- function(proporcao, moedas){ > ###variando o "n" e o "p" > for (i in 1:500){ > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... > n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > } > > return(a) > } > ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas > ###comeco das chamas > #########################################prob <- 0.4 > moeda(prob) > #########################################n_1 <- 100 > moedas(n_1, prob) > ######################################### > vetor_moedas <- moedas(n_1, prob) > ######################################### > proporcao(vetor_moedas) > ######################################### > simulacao(proporcao,moedas) > #################################### > Um grande abraço a todos!! > > Paulo Henrique de A. S. Pimenta. > Graduando em Meteorologia (Bacharelado) - IAG/USP. > Fone: +5511981318435. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/683e3d06/atta…> > > ------------------------------ > > Message: 12 > Date: Wed, 27 Jan 2016 10:57:57 -0300 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida em looping e data.frame > Message-ID: > <CABmC8gmxHtjGenrn6us=VnsgRE2E1B_zCY_nBZ+Vnyp1CZX5nw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Paulo, bom dia! > > A instrução for (i in 1:500) não parece coerente, pois me parece que você > deveria variar em função da iteração n x p. Ao variar n você irá criar > vetores de tamanhos diferentes, impossibilitando o armazenamento direto em > um data.frame(). Poderia utilizar uma lista. > > Você pode tentar algo como: > simulacao <- function(proporcao, moedas){for (n in c(10, 20, 30)) for (p in > c(0.3, 0.4, 0.5)) print(moedas(n,p))} > simulacao(proporcao, moedas) > > Outro ponto é que você poderia usar rbinom() para os sorteios, enxugando > seu código. Deixo uma ideia abaixo: > > ### <code r> > args(rbinom) ### function (n, size, prob) > vetor_moedas <- rbinom(100, 1, .4); vetor_moedas > sum(vetor_moedas)/length(vetor_moedas) ### proporção final > > for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) print(rbinom(n, 1,p)) > # [1] 0 1 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 0 0 1 > # [1] 0 1 1 0 0 0 0 1 0 1 > # [1] 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 > # [1] 0 1 0 0 0 1 0 1 0 0 0 1 1 0 0 1 0 0 0 0 > # [1] 0 0 0 1 0 1 1 1 0 1 1 0 1 0 0 1 1 1 1 1 > # [1] 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 1 1 1 0 1 0 1 0 0 1 1 1 0 0 0 > # [1] 1 0 0 1 1 1 1 1 1 1 0 0 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 1 0 > ### </code> > > > > ================================================ > Éder Comunello > PhD Student in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > > > > Em 27 de janeiro de 2016 03:33, Paulo Henrique Pimenta < > paulopimenta6(a)hotmail.com> escreveu: > > > Olá a todos > > > > Estou escrevendo um código que calcula o lance de uma moeda jogada e a > > probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, > > num bloco do meu código fiquei com uma dúvida relacionada a vetorização. > > Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > > > ####Laboratorio 1 de R > > > > ###cara eh 1 > > ###coroa eh 0 > > vetor_moedas <- NULL > > vetor <- NULL > > a <- NULL > > > > moeda <- function(p){ > > > > if(runif(1) < p){ > > > > y <- 1 > > > > }else{ > > > > y <- 0 > > > > } > > return(y) > > > > } > > > > ####separando uma funcao da outra#### > > > > moedas <- function(numeros,probabilidade){ > > > > for (i in 1:numeros){ > > > > vetor[i] <- > > moeda(probabilidade) > > > > } > > > > return(vetor) > > > > } > > > > > > ####separando uma funcao da outra#### > > > > proporcao <- function(moedas){ > > > > numeros_1 <- sum(vetor_moedas) > > > > numeros_totais <- length(vetor_moedas) > > > > proporcao_final <- numeros_1/numeros_totais > > > > return(proporcao_final) > > > > } > > > > ####separando uma funcao da outra#### > > > > simulacao <- function(proporcao, moedas){ > > > > ###variando o "n" e o "p" > > > > for (i in 1:500){ > > > > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como > > data.frame todos os dados feitos pela função moedas, como fazer com que a2 > > também aloque um outro data.frame que aloque informações da função moedas e > > assim por diante... > > > > n <- c(10, 20, 30, 40, 50) > > p <- c(0.2, 0.3, 0.4, 0.5) > > > > a[i] <- data.frame(moedas(n[i],p[i])) > > > > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e > > recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > > > > } > > > > > > return(a) > > > > } > > > > ############################## > > ############################## > > ###Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda > > ###lancamento de 10 moedas - funcao moedas > > > > ###comeco das chamas > > > > ######################################### > > prob <- 0.4 > > > > moeda(prob) > > > > ######################################### > > n_1 <- 100 > > > > moedas(n_1, prob) > > > > ######################################### > > > > vetor_moedas <- moedas(n_1, prob) > > > > ######################################### > > > > proporcao(vetor_moedas) > > > > ######################################### > > > > simulacao(proporcao,moedas) > > > > ### > > ### > > ############################## > > > > > > Um grande abraço a todos!! > > > > Paulo Henrique de A. S. Pimenta. > > > > Graduando em Meteorologia (Bacharelado) - IAG/USP. > > Fone: +5511981318435.

2 2

Re: [R-br] Digest R-br, volume 61, assunto 17
by Paulo Henrique Pimenta 27 Jan '16

27 Jan '16

Usei a função simulacao <- function(proporcao, moedas){ for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) a <- print(moedas(n,p)) return(a) } b <- simulacao(proporcao, moedas) No final ele possui somente um vetor com os dados gerados com n = 30 e p = 0.5 e não todos os gerados pelas combinações. Como eu poderia ter uma vetor ou alguma outra forma de armazenar estes dados? Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435. > > Message: 11 > Date: Wed, 27 Jan 2016 09:33:10 +0300 > From: Paulo Henrique Pimenta <paulopimenta6(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Dúvida em looping e data.frame > Message-ID: <BAY182-W241C2B8A79DFC8287003AFF4D90(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos > Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > ####Laboratorio 1 de R > ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL > moeda <- function(p){ > if(runif(1) < p){ > y <- 1 > }else{ > y <- 0 } return(y) } ####separando uma funcao da outra#### > moedas <- function(numeros,probabilidade){ > for (i in 1:numeros){ > vetor[i] <- moeda(probabilidade) > } > return(vetor) > } > > ####separando uma funcao da outra#### > proporcao <- function(moedas){ > numeros_1 <- sum(vetor_moedas) > numeros_totais <- length(vetor_moedas) > proporcao_final <- numeros_1/numeros_totais > return(proporcao_final) > } > ####separando uma funcao da outra#### > simulacao <- function(proporcao, moedas){ > ###variando o "n" e o "p" > for (i in 1:500){ > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... > n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > } > > return(a) > } > ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas > ###comeco das chamas > #########################################prob <- 0.4 > moeda(prob) > #########################################n_1 <- 100 > moedas(n_1, prob) > ######################################### > vetor_moedas <- moedas(n_1, prob) > ######################################### > proporcao(vetor_moedas) > ######################################### > simulacao(proporcao,moedas) > #################################### > Um grande abraço a todos!! > > Paulo Henrique de A. S. Pimenta. > Graduando em Meteorologia (Bacharelado) - IAG/USP. > Fone: +5511981318435. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/683e3d06/atta…> > > ------------------------------ > > Message: 12 > Date: Wed, 27 Jan 2016 10:57:57 -0300 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida em looping e data.frame > Message-ID: > <CABmC8gmxHtjGenrn6us=VnsgRE2E1B_zCY_nBZ+Vnyp1CZX5nw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Paulo, bom dia! > > A instrução for (i in 1:500) não parece coerente, pois me parece que você > deveria variar em função da iteração n x p. Ao variar n você irá criar > vetores de tamanhos diferentes, impossibilitando o armazenamento direto em > um data.frame(). Poderia utilizar uma lista. > > Você pode tentar algo como: > simulacao <- function(proporcao, moedas){for (n in c(10, 20, 30)) for (p in > c(0.3, 0.4, 0.5)) print(moedas(n,p))} > simulacao(proporcao, moedas) > > Outro ponto é que você poderia usar rbinom() para os sorteios, enxugando > seu código. Deixo uma ideia abaixo: > > ### <code r> > args(rbinom) ### function (n, size, prob) > vetor_moedas <- rbinom(100, 1, .4); vetor_moedas > sum(vetor_moedas)/length(vetor_moedas) ### proporção final > > for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) print(rbinom(n, 1,p)) > # [1] 0 1 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 0 0 1 > # [1] 0 1 1 0 0 0 0 1 0 1 > # [1] 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 > # [1] 0 1 0 0 0 1 0 1 0 0 0 1 1 0 0 1 0 0 0 0 > # [1] 0 0 0 1 0 1 1 1 0 1 1 0 1 0 0 1 1 1 1 1 > # [1] 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 1 1 1 0 1 0 1 0 0 1 1 1 0 0 0 > # [1] 1 0 0 1 1 1 1 1 1 1 0 0 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 1 0 > ### </code> > > > > ================================================ > Éder Comunello > PhD Student in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > > > > Em 27 de janeiro de 2016 03:33, Paulo Henrique Pimenta < > paulopimenta6(a)hotmail.com> escreveu: > > > Olá a todos > > > > Estou escrevendo um código que calcula o lance de uma moeda jogada e a > > probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, > > num bloco do meu código fiquei com uma dúvida relacionada a vetorização. > > Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > > > ####Laboratorio 1 de R > > > > ###cara eh 1 > > ###coroa eh 0 > > vetor_moedas <- NULL > > vetor <- NULL > > a <- NULL > > > > moeda <- function(p){ > > > > if(runif(1) < p){ > > > > y <- 1 > > > > }else{ > > > > y <- 0 > > > > } > > return(y) > > > > } > > > > ####separando uma funcao da outra#### > > > > moedas <- function(numeros,probabilidade){ > > > > for (i in 1:numeros){ > > > > vetor[i] <- > > moeda(probabilidade) > > > > } > > > > return(vetor) > > > > } > > > > > > ####separando uma funcao da outra#### > > > > proporcao <- function(moedas){ > > > > numeros_1 <- sum(vetor_moedas) > > > > numeros_totais <- length(vetor_moedas) > > > > proporcao_final <- numeros_1/numeros_totais > > > > return(proporcao_final) > > > > } > > > > ####separando uma funcao da outra#### > > > > simulacao <- function(proporcao, moedas){ > > > > ###variando o "n" e o "p" > > > > for (i in 1:500){ > > > > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como > > data.frame todos os dados feitos pela função moedas, como fazer com que a2 > > também aloque um outro data.frame que aloque informações da função moedas e > > assim por diante... > > > > n <- c(10, 20, 30, 40, 50) > > p <- c(0.2, 0.3, 0.4, 0.5) > > > > a[i] <- data.frame(moedas(n[i],p[i])) > > > > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e > > recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > > > > } > > > > > > return(a) > > > > } > > > > ############################## > > ############################## > > ###Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda > > ###lancamento de 10 moedas - funcao moedas > > > > ###comeco das chamas > > > > ######################################### > > prob <- 0.4 > > > > moeda(prob) > > > > ######################################### > > n_1 <- 100 > > > > moedas(n_1, prob) > > > > ######################################### > > > > vetor_moedas <- moedas(n_1, prob) > > > > ######################################### > > > > proporcao(vetor_moedas) > > > > ######################################### > > > > simulacao(proporcao,moedas) > > > > ### > > ### > > ############################## > > > > > > Um grande abraço a todos!! > > > > Paulo Henrique de A. S. Pimenta. > > > > Graduando em Meteorologia (Bacharelado) - IAG/USP. > > Fone: +5511981318435. > > > > _______________________________________________ > > 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/20160127/515a392f/atta…> > > ------------------------------ > > 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 61, assunto 17 > *****************************************

1 1

Validando o CPF
by Wagner Tassinari 27 Jan '16

27 Jan '16

Olá pessoal, Alguém conhece algum código em R para verificar a consistência do número do CPF ? Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

2 1

Dúvida em looping e data.frame
by Paulo Henrique Pimenta 27 Jan '16

27 Jan '16

Olá a todos Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: ####Laboratorio 1 de R ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL moeda <- function(p){ if(runif(1) < p){ y <- 1 }else{ y <- 0 } return(y) } ####separando uma funcao da outra#### moedas <- function(numeros,probabilidade){ for (i in 1:numeros){ vetor[i] <- moeda(probabilidade) } return(vetor) } ####separando uma funcao da outra#### proporcao <- function(moedas){ numeros_1 <- sum(vetor_moedas) numeros_totais <- length(vetor_moedas) proporcao_final <- numeros_1/numeros_totais return(proporcao_final) } ####separando uma funcao da outra#### simulacao <- function(proporcao, moedas){ ###variando o "n" e o "p" for (i in 1:500){ \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! } return(a) } ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas ###comeco das chamas #########################################prob <- 0.4 moeda(prob) #########################################n_1 <- 100 moedas(n_1, prob) ######################################### vetor_moedas <- moedas(n_1, prob) ######################################### proporcao(vetor_moedas) ######################################### simulacao(proporcao,moedas) #################################### Um grande abraço a todos!! Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

2 1

Criar um convex hull para um objeto shapefile
by ASANTOS 26 Jan '16

26 Jan '16

Prezados Listeiros, Eu gostaria de criar um convex hull que envolvesse 4 polígonos, representados por um objeto oriundo da leitura de um readShapeLines(), sem ter que transformar o objeto em SpatialPolygonsDataFrame, isso é possível? Meu CRM: #Pacotes require(rgdal) require(maptools) #------------------------------------------------------------------------------- #Crio 4 polígonos sr <- SpatialPolygons(list( Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409, 180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676, 332618, 332413, 332349)))),'1'), Polygons(list(Polygon(cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437, 179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683, 331133, 331623, 332152, 332357, 332373)))),'2'), Polygons(list(Polygon(cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875, 179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110), c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004, 329783, 329665, 329720, 329933, 330478, 331062, 331086)))),'3'), Polygons(list(Polygon(cbind(c(180304, 180403,179632,179420,180304), c(332791, 333204, 333635, 333058, 332791)))),'4'))) plot(sr) #Converto em polígono espacial srdf=SpatialPolygonsDataFrame(sr, data.frame(row.names=c('1','2','3','4'), PIDS=1:4)) srdf@data #Criação do shapefile para ilustrar o problema writeOGR(srdf, getwd(), 'POLY', 'ESRI Shapefile') #Leitura do shapefile contorno_line_X <- readShapeLines ("POLY.shp") #Plot plot(contorno_line_X) #Tentativa de criar o convex hull df.data = as.data.frame(contorno_line_X) ch <- chull(df.data) lines(ch, col="red") ##Não funciona Obrigado, -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 5

Pacote rJava
by Gislaine Batistela 26 Jan '16

26 Jan '16

Olá! Estou com dificuldade para carregar o pacote rJava. Por favor, alguém poderia me ajudar? Vejam as mensagens: > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- trying URL 'http://www.vps.fmvz.usp.br/CRAN/bin/windows/contrib/3.0/rJava_0.9-6.zip' Content type 'application/zip' length 758820 bytes (741 Kb) opened URL downloaded 741 Kb package ‘rJava’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\Unesp\AppData\Local\Temp\Rtmpsx6PLo\downloaded_packages > local({pkg <- select.list(sort(.packages(all.available = TRUE)),graphics=TRUE) + if(nchar(pkg)) library(pkg, character.only=TRUE)}) Error : .onLoad failed in loadNamespace() for 'rJava', details: call: fun(libname, pkgname) error: JAVA_HOME cannot be determined from the Registry Error: package or namespace load failed for ‘rJava’ -- Gislaine Cristina Batistela

3 3