[Biojava-l] Create DP object from profileHMM class file
Tisanai
T.Krisanathamakul@dcs.shef.ac.uk
Mon, 7 Oct 2002 14:39:50 +0100
Hi
By this code I would like to create DP object from several phmm file.
for(int i=0;i<md_out_lst.align.length;i++){
String model_out_name = md_out_lst.align[i];
File md_file = new File(model_out_name);
FileInputStream fis_md = new FileInputStream(md_file);
ObjectInputStream ois_md = new ObjectInputStream(fis_md);
ProfileHMM model = (ProfileHMM) ois_md.readObject();
ois_md.close();
dp[i] = DPFactory.DEFAULT.createDP(model);
}
I found that it always stuck at the second file (i=2). If there is only one file in my list this code will work fine. But if there is more than one file in the list when it try to create the second dp object (dp[1]). This kind of error will shown out:
org.biojava.bio.BioError: State d-15 is known in states but is not listed in the transFrom table
at org.biojava.bio.dp.SimpleMarkovModel.transitionsFrom(SimpleMarkovModel.java:227)
at org.biojava.bio.dp.DP$HMMOrderByTransition.transitionsTo(DP.java:599)
at org.biojava.bio.dp.DP$HMMOrderByTransition.compare(DP.java:586)
at org.biojava.bio.dp.DP.stateList(DP.java:123)
at org.biojava.bio.dp.DP.update(DP.java:353)
at org.biojava.bio.dp.onehead.SingleDP.update(SingleDP.java:49)
at org.biojava.bio.dp.DP.<init>(DP.java:377)
at org.biojava.bio.dp.onehead.SingleDP.<init>(SingleDP.java:41)
at org.biojava.bio.dp.DPFactory$DefaultFactory.createDP(DPFactory.java:53)
at T_Zscore.main(T_Zscore.java:64)
How can I fix my code?
Thank
Tisanai