{"action":"opened","issue":{"url":"https://api.github.com/repos/IgnasiLucas/mosquito/issues/8","repository_url":"https://api.github.com/repos/IgnasiLucas/mosquito","labels_url":"https://api.github.com/repos/IgnasiLucas/mosquito/issues/8/labels{/name}","comments_url":"https://api.github.com/repos/IgnasiLucas/mosquito/issues/8/comments","events_url":"https://api.github.com/repos/IgnasiLucas/mosquito/issues/8/events","html_url":"https://github.com/IgnasiLucas/mosquito/issues/8","id":168759534,"number":8,"title":"Increase the number of reads in the sample fastq in 2016-07-26","user":{"login":"IgnasiLucas","id":16137858,"avatar_url":"https://avatars.githubusercontent.com/u/16137858?v=3","gravatar_id":"","url":"https://api.github.com/users/IgnasiLucas","html_url":"https://github.com/IgnasiLucas","followers_url":"https://api.github.com/users/IgnasiLucas/followers","following_url":"https://api.github.com/users/IgnasiLucas/following{/other_user}","gists_url":"https://api.github.com/users/IgnasiLucas/gists{/gist_id}","starred_url":"https://api.github.com/users/IgnasiLucas/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/IgnasiLucas/subscriptions","organizations_url":"https://api.github.com/users/IgnasiLucas/orgs","repos_url":"https://api.github.com/users/IgnasiLucas/repos","events_url":"https://api.github.com/users/IgnasiLucas/events{/privacy}","received_events_url":"https://api.github.com/users/IgnasiLucas/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2016-08-01T22:44:38Z","updated_at":"2016-08-01T22:44:38Z","closed_at":null,"body":"Currently, there are only 10 sequences sampled from each individual, which is the default in sample_seqs. That's why the same number of clusters are generated under all values of wclust. We need at least 100000 sequences per sample. Get more (500000) if that does not produce a meaningful graph. Then plot the total number of clusters (Y axis) against the wclust value (X axis). To do this, I would write an R script (plot.R) and execute it from the README.sh. The plot.R script could be something like this:\r\n\r\n# This reads the table in a data frame:\r\nclusters <- read.table('summary_clust.txt', header=TRUE)\r\n# This makes the variables in the data frame directly available in the work space:\r\nattach(clusters)\r\n# This opens the png file where to send the plot:\r\npng(filename='molestus_clusters.png')\r\n# This plots the first line:\r\nplot(Wclust[muestra == 'Mol01_se'], total[muestra == 'Mol01_se'], type='l')\r\n# This adds the second line. Do the same for subsequent ones:\r\nlines(Wclust[muestra == 'Mol02_se'], total[muestra == 'Mol02_se'])\r\n...\r\n# Finally, it is necessary to close the png file:\r\ndev.off()\r\n\r\nCheck help(plot) in R for details. Please, separate the pipiens plot in a different png file. @mgmorant "}}
{"actor":{"display_login":"IgnasiLucas"}}