tag:blogger.com,1999:blog-4946490806848569840.comments2018-10-22T07:28:29.626-07:00R snippetsBogumił Kamińskinoreply@blogger.comBlogger213125tag:blogger.com,1999:blog-4946490806848569840.post-19905726505208403132018-10-22T07:28:29.626-07:002018-10-22T07:28:29.626-07:00For larger data probably R is not the best tool fo...For larger data probably R is not the best tool for the job.Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-43555592932145209742018-10-22T07:22:55.818-07:002018-10-22T07:22:55.818-07:00Can you help me with visualizing this model on zip...Can you help me with visualizing this model on zip code level data? thanksUnknownhttps://www.blogger.com/profile/00170846508259360447noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-87470814008191401422018-06-29T00:36:57.454-07:002018-06-29T00:36:57.454-07:00Thanks for sharing about stack overflow survey. w...Thanks for sharing about stack overflow survey. we are clearly understand this concept.we do have dedicated and experience developer.<a href="https://www.developeronrent.com/hire-python-developers.html" rel="nofollow">Hire Python Developer</a>developeronrent dorhttps://www.blogger.com/profile/01234746697832066061noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-75123850753464539102018-06-19T07:01:01.002-07:002018-06-19T07:01:01.002-07:00Spot on with this write-up, I actually suppose thi...Spot on with this write-up, I actually suppose this web site wants far more consideration. I’ll probably be once more to learn far more, thanks for that info. <a href="https://online-casinos.us.org" rel="nofollow">usa online casino</a>Fillip Kirkorovhttps://www.blogger.com/profile/04617626994297367122noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-81002622042197615502018-06-16T08:09:52.489-07:002018-06-16T08:09:52.489-07:00Once I initially commented I clicked the -Notify m...Once I initially commented I clicked the -Notify me when new comments are added- checkbox and now every time a comment is added I get 4 emails with the identical comment. Is there any method you can remove me from that service? Thanks! <a href="https://online-casinos.us.org" rel="nofollow">casino games</a>Yose Abramovichttps://www.blogger.com/profile/11140186761749280113noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-75738759057591525122018-04-18T13:05:40.205-07:002018-04-18T13:05:40.205-07:00The answer is yes to both questionsThe answer is yes to both questionsBogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-17735139959923283802018-04-18T07:14:41.211-07:002018-04-18T07:14:41.211-07:00HI there: Two questions: 1) do the permutations ne...HI there: Two questions: 1) do the permutations need to match the number of constituencies in the variable constit? 2) Is the number 107 the number required to win the election? Simonhttps://www.blogger.com/profile/05981296362335774054noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-89503172257600969522017-06-28T01:59:04.178-07:002017-06-28T01:59:04.178-07:00This is exactly how I understand it :) and that is...This is exactly how I understand it :) and that is why I show what programming languages Julia users use.Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-87578045528702245932017-06-28T01:28:03.782-07:002017-06-28T01:28:03.782-07:00Here we have a great example of Simpson's para...Here we have a great example of Simpson's paradox ;).<br />R is used almost 100% by scientific oriented programmers (statistical research, data science...) while python is only used by a fraction as a scientific language, as many people uses as a general purpose language (web programming, OS scripting, apps..). I am almost positive that comparing the people who uses both languages as a scientific oriented language you will obtain similar results (numpy, pandas...).Gabrielhttps://www.blogger.com/profile/09238993347653608368noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-77979245212621100552017-06-27T22:08:48.505-07:002017-06-27T22:08:48.505-07:00I wanted to say that among all languages R is the ...I wanted to say that among all languages R is the language who's users user Julia most. I agree that the wording is confusing and will correct it.Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-86405451118110973742017-06-27T16:55:45.530-07:002017-06-27T16:55:45.530-07:00I'm confused. It says "most" R users...I'm confused. It says "most" R users also use Julia then it says 2.26%. 2% of R users are "most" R users?Unknownhttps://www.blogger.com/profile/10525668526999547235noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-68266099081588945922017-01-16T04:43:49.113-08:002017-01-16T04:43:49.113-08:00great articlegreat articleAbhijit Chandahttps://www.blogger.com/profile/09188944739603626888noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-31156634314219241572016-12-18T22:07:06.510-08:002016-12-18T22:07:06.510-08:00Thank you very much for your help! I will try to u...Thank you very much for your help! I will try to understand this function.Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-13885661492251314632016-12-18T10:04:06.686-08:002016-12-18T10:04:06.686-08:00My function: unfortunately :(, I think that the be...My function: unfortunately :(, I think that the best advice is to read documentation of all functions I use in R help.<br /><br />Generators: https://en.wikipedia.org/wiki/Generator_(computer_programming)Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-41872857698365228842016-12-18T04:47:53.912-08:002016-12-18T04:47:53.912-08:00Thank you very much! I had attempted to generate a...Thank you very much! I had attempted to generate all subsets of {1,...,200} to test an assumption of EBIC, but I just found it would take at least 2^200 bytes to store its power sets. What should I read to understand your all.subsets.fast() function? What should I read to understand "generator"? I am a grad student in stats, have a bsc in math. Thank you very much!Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-9285473693288582672016-12-08T00:39:17.578-08:002016-12-08T00:39:17.578-08:00Such a large structure simply does not fit into me...Such a large structure simply does not fit into memory.<br />You would have to use generator rather than materialized structure. Anyway 2^50=1125899906842624 so you will not be able to iterate over such a large number of elements anyway.Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-7839262773817293912016-12-07T19:02:08.590-08:002016-12-07T19:02:08.590-08:00Hi, your code is much faster than the set_power() ...Hi, your code is much faster than the set_power() in sets package, but when I try all.subsets.fast(seq(1:50)), it says:<br />Error: cannot allocate vector of size 4194304.0 Gb. How can I fix it? Thanks!Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-58607957929094905012016-12-06T22:45:44.517-08:002016-12-06T22:45:44.517-08:00Great post. Thanks. Great post. Thanks. 常青树https://www.blogger.com/profile/09499813030229291918noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-17136789747140136352016-12-06T08:56:54.244-08:002016-12-06T08:56:54.244-08:00Thanks for a nicer version of gen_data func!
Regar...Thanks for a nicer version of gen_data func!<br />Regarding pam. The post is about optimization in general. pam is limited in functionality and do not handle side constraints. For a real life problem you can check my presentation I gave at EARL 2014 (https://goo.gl/YY1d31). Wit Jakuczunhttps://www.blogger.com/profile/11026188487219287866noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-82417787117758010962016-12-06T05:47:53.382-08:002016-12-06T05:47:53.382-08:00An absolutely killer-feature of Julia JuMP is the ...An absolutely killer-feature of Julia JuMP is the ability to define lazy-constraints via callbacks in a solver-independent manner. Many of the best algorithms for problems make use of such techniques (such as subtour-elimination in Travelling Salesman) and the ability to easily program them is something that I have not seen in any library other than JuMP.Adam Sardarhttps://www.blogger.com/profile/07517057242770665359noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-14996400736357874222016-12-06T04:22:54.917-08:002016-12-06T04:22:54.917-08:00Thank you for the comment about ompr package. This...Thank you for the comment about ompr package. This is exactly what I have been looking for! I will have to check it out in detail, but it looks very good.<br /><br />As for the speed - most of the time is spent in solver and the solver is the same in both cases.Bogumił Kamińskihttps://www.blogger.com/profile/06250268799809238730noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-73424061757968941182016-12-06T04:09:11.603-08:002016-12-06T04:09:11.603-08:00I really understand that glpkAPI has poor designed...I really understand that glpkAPI has poor designed interface in comparison with JuMP. But why R code is so unnecessary obfuscated and complicated in your example? Especially confusing is two type-castings to data.frame in 'gen_data'. In Julia you use matrix type for your data. Why don't you use matrix in R? Just for completeness - in real life R code for this tasks looks like this:<br /><br />gen_data = function(N, K, scaling) {<br /> alpha = 2 * pi * (1:N) / K<br /> sin_pi_K = 2 * sin(pi / K)<br /> x = matrix(rnorm(n = 2*N), nrow = N)<br /> x + cbind(cos(alpha), sin(alpha)) * scaling / sin_pi_K<br />}<br /><br />set.seed(1)<br />x = gen_data(200, 4, 5)<br /><br />library(cluster)<br />lapply(2:6, function(k) pam(x, k = k))Gregory Deminhttps://www.blogger.com/profile/09827917062640566527noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-42416569071838202612016-12-06T02:29:01.511-08:002016-12-06T02:29:01.511-08:00Great post! Nice to have the comparisons, pretty c...Great post! Nice to have the comparisons, pretty clear that Julia's syntax really wins the day here, is there much of a speed difference? <br /><br />Also, have you seen the ompr package (https://github.com/dirkschumacher/ompr) that Dirk Schumacher is working on? It's inspired by the JuMP project.Surferhttps://www.blogger.com/profile/10433810043220326846noreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-50916168481686147142016-09-20T00:08:02.743-07:002016-09-20T00:08:02.743-07:00Use ROCR package in R
use function performance
ex...Use ROCR package in R<br />use function performance<br /><br />example<br /><br />library(ROCR)<br />pred<-prediction(actual,predicted)<br />perf<-performance(pred,"tpr","fpr")<br />plot(perf,col="red")<br />abline(0,1, lty = 8, col = "grey")<br /><br />auc<-performance(pred,"auc")<br />unlist(auc@y.values)Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4946490806848569840.post-76971942695951151442016-08-02T18:37:03.730-07:002016-08-02T18:37:03.730-07:00I was looking for something like a singleton patte...I was looking for something like a singleton pattern in R to keep objects that are expensive to load and this is the best approach I found so far. Thanks for that! <br />Previously I was playing with globals but this is messy and lintr was complaining about it.<br /><br />This is how I used it to load the biomaRt object.<br /><br />load_biomart <- function() {<br /> ensembl <- attr(load_biomart, "cached_ensembl")<br /> if (is.null(ensembl)) {<br /> ## Loads ensembl biomart for Homo sapiens ensembl = useMart('ensembl',dataset='hsapiens_gene_ensembl')<br /> futile.logger::flog.info("Loading biomart for the first time")<br /> ensembl_mart <- biomaRt::useMart("ENSEMBL_MART_ENSEMBL", host = "www.ensembl.org")<br /> dataset <- "hsapiens_gene_ensembl"<br /> ensembl <- biomaRt::useDataset(dataset, mart = ensembl_mart)<br /> attr(load_biomart, "cached_ensembl") <<- ensembl<br /> } else {<br /> futile.logger::flog.info("Returning cached biomart")<br /> }<br /> ensembl<br />}Chinohttps://www.blogger.com/profile/14034371996494176728noreply@blogger.com