- R-br - DINF UFPR

subset de observações específicas de uma variável
by cayanportela＠hotmail.com 08 Jul '16

08 Jul '16

is.factor(dados1$HScode) [1] TRUE dados1$HScode <- as.character(dados1$HScode) dados2 <- subset(dados1,nchar(dados1$HScode)==4) Caso a sua variavel esteja definida como fator, transforme para character e use o nchar para contagem do numero de caracteres. <quote author='R-br mailing list'> Boa noite! Tenho uma base de dados em que uma das variáveis (HScode) possui três níveis. Por exemplo, 01 0101 010122 0102 010233 02 0201 020179 Contudo preciso de um subset apenas com as observações com quatro dígitos. Ou seja, 0101 0102 0201 structure(list(reporter = c("Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania"), HScode = c("01", "0101", "010110", "010190", "0102", "010210", "010290", "0103", "010310", "010391", "010392", "0104", "010410", "010420", "0105", "010511", "010512", "010519", "010594", "010599", "0106", "010611", "010612", "010619", "010620", "010631", "010632", "010639", "010690", "02", "0201", "020110", "020120", "020130", "0202"), boundtariff = c("6.571969696969697", "8.75", "7.5", "10", "7.291666666666667", "5", "9.583333333333334", "5", "5", "5", "5", "5", "5", "5", "9.5", "5", "5", "5", "12.5", "20", "5", "5", "5", "5", "5", "5", "5", "5", "5", "10.27914614121511", "10", "10", "10", "10", "10")), datalabel = "", time.stamp = " 7 Jul 2016 22:53", .Names = c("reporter", "HScode", "boundtariff"), formats = c("%54s", "%9s", "%20s"), types = c(54L, 7L, 20L), val.labels = c("", "", ""), var.labels = c("", "", ""), version = 12L, row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35"), class = "data.frame") Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE Site: www.ccsa.ufpb.br/ppge orcid www.orcid.org/0000-0002-1288-0135 _______________________________________________ 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. iD icon.gif (1K) <http://r-br.2285057.n4.nabble.com/attachment/4666373/0/iD%20icon.gif> </quote> Quoted from: http://r-br.2285057.n4.nabble.com/R-br-subset-de-observacoes-especificas-de… _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

subset de observações específicas de uma variável
by Alexandre Loures 08 Jul '16

08 Jul '16

Boa noite! Tenho uma base de dados em que uma das variáveis (HScode) possui três níveis. Por exemplo, 01 0101 010122 0102 010233 02 0201 020179 Contudo preciso de um subset apenas com as observações com quatro dígitos. Ou seja, 0101 0102 0201 structure(list(reporter = c("Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania"), HScode = c("01", "0101", "010110", "010190", "0102", "010210", "010290", "0103", "010310", "010391", "010392", "0104", "010410", "010420", "0105", "010511", "010512", "010519", "010594", "010599", "0106", "010611", "010612", "010619", "010620", "010631", "010632", "010639", "010690", "02", "0201", "020110", "020120", "020130", "0202"), boundtariff = c("6.571969696969697", "8.75", "7.5", "10", "7.291666666666667", "5", "9.583333333333334", "5", "5", "5", "5", "5", "5", "5", "9.5", "5", "5", "5", "12.5", "20", "5", "5", "5", "5", "5", "5", "5", "5", "5", "10.27914614121511", "10", "10", "10", "10", "10")), datalabel = "", time.stamp = " 7 Jul 2016 22:53", .Names = c("reporter", "HScode", "boundtariff"), formats = c("%54s", "%9s", "%20s"), types = c(54L, 7L, 20L), val.labels = c("", "", ""), var.labels = c("", "", ""), version = 12L, row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35"), class = "data.frame") Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE Site: www.ccsa.ufpb.br/ppge orcid www.orcid.org/0000-0002-1288-0135

3 2

Extrair dados de ajuste
by Pedro Campelo 05 Jul '16

05 Jul '16

Pessoal, Estou rodando um script e preciso extrair os dados dos coeficientes ajustados para um modelo. No meu exemplo, quero extrair e separar os valores "k" e "n". *> data.frame(A)* * M R* *1 10 24.61* *2 20 29.41* *3 30 33.07* *4 40 37.37* *5 50 38.22* *6 60 41.14* *7 70 41.61* *8 80 43.41* *9 90 47.18* *10 100 48.77* *11 110 51.02* *12 120 52.95* *13 130 52.98* *14 140 52.92* *> model<-nls(R~k*M^n, start=list(k=0.5,n=1))* *> summary(model)* Obrigado *----* *Pedro Henrique Campelo Felix* *Bachelor of Food Engineering (UNI-BH)* *Msc. Chemical Engineering (UFMG)* *PhD Student of Food Science (UFLA)* *+55 (31) 97560-1545 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*

3 3

Frequência de um registro específico
by Wagner Tassinari 01 Jul '16

01 Jul '16

Olá pessoal, Existe algum comando para verificarmos a frequencia de um registro em específico, por exemplo: Meu banco está dessa forma: var1 var2 .... vark 1 2 0 4 5 6 0 6 0 4 ... Gostaria de saber a frequencia de zeros da var1 obrigadoi ----- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Cel: (21) 96488-5982 (WhatsApp) Skype: wagner.tassinari wtassinari(a)gmail.com ----------------------------------------------------------------------

5 4

Inputando registros
by Wagner Tassinari 30 Jun '16

30 Jun '16

Olá pessoal, tenho o seguinte banco de dados: var1 var2 20 25 40 40 45 40 50 55 50 60 60 Quero criar uma var3 que seja igual a var2 mas onde está em branco impute os valores de var1, ex: var3 25 40 45 40 55 50 60 Estou fazendo da seguinte forma: banco$var1=as.character(banco$var1) banco$var2=as.character(banco$var2) banco$var3= banco$var2 banco$var3[banco$var2 == " "] = banco$var1 Mas não esta dando certo !!!!!!!!!! Obrigado pessoal ----- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Cel: (21) 96488-5982 (WhatsApp) Skype: wagner.tassinari wtassinari(a)gmail.com ----------------------------------------------------------------------

4 3

Dúvidas - Scott-Knott
by Pedro Campelo 28 Jun '16

28 Jun '16

Pessoal, Preciso de uma ajuda na execução desse script abaixo: *#Data.Frame* *> data.frame(A)* * M R* *1 1 0.150* *2 1 0.117* *3 1 0.125* *4 2 0.312* *5 2 0.288* *6 2 0.262* *7 3 0.318* *8 3 0.315* *9 3 0.299* *10 4 0.395* *11 4 0.312* *12 4 0.365* *#ANOVA One Way* *av<-aov(R ~ M, data.frame(A))* *summary(av)* *#Teste de Scott Knott* *sk <- SK(x=av, y=A$R, model='y~x', which='x', dispersion='se')* *summary(sk)* Quando vou executar a rotina, aparece esta mensagem: *> sk <- SK(x=av, y=A$R, model='y~x', which='x', dispersion='se')* *Error in SK.aov(x = av, y = A$R, model = "y~x", which = "x", dispersion = "se") : * * No factors in the fitted model!* *In addition: Warning message:* *In replications(paste("~", xx), data = mf) : non-factors ignored: M* Como poderia resolver isto? Obrigado *----* *Pedro Henrique Campelo Felix* *Bachelor of Food Engineering (UNI-BH)* *Msc. Chemical Engineering (UFMG)* *PhD Student of Food Science (UFLA)* *+55 (31) 97560-1545 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*

3 5

obter data.frames com mesmas dimensões
by Heloíse Pavanato 28 Jun '16

28 Jun '16

Boa tarde, Eu tenho dois data.frames com número de linhas diferentes (nrow(rep1C) = 73, nrow(rep2C) = 45), e colunas iguais (19). Eu quero que ambos data.frames fiquem com a mesma dimensão. Para isso, rep2C deve ser aumentado com "NA". Porém, quero manter os mesmos índices "area", "season", "month" para ambos data.frames. Alguém tem alguma idéia de que função usar? Segue CMR: #---------------------------------------------------------------------- ## Load data # file: *https://www.datafilehost.com/d/64a6a8d9 <https://www.datafilehost.com/d/64a6a8d9>* load('dataC.inia.rda') ## Building matrix for replicates # Sampling 1 rep1C <- dataC.inia[dataC.inia$replicate == 1,] table(rep1C$area, rep1C$season, rep1C$transect) dim(rep1C) # Sampling 2 rep2 <- dataC.inia[dataC$replicate == 2,] table(rep2$area, rep2$season, rep2$transect) dim(rep2) #---------------------------------------------------------------------- Obrigada pela ajuda! Heloise

3 6

Inversa condicional
by Andre Oliveira 27 Jun '16

27 Jun '16

Bom dia, além da função ginv do pacote do MASS para calcula inversa condicional. Qual outro pacote faz esta mesma função só que outras inversas sem ser a de apenas a inversa condicional de Moore-Penrose? obrigado. 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

2 1

Normalização de dados
by Elias Carvalho 24 Jun '16

24 Jun '16

Pessoal, Eu preciso fazer uma análise para verificar a relação das variáveis com a morte e meus dados são semelhantes a tabela abaixo: *CD_MUN* *MUNIC* *UF* *POPULAÇÃO* *UTI* *PIB* *AUTO* *MOTO* *ONIBUS* *CAMINHÃO* *MORTE* 1 A X 3577 10 18533,94 1549 661 40 360 2 *2* *B* *Y* *11958* *2* *12487,36* *5402* *2234* *185* *1401* *11* 3 C Z 6936 8 9947 2187 1478 60 666 6 Nessa tabela o município da linha 2 tem mais população, mais PIB, mais veículos e mais morte. Isso é óbvio porque a sua população e maior, mas isso não necessariamente quer dizer que esse município tenha mais mortes ou que elas são causadas mais provavelmente pela frota de carros ou motos. Então pensando em evitar viéses, e por isso eu pensei em normalizar os dados dividindo as colunas UTI, PIB, AUTO, MOTO, ONIBUS, CAMINHÃO e MORTE pela coluna da POPULAÇÃO e dividir por 10.000, por exemplo para MORTE da linha 2, ficaria: (11/11958)*1000 = 9,1989. Sessa forma a tabela ficaria assim: *CD_MUN* *MUNIC* *UF* *POPULAÇÃO* *UTI* *PIB* *AUTO* *MOTO* *ONIBUS* *CAMINHÃO* *MORTE* 1 A X 3577 27,9564 51814,2018 4330,4445 1847,9172 111,8256 1006,4300 5,5913 *2* *B* *Y* *11958* *1,6725* *10442,6827* *4517,4778* *1868,2054* *154,7081* *1171,6006* *9,1989* 3 C Z 6936 11,5340 14341,1188 3153,1142 2130,9112 86,5052 960,2076 8,6505 Ou ainda poderia usar alguma ferramento do R que fosse mais confiável ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

2 2

Operação entre datas
by Amikobh 24 Jun '16

24 Jun '16

Prezados, Desde já, obrigado pela ajuda. Tendo duas variáveis com datas referentes a determinada pessoa: NASCIMENTO DATA_FATO 26/11/97 23/02/15 20/03/99 28/03/15 14/05/99 08/04/15 06/05/98 08/04/15 03/01/00 08/04/15 19/05/97 08/05/15 01/02/01 17/05/15 28/11/97 03/06/15 10/02/00 03/06/15 Como procedo para criar mais uma coluna que calcula a idade no dia da ocorrência de determinado fato? De tal forma que obtenha: NASCIMENTO DATA_FATO IDADE 26/11/97 23/02/15 17 20/03/99 28/03/15 16 14/05/99 08/04/15 15 06/05/98 08/04/15 16 03/01/00 08/04/15 15 19/05/97 08/05/15 17 01/02/01 17/05/15 14 28/11/97 03/06/15 17 10/02/00 03/06/15 15 Valeu, gente!

4 7