r - Subset data frame based on unique combination of multiple conditions -


i cant seem find answer through search on so. i'm trying select subset of data.frame based on 4 conditions (lon1, lon2, lat1 , lat2). have huge dissimilarity matrix has been vectorized , sites (lon1, lon2, lat1 , lat2) cbind it. here example data frame: 

out1 <- data.frame(lon1 = sample(1:10), lon2 = sample(1:10),                     lat1 = sample(1:10), lat2 = sample(1:10),                     dissimilarity = sample(seq(0,1,.1),10)) > out1      lon1   lon2    lat1 lat2 dissimilarity 1     2      6      4      4           0.6 2     4      2      1      3           1.0 3    10      9      2      6           0.0 4     3      1     10      8           0.5 5     9      5      9      1           0.8 6     5      7      5      9           0.9 7     1      8      6      7           0.2 8     8      3      8      5           0.7 9     7      4      3     10           0.3 10    6     10      7      2           0.1  out2 <- out1[c(2,5,6,8),]     lon1 lon2 lat1 lat2 dissimilarity 1     4   2   1      3           1.0 2     9   5   9      1           0.8 3     5   7   5      9           0.9 4     8   3   8      5           0.7

i tried using %in% function few times in manner:

test <- out1[(out1$lon1 %in% out2$lon1) & (out1$lon2 %in% out2$lon2) &               (out1$lat1 %in% out2$lat1) & (out1$lat2 %in% out2$lat2), ]

this seems work basic example provide here. but, when apply huge data frame (with many lat , lons repeated) larger subset unqiue combinations require. assume because match function in %in% can match vector. it's matching condition1 & condition2 & condition3 & condition4 , returning results gives subset same orginal out1. want case when 4 values same row. way i'll subset of data pairwise dissimilarities i'm interested in.

any ideas on how subset row based on unique combination of 4 variables appreciated.

i think you're looking for. want duplicated function returns you're expecting. 

out1[duplicated(rbind(out2, out1)[, 1:4])[-seq_len(nrow(out2))], ]

how work? first rbind out2 , out1. call duplicated on it. columns in out2 , in out1 marked true in out1. because first occurrence on out2 , not duplicated there. second time finds entry, in out1 , it'll know there has been row before. so, it'll mark duplicated. have duplicated entries. subset elements of out1 removing first n elements n = nrow(out1). subset using logical vector on out1.

you can go through explanation , run code step step follow-up. here's break-down version working out logic.

tt <- rbind(out2, out1) tt.dup <- duplicated(tt[, 1:4)] # marks duplicate rows in out1 1st 4 cols tt.dup <- tt.dup[-seq_len(nrow(out2))] # remove out2 entries (first n) out1[tt.dup, ] # index true/duplicated elements out1

Comments

Post a Comment

Popular posts from this blog

php - cannot display multiple markers in google maps v3 from traceroute result -

i'm newbie in php mysql , google maps v3. have problems in google maps v3. made traceroute web based , map traceroute's result in google maps v3 final project of college education. have try many tutorials didn't work. have tried tutorials in forum nothing worked. problem cannot display multiple markers got location data database. maybe because i'm newbie , not know how do. confuse looping show markers. here code <?php error_reporting(e_all ^ (e_notice)); ini_set('max_execution_time', 360); $enable_log_user = false; global $ip, $host_name, $host_ip; $host = @$_post['host']; $trace = @$_post['trace']; $self = $_server['php_self']; include("phpsqlajax_dbinfo.php"); $connection = mysql_connect ('127.0.0.1', $username, $password); if (!$connection) { die('not connected : ' . mysql_error());} $db_selected = mysql_select_db($database, $connection); if (!$db_selected) { die ('can\'t use db : ...
Read more

c# - DetailsView in ASP.Net - How to add another column on the side/add a control in each row? -

i have detailsview control on page used edit various fields of record, works in respect. i looking way add 1 column (and if works, why not more) right, absolutely read-only, show same fields of record comparison purposes. i aware there no obvious way such thing out of box detailsview . have looked other controls (transposed gridview , recommended formview , listview ), nothing satisfies. have special data binding setup using detailsview , can't out of without losing features. anyone on how "hack in" additional columns (for display only) on detailsview ? the solution have now, use second detailsview , visible set false in aspx. in code, make sure databind hidden detailsview hosts data third column first, initial detailsview named itemdetails . and in item created event, pass third column html rendering of hidden controls (in last code block) : protected void itemdetails_itemcreated(object sender, eventargs e) { if (dataite...
Read more

javascript - firefox memory leak -

i have following script (i use raphaeljs also) draws circle, waits 1 second, draws circle again, etc (you can copy paste - functional, need include raphael raphaeljs.com): <div id="holder"> </div> <script src="raphael.js" type="text/javascript"></script> <script> var paper = raphael(10, 50, 320, 200); var rad = math.pi / 180; var finangle = 0; function sector(cx, cy, r, startangle, endangle, params) { var x1 = cx + r * math.cos(-startangle * rad), x2 = cx + r * math.cos(-endangle * rad), y1 = cy + r * math.sin(-startangle * rad), y2 = cy + r * math.sin(-endangle * rad); return ["m", cx, cy, "l", x1, y1, "a", r, r, 0, +(endangle - startangle > 180), 0, x2, y2]; } var path = sector(100,100,50,0,180); var arc = paper.path(path); function draw() { path = sector(100,100,50,0,finangle); ...
Read more