l******s
发帖数: 3045 | 1 谢谢分享。不过会有time O(n)的可能么？
写了一个O(nln(n))的对付一下。
private static int jumpTimes(int[] leaves, int distance, int maxStep){
SortedDictionary dict = new SortedDictionary();
int i = 0;
for(i = 0; i < leaves.Length; i++)
if(!dict.ContainsKey(leaves[i])) dict[leaves[i]] = i;
int[,] lt = new int[dict.Count, 2];
i = 0;
foreach(var d in dict)
{ lt[i, 0] = d.Key; lt[i++, 1] = d.Value; }
for(i = 0; i < lt.Length; i++){
if(lt[i, 0] - (i == 0 ? 0 :... 阅读全帖 |
|
l******g
发帖数: 3 | 2 在用NDSolve解DE的时候, 那本厚厚的手册里给出一些参数选项, 比如: AccuracyGoal,
PrecisionGoal, MaxSteps, StartingStepSize, MaxStepSize.
我在用NDSolve求解一个Time Dependent Defferential Equation, 需要对MaxSteps和
StartingStepSize选择合适的值才能给出正确的结果. 但是我不理解MaxSteps,
StartingStepSize的含义和他们之间的关系. 因此总是做不对. 恳请各位高手帮我,
不胜感激 |
|
i*********7
发帖数: 348 | 3 day 8是不行的。
你可以数数。
Day maxstep
1 2
2 3
3 4
4 5
5 6
6 7
7 8
8 9
day 8最大只能到第九步 |
|
l**s
发帖数: 1 | 4 没有用过mathmatica,不负责任的泛泛的讲讲：）
因为有时间微分项，所以时间步长必然会影响效率
但很多问题在最初的瞬时解变化很大
所以最初的时候需要非常小的时间步长
随着解趋于稳定，后面的时间步长可以很大很大而不影响精度
但有时候解是振荡的，或者是个移动的波
所以最大时间步长始终不能太大
startingstepsize和maxsteps大概就起这么个控制作用
一般的ODEs和DAEs solver都有这个选项
著名的如dassl（fortran）和ode15s(matlab) |
|
n*****t
发帖数: 41 | 5 谢谢，文件如下。我只要读取最后一个数，就是那个likelihood的数值。
MERLIN 1.1.2 - (c) 2000-2007 Goncalo Abecasis
References for this version of Merlin:
Abecasis et al (2002) Nat Gen 30:97-101 [original citation]
Fingerlin et al (2004) AJHG 74:432-43 [case selection for
association studies]
Abecasis and Wigginton (2005) AJHG 77:754-67 [ld modeling, parametric
analyses]
Fingerlin et al (2006) Gen Epidemiol 30:384-96 [sex-specific maps]
Chen and Abecasis (2007) AJHG 81:913-26 [qtl assoc... 阅读全帖 |
|
