Maio 2016 - R-br - DINF UFPR

ajuda com implementação de Programação Linear com o R - constrOptim
by Cleber N.Borges 26 Mai '16

26 Mai '16

Prezados, Bom dia! Como na lista tem pessoas com conhecimento e experiência, peço algumas orientações... :-) Imagino que o método de "Programação Linear (LP)" deve servir para resolver meu problema de achar os parâmetros de uma dada função objetivo. Sei que no R tem a opção nativa "constrOptim" (poderia ser outra, posso instalar pacotes) Entretanto eu não sei como elaborar minhas condições de contorno em formato matricial que é o modo como o comando constrOptim trabalha (até onde consegui entender) Agradeceria muito se alguém puder me ajudar, respondendo se eu realmente estou certo em pensar que a LP é a ferramenta certa para meu caso e /(ou que) me indique um material básico (introdutório). Abaixo segue um exemplo simplificado do tipos de função que pretendo obter parâmetros. Agradeço antecipadamente pela ajuda e por qualquer dica Cleber Borges ##################### # fobjetivo <- function( x = 0.1 ){ k1 <- 0.007585776 k2 <- 6.16595e-08 k3 <- 2.137962e-13 k4 <- 1e-14 k <- k1 * k2 * k3 * k4 # Determinar as 6 variáveis: a, b, c, d, e, f ################### # 6 Condições de contorno: ( Igualdades ) ################### k1 == a * d / c # 1 k2 == a * e / d # 2 k3 == a * f / e # 3 k4 == a * b # 4 4 * a == b + d + e + f # 5 x == c + d + e + f # 6 return( c( a, b, c, d, e, f ) ) } --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

3 4

Substituindo caracteres em palavras
by Wagner Tassinari 25 Mai '16

25 Mai '16

Olá pessoal, Tenho a seguinte amostra com os nomes: [1] "PIETRO LUCAS DA SILVA NAO MORA" "ROSECLEIA DE JESUS PORTUGAL NAO MORA MA" [3] "TAINA OLIVEIRA NAO MORA" "DANIELLE DE ALMEIDA MELO DA SILVA NAO E" [5] "DAVI FERREIRA DE CARVALHO NAO RESIDE NE" "EDNALDO DA SILVA NAO RESIDE MAIS NA ARE" [7] "MARCIA FERREIRA DE MORAIS NAO RESIDE NE" "MARIA JOANA LUDGERIO DE PAULO NAO MORA" [9] "ROSANGELA DE BRITO NAO R" "NATANAEL DOS SANTOS MACENA NAO MORA MAIS" Esses nomes apresentam imperfeições no preenchimento, daí eu gostaria de excluir tudo que está após a palavra NAO e criar outra variável para os nomes ficarem limpos dessa forma: [1] "PIETRO LUCAS DA SILVA" "ROSECLEIA DE JESUS PORTUGAL" [3] "TAINA OLIVEIRA" "DANIELLE DE ALMEIDA MELO DA SILVA" [5] "DAVI FERREIRA DE CARVALHO" "EDNALDO DA SILVA" [7] "MARCIA FERREIRA DE MORAIS" "MARIA JOANA LUDGERIO DE PAULO" [9] "ROSANGELA DE BRITO " "NATANAEL DOS SANTOS MACENA" 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 ----------------------------------------------------------------------

2 1

como fazer plot de componentes principais em cores?
by Jackson Rodrigues 24 Mai '16

24 Mai '16

Oi pessoal, Eu preciso fazer um PCA que varie as cores em função dos anos. O mais antigo em vermelho e a medida que chega as datas mais recentes vai ficando roxo (uma escala de cor). eu tenho esse exemplo hipotético. mat5 <- matrix(rnorm(2000), ,4) # minha matriz mat5 ages<-seq(1, 50000, by = 100) #idades ages library(plotrix) agecol<-color.scale(ages,extremes=c("purple","red")) # encaixando a escala de cor nas idades rownames(mat5) <- c(seq ) # fazendo as idades os nomes da linhas da matriz library("vegan") test.pca<-rda(mat5) # calculando componentes principais plot(test.pca,display="sites", cex=1, type="p",scaling=-3) # plotando componentes princiapais. O problema é: Como fazer os pontos do gráfico ficarem coloridos como estabelecido em " agecol<-color.scale(ages,extremes=c("purple","red")) " Agradeço desde já a ajuda Abraços Jackson

3 2

Subistituindo Caracteres em palavras
by Wagner Tassinari 24 Mai '16

24 Mai '16

Olá pessoal, Tenho a seguinte amostra com os nomes: [1] "PIETRO LUCAS DA SILVA NAO MORA" "ROSECLEIA DE JESUS PORTUGAL NAO MORA MA" [3] "TAINA OLIVEIRA NAO MORA" "DANIELLE DE ALMEIDA MELO DA SILVA NAO E" [5] "DAVI FERREIRA DE CARVALHO NAO RESIDE NE" "EDNALDO DA SILVA NAO RESIDE MAIS NA ARE" [7] "MARCIA FERREIRA DE MORAIS NAO RESIDE NE" "MARIA JOANA LUDGERIO DE PAULO NAO MORA" [9] "ROSANGELA DE BRITO NAO R" "NATANAEL DOS SANTOS MACENA NAO MORA MAIS" Esses nomes apresentam imperfeições no preenchimento, daí eu gostaria de excluir tudo que está após a palavra NAO e criar outra variável para os nomes ficarem limpos dessa forma: [1] "PIETRO LUCAS DA SILVA" "ROSECLEIA DE JESUS PORTUGAL" [3] "TAINA OLIVEIRA" "DANIELLE DE ALMEIDA MELO DA SILVA" [5] "DAVI FERREIRA DE CARVALHO" "EDNALDO DA SILVA" [7] "MARCIA FERREIRA DE MORAIS" "MARIA JOANA LUDGERIO DE PAULO" [9] "ROSANGELA DE BRITO " "NATANAEL DOS SANTOS MACENA" 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 ----------------------------------------------------------------------

3 2

invalid multibyte string
by Wagner Tassinari 24 Mai '16

24 Mai '16

Estou tentando verificar o número de caracteres de uma variável em doi bancos de dados. Um banco eu consigo e no outro me dá o erro: invalid multibyte string, element 593 Alguém saberia como posso lidar com isso ? Obrigado e abs --- Segue o comando reproduzido abaixo: > table(nchar(bancolink3$nom_pes)) # Observando o numero de caracteres do nome da pessoa 3 4 5 6 7 8 9 10 11 12 13 14 15 16 2 1 1 4 6 21 209 843 2156 4112 6802 9330 11849 15355 17 18 19 20 21 22 23 24 25 26 27 28 29 30 20044 32021 49598 76109 105911 129330 143014 144552 131724 117983 103914 93237 85936 78879 31 32 33 34 35 36 37 38 39 40 41 42 43 44 70489 60862 49874 37813 28071 20123 14360 10916 7915 5641 3430 2594 1868 1225 45 46 47 48 49 50 51 52 53 54 55 56 57 58 860 583 392 242 185 118 63 51 20 18 14 6 6 6 59 60 61 62 63 69 70 1 5 1 1 1 1 2 > table(nchar(fichaAlink3$NOME)) # Observando o numero de caracteres do nome da pessoa Error in nchar(fichaAlink3$NOME) : invalid multibyte string, element 593 ----- 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 ----------------------------------------------------------------------

2 1

Re: [R-br] out of topic
by sznelwar＠uol.com.br 23 Mai '16

23 Mai '16

2 1

ANALISe DE EXPERIMENTO DE ADUBAÇÃO COM COLHEITAS ANUAIS
by Wenceslau Teixeira 23 Mai '16

23 Mai '16

Caros um colega me pediu auxilio para analisar um experimento cujas colheitas foram repetidas por nove anos. O experimento tinha cinco tratamentos (não eram doses) de fertilizantes, que foram aplicados anualmente nas mesmas parcelas, em um delinamento em blocos. Era feito plantio e colheita sucessivas (em todos os anos) de milho e de soja (por nove anos). Os parâmetros medidos foram colheitas anuais (kg ha-1), teores foliares, analise de solos. Alguem tem alguma literatura ou exemplo (script) sobre experimentos com repetição no tempo ? Grato, Wenceslau Teixeira Embrapa Solos -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

Enc: Re: Eixo x
by suzabio 22 Mai '16

22 Mai '16

Oi Éder, Mais cedo, eu mandei esse email com um pedaço do banco de dados em txt, mas não foi. Eu estou reenviando por que o comando labels parece suficiente para mim, mas não funcionou. Eu coloquei os comandos que vc mandou antes do comando do axis 1 e dentro dele coloquei labels = T, mas não deu certo. O outro comando deu certo, mas o gráfico ficou estranho. No comando tá dizendo que os dados têm distribuição normal? Pq eles não têm. Achei que poderia ser por isso que o gráfico ficou estranho. Obrigada e um abraço,Suzana Enviado do meu dispositivo Samsung -------- Mensagem original -------- De : suzana de souza rodrigues <suzabio(a)yahoo.com.br> Data: 21/05/2016 09:24 (GMT-03:00) Para: r-br(a)listas.c3sl.ufpr.br, comunello.eder(a)gmail.com Assunto: Re: [R-br] Eixo x Oi Éder,Vou mandar em anexo um pedaço do banco em txt. Em Sábado, 21 de Maio de 2016 8:29, Éder Comunello <comunello.eder(a)gmail.com> escreveu: Suzana, bom dia! É o tipo de situação em que você deveria fornecer uma amostra dos seus dados ou pelo menos informar a natureza deles. A julgar pelo gráfico, seu objeto "Semana" é um vetor de inteiros. Nesse caso, uma opção seria criar um novo objeto do tipo "Date", onde cada observação poderia ser a data central da semana em questão. Nesse caso, acredito que baste utilizar esse novo objeto no lugar de "Semana" com a opção axes=T, sem precisar recorrer ao comando axis(). Mas é só um palpite. ================================================ Éder ComunelloAgronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa)Dourados, MS, Brazil |<O>|================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 Em 20 de maio de 2016 21:22, suzabio <suzabio(a)yahoo.com.br> escreveu: Olá,Eu criei um gráfico com esses comandos no R studio:plot (Semana, RESPMe5, xlab = "Year", ylab = "Admissions", main = "Respiratory diseases under 5 years", axes = F) lines ((lowess(Semana, RESPMe5, f = 0.05)),lwd = 2) axis(1, at = seq(0,370, by = 52), pos = 0) axis(2, at = seq(0, 1100, by = 100), pos = 0) Ao invés de aparecer o número de semanas no eixo x, gostaria de escrever o nome dessas semanas: jan-2009, jan-2010, e assim por diante, como eu faço isso?Obrigada!Suzana Enviado do meu dispositivo Samsung _______________________________________________ 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 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� c�igo m�imo reproduz�el.

2 1

Criar duas novas variáveis a partir de outra
by Alexandre Loures 21 Mai '16

21 Mai '16

Boa tarde! Tenho uma variável como abaixo e quero criar duas novas variáveis com a seguinte condição: Primeira variável a ser criada (exp): para todo F_3idoXidd_2_ * a nova variável assumirá BRA e para todo F_3idoXidd_3_* a nova variável assumirá ARG e assim por diante. Segundo variável a ser criada (imp): para todo F_3idoXidd_*_2 a nova variável assumirá BRA e para todo F_3idoXidd_*_3 a nova variável assumirá ARG e assim por diante. P.S.: O asterisco (*) indica qualquer outro número, ou seja, para todo F_3idoXidd_2_* sempre a nova variável será BRA. F_3idoXidd_2_2 F_3idoXidd_2_3 F_3idoXidd_2_4 F_3idoXidd_2_5 . . . F_3idoXidd_2_27 F_3idoXidd_2_28 F_3idoXidd_2_29 F_3idoXidd_2_30 F_3idoXidd_3_2 F_3idoXidd_3_3 F_3idoXidd_3_4 F_3idoXidd_3_5 . . . F_3idoXidd_3_27 F_3idoXidd_3_28 F_3idoXidd_3_29 F_3idoXidd_3_30 . . . Alguém poderia me ajudar? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

1 0

Re: [R-br] Eixo x
by suzabio 21 Mai '16

21 Mai '16

Bom dia Éder, Realmente Semana é um vetor de números inteiros. O gráfico é número de semanas por número de internações por problemas respiratórios. Vc acha que eu devo criar esse novo objeto com o comando as.Date? É a melhor solução? Um abraço e obrigada! Enviado do meu dispositivo Samsung -------- Mensagem original -------- De : Éder Comunello <comunello.eder(a)gmail.com> Data: 21/05/2016 08:28 (GMT-03:00) Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Eixo x Suzana, bom dia! É o tipo de situação em que você deveria fornecer uma amostra dos seus dados ou pelo menos informar a natureza deles. A julgar pelo gráfico, seu objeto "Semana" é um vetor de inteiros. Nesse caso, uma opção seria criar um novo objeto do tipo "Date", onde cada observação poderia ser a data central da semana em questão. Nesse caso, acredito que baste utilizar esse novo objeto no lugar de "Semana" com a opção axes=T, sem precisar recorrer ao comando axis(). Mas é só um palpite. ================================================ Éder ComunelloAgronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa)Dourados, MS, Brazil |<O>|================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 Em 20 de maio de 2016 21:22, suzabio <suzabio(a)yahoo.com.br> escreveu: Olá,Eu criei um gráfico com esses comandos no R studio:plot (Semana, RESPMe5, xlab = "Year", ylab = "Admissions", main = "Respiratory diseases under 5 years", axes = F) lines ((lowess(Semana, RESPMe5, f = 0.05)),lwd = 2) axis(1, at = seq(0,370, by = 52), pos = 0) axis(2, at = seq(0, 1100, by = 100), pos = 0) Ao invés de aparecer o número de semanas no eixo x, gostaria de escrever o nome dessas semanas: jan-2009, jan-2010, e assim por diante, como eu faço isso?Obrigada!Suzana Enviado do meu dispositivo Samsung _______________________________________________ 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.

2 2