Como consta no link da pergunta, um exemplo mínimo reproduzível deve ter os seguintes conteúdos:
- Um pequeno conjunto de dados;
- O menor código possível que seja executável e que reproduza o erro no pequeno conjunto de dados mencionado;
- As informações sobre a versão do
R
e do sistema em que está rodando o código, bem como dos pacotes utilizados;
- Se for utilizar dados aleatórios, garantir que os resultados sejam os mesmos;
Nessa resposta vou elencar algumas das principais funções em R
para cumprir essas tarefas.
Vale lembrar que os exemplos das páginas de ajuda das funções do R
podem ser de grande valia para ter uma noção da estrutura de um exemplo mínimo reproduzível. Em geral, os códigos dos exemplos da ajuda do R
satisfazem a esses requisitos.
Produzindo o conjunto de dados
Para usar seu próprio conjuto de dados, a função dput()
, juntamente com head()
pode ser bastante útil. Por exemplo o código abaixo fornece as 10 primeiras observações da base de dados iris
já com a estrutura necessária para "remontar" a base de dados. Assim, para quem for tentar responder a sua pergunta, basta copiar e colar o código em structure()
.
dput(head(iris, 10))
#> structure(list(Sepal.Length = c(5.1, 4.9, 4.7, 4.6, 5, 5.4, 4.6,
#> 5, 4.4, 4.9), Sepal.Width = c(3.5, 3, 3.2, 3.1, 3.6, 3.9, 3.4,
#> 3.4, 2.9, 3.1), Petal.Length = c(1.4, 1.4, 1.3, 1.5, 1.4, 1.7,
#> 1.4, 1.5, 1.4, 1.5), Petal.Width = c(0.2, 0.2, 0.2, 0.2, 0.2,
#> 0.4, 0.3, 0.2, 0.2, 0.1), Species = structure(c(1L, 1L, 1L, 1L,
#> 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("setosa", "versicolor", "virginica"
#> ), class = "factor")), .Names = c("Sepal.Length", "Sepal.Width",
#> "Petal.Length", "Petal.Width", "Species"), row.names = c(NA,
#> 10L), class = "data.frame")
Reproduzindo os dados:
dados <- structure(list(Sepal.Length = c(
5.1, 4.9, 4.7, 4.6, 5, 5.4, 4.6,
5, 4.4, 4.9
), Sepal.Width = c(
3.5, 3, 3.2, 3.1, 3.6, 3.9, 3.4,
3.4, 2.9, 3.1
), Petal.Length = c(
1.4, 1.4, 1.3, 1.5, 1.4, 1.7,
1.4, 1.5, 1.4, 1.5
), Petal.Width = c(
0.2, 0.2, 0.2, 0.2, 0.2,
0.4, 0.3, 0.2, 0.2, 0.1
), Species = structure(c(
1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L
), .Label = c("setosa", "versicolor", "virginica"), class = "factor")), .Names = c(
"Sepal.Length", "Sepal.Width",
"Petal.Length", "Petal.Width", "Species"
), row.names = c(
NA,
10L
), class = "data.frame")
Uma solução menos ideal do que esta seria fornecer os dados em formato texto, como, por exemplo, no caso abaixo:
texto <- "Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1 5.1 3.5 1.4 0.2 setosa
2 4.9 3.0 1.4 0.2 setosa
3 4.7 3.2 1.3 0.2 setosa
4 4.6 3.1 1.5 0.2 setosa
5 5.0 3.6 1.4 0.2 setosa
6 5.4 3.9 1.7 0.4 setosa"
Neste caso, o usuário que for responder sua pergunta pode remontar a base de dados utilizando a função read.table()
:
dados <- read.table(text=texto)
Outra forma de produzir um conjunto de dados é gerando valores aleatórios, por exemplo, com a função rnorm()
(você também pode gerar de outras distribuições sem ser a normal, caso seja pertinente) ou com a função sample()
para uma amostragem de valores de algum vetor. Um caso útil pode ser a função letters()
, para gerar caracteres ou fatores. Neste caso, não esqueça de fornecer a seed
para que o exemplo seja reproduzível.
Exemplo:
set.seed(1) # garantir reproducibilidade
dados <- data.frame(x = rnorm(10), y = sample(letters, 10))
dados
#> x y
#> 1 -0.6264538 y
#> 2 0.1836433 f
#> 3 -0.8356286 p
#> 4 1.5952808 c
#> 5 0.3295078 z
#> 6 -0.8204684 i
#> 7 0.4874291 a
#> 8 0.7383247 h
#> 9 0.5757814 x
#> 10 -0.3053884 v
Outras funções interessantes neste caso são as funções do tipo as
, como as.factor()
, as.Date()
etc, para você converter os dados para o formato necessário.
Produzindo o código mínimo
Tente identificar a menor parte necessária do seu código que gere o erro ou a dúvida que você tem. Antes de enviar o código, certifique-se de que você elencou os pacotes necessários para que ele seja reproduzível. Para isso, é bom testar o seu código após reiniciar a sessão do R
, para ter certeza de que tudo o que é necessário está lá.
Exemplo:
library(lattice) # a biblioteca utilizada
set.seed(1) # a seed
dados <- data.frame(x = as.character(rnorm(10)), y = sample(letters, 10)) # o conjunto de dados
densityplot(as.numeric(dados$x))
as.numeric(dados$x)
#> [1] 2 5 4 10 6 3 7 9 8 1
Este exemplo corresponderia a uma pergunta do tipo: "Estou tentando fazer um gráfico de densidade com o lattice
como no código acima, porque quando converto os dados para numeric eles viram 2, 5, 4 ... e não permanecem como os dados originais do rnorm
?"
Informações do sistema
Por fim, quando necessário, você pode fornecer as informações do seu sistema com sessionInfo()
, que dá informações detalhadas da sua seção. No meu caso, essas informações foram:
R version 3.0.1 (2013-05-16)
Platform: i386-w64-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252
[3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C
[5] LC_TIME=Portuguese_Brazil.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] lattice_0.20-15
loaded via a namespace (and not attached):
[1] grid_3.0.1 tools_3.0.1
Pacote reprex
Para ajudar a criar o exemplo reproduzível o pacote reprex pode ser bastante útil, inclusive os exemplos anteriores foram gerado nele. Este é um pacote feito especificamente para ajudar a criar e executar exemplos reproduzíveis (o nome reprex é a abreviação de Reproducible Example), já com formatação para sites como GitHub e StackOverflow.
Uma forma simples de criar um exemplo reproduzível com o pacote é copiar o código em R
para sua área de transferência. Depois basta carregar o pacote com library(reprex)
e rodar o comando reprex(venue = "so")
que o código com os resultados comentados já formatados estarão disponíveis para ser colados na venue escolhida (nesse exemplo "so" é a venue stackoverflow). Todas as imagens geradas são colocadas no imgur e o link é gerado automaticamente para postagem, bastando colar o resultado.
O pacote tem outras funções bastante úteis. Por exemplo, você pode automaticamente incluir as informações do sistema com o argumento si = TRUE
e também automaticamente formatar o seu código usando o estilo sugerido pelo Hadley com o argumento style = TRUE
. Para mais informações veja a página do pacote.