我正在存储不同元素之间的联系。我想消除某种类型的元素并存储由消除的元素互连的元素的新联系。
想象一下这个问题。你有一个水分子与其他分子接触(如果接触是氢键，我的水周围可能还有 4 个其他分子)。如下图所示(A、B、C、D 是其他一些原子，点表示接触)。
A B
| |
H H
. .
O
/ \
H H
. .
C D
我有所有点的信息，我需要去除中心的水并创建描述 A-C、A-D、A-B、B-C、B-D 和 C-D 接触的记录。
目前，我在数据库中有以下结构:
表原子:
“id”整数主键，"amino"char(3) NOT NULL,(HOH 表示水或其他值)表联系人:
"acceptor_id"integer NOT NULL,(靠近我的氢的原子，这里是 C 或 D)"donor_id"integer NOT NULL,(此处为 A 或 B)"directness"char(1) NOT NULL,(这应该是 D 表示直接，W 表示水介导)编辑: 前面描述的案例中的数据看起来如何。
原子:
id|atom|amino
1 | O | HOH
2 | N | ARG <- atom A from image
3 | S | CYS <- B
4 | O | SER <- C
5 | N | ARG <- D
联系人:
donor_id|acceptor_id|directness
1 4 D
1 5 D
2 1 D
3 1 D
我需要从中制作
联系人:
donor_id|acceptor_id|directness
3 4 W <- B-C
3 5 W <- B-D
2 4 W <- A-C
2 5 W <- A-D
2 3 X <- A-B (These last two rows are escaping me,
4 5 X <- C-D there could be also row D-C, but not
both C-D and D-C. A char 'X' could
be used to mark "no donor/acceptor")
现在，我将查看所有具有 donor.amino = "HOH" 的联系人。在此示例中，这将从 C 和 D 中选择联系人。对于这些选定的联系人中的每一个，我查找具有与当前 donor_id 相同的 acceptor_id 的联系人选定的联系人。根据这些信息，我创建了新联系人。最后，我删除了所有进出 HOH 的联系人。
这样，我显然无法创建 C-D 和 A-B 联系人(其他 4 个都可以)。
如果我尝试类似的方法 - 尝试找到具有相同 donor_id 的两个联系人，我最终会得到重复的联系人(C-D 和 D-C)。
我梦想着一些一页长的 SQL 查询，它只检索这六个需要的行。 :-)
但是，欢迎任何其他想法。
最好保留有关谁是捐赠者的信息(如果可能)，但并非绝对必要。
非常感谢所有阅读此问题的人。
最佳答案
您的解释有一个困难。
您从一个有向图开始，其中每条边代表连接 X=>Y，其中 X 是供体，Y 是受体。 atoms 表是该图的 SQL 表示。
您似乎想要的是未定向的东西。因此，链接 X-Y 意味着 X 和 Y 通过水分子(当然或其他一些物种)连接，但 X 和 Y 都可以是供体或受体。出于这个原因，您的最后一张表有歧义(您注意到)，因此某些链接可能会以任何一种方式出现。在我看来，这意味着您最后一列中的列标题 donor_id 和 acceptor_id 没有您所解释的任何含义。这当然可能是我的困惑。
如果你想要的只是一个包含所有 6 个链接的表格，每个链接为一行，但不要太担心跟踪捐赠者/接受者的事情，那么这在 sqlite3 中对我有用:
create temporary view hoh_view as
select donor_id as id, atoms.id as hoh_id from contacts, atoms
where acceptor_id=atoms.id and atoms.amino='HOH'
union select acceptor_id as id, atoms.id as hoh_id from contacts, atoms
where donor_id=atoms.id and atoms.amino='HOH';
select a.id, b.id from hoh_view as a, hoh_view as b
where a.id > b.id and a.hoh_id=b.hoh_id;
我使用临时 View 使事情更清楚的地方。如果您愿意，可以将对 hoh_view 的每个引用替换为第一个查询，从而将所有这些放入一个大查询中。我觉得有点讨厌，可能有办法整理它。
如果您确实想跟踪供体/受体关系，您需要解释当两个氨基酸都是受体或供体时您如何决定做什么(即示例中的最后两行)。
如果那不符合您的要求，那么也许我可以修复它，让它符合您的要求。
关于python - 优化存储在 SQLite 中的数据 - 如何加入多个联系人？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/2904205/
我正在学习如何使用Nokogiri，根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div
总的来说，我对ruby还比较陌生，我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础，我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用
我试图在一个项目中使用rake，如果我把所有东西都放到Rakefile中，它会很大并且很难读取/找到东西，所以我试着将每个命名空间放在lib/rake中它自己的文件中，我添加了这个到我的rake文件的顶部:Dir['#{File.dirname(__FILE__)}/lib/rake/*.rake'].map{|f|requiref}它加载文件没问题，但没有任务。我现在只有一个.rake文件作为测试，名为“servers.rake”，它看起来像这样:namespace:serverdotask:testdoputs"test"endend所以当我运行rakeserver:testid时
作为我的Rails应用程序的一部分，我编写了一个小导入程序，它从我们的LDAP系统中吸取数据并将其塞入一个用户表中。不幸的是，与LDAP相关的代码在遍历我们的32K用户时泄漏了大量内存，我一直无法弄清楚如何解决这个问题。这个问题似乎在某种程度上与LDAP库有关，因为当我删除对LDAP内容的调用时，内存使用情况会很好地稳定下来。此外，不断增加的对象是Net::BER::BerIdentifiedString和Net::BER::BerIdentifiedArray，它们都是LDAP库的一部分。当我运行导入时，内存使用量最终达到超过1GB的峰值。如果问题存在，我需要找到一些方法来更正我的代
关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。
Rails2.3可以选择随时使用RouteSet#add_configuration_file添加更多路由。是否可以在Rails3项目中做同样的事情？ 最佳答案 在config/application.rb中:config.paths.config.routes在Rails3.2(也可能是Rails3.1)中，使用:config.paths["config/routes"] 关于ruby-on-rails-Rails3中的多个路由文件，我们在StackOverflow上找到一个类似的问题
给定这段代码defcreate@upgrades=User.update_all(["role=?","upgraded"],:id=>params[:upgrade])redirect_toadmin_upgrades_path,:notice=>"Successfullyupgradeduser."end我如何在该操作中实际验证它们是否已保存或未重定向到适当的页面和消息？ 最佳答案 在Rails3中，update_all不返回任何有意义的信息，除了已更新的记录数(这可能取决于您的DBMS是否返回该信息)。http://ar.ru
我在我的项目目录中完成了compasscreate.和compassinitrails。几个问题:我已将我的.sass文件放在public/stylesheets中。这是放置它们的正确位置吗？当我运行compasswatch时，它不会自动编译这些.sass文件。我必须手动指定文件:compasswatchpublic/stylesheets/myfile.sass等。如何让它自动运行？文件ie.css、print.css和screen.css已放在stylesheets/compiled。如何在编译后不让它们重新出现的情况下删除它们？我自己编译的.sass文件编译成compiled/t
我有多个ActiveRecord子类Item的实例数组，我需要根据最早的事件循环打印。在这种情况下，我需要打印付款和维护日期，如下所示:ItemAmaintenancerequiredin5daysItemBpaymentrequiredin6daysItemApaymentrequiredin7daysItemBmaintenancerequiredin8days我目前有两个查询，用于查找maintenance和payment项目(非排他性查询)，并输出如下内容:paymentrequiredin...maintenancerequiredin...有什么方法可以改善上述(丑陋的)代
我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚