ChemmineR Help
0
0
Entering edit mode
Thomas Girke ★ 1.7k
@thomas-girke-993
Last seen 8 months ago
United States
Dear Danishuddin, As far as I can tell, the error points to a typo you have in your R code: it cannot find "sdfap" since the correct function name for computing atom pair descriptors is "sdf2ap". Also, if you have only 5000 molecules in your SD file then it may not be necessary to use the sdfStream function which is for processing very large numbers of molecules while minimizing the memory footprint. Importing your SD file directly into an SDFset object with "read.SDFset" will be sufficient. A simple clustering routine with atom pairs/fingerprints and the binning clustering algorithm would look like this: ## Import your own file library(ChemmineR) sdfset <- read.SDFset("your_SD_file") ## Or use sample set provided by ChemmineR for testing data(sdfsample); sdfset <- sdfsample ## Compute atom pairs apset <- sdf2ap(sdfset) # Will run for some minutes ## Or compute atom pair fingerprints fpset <- desc2fp(apset) ## Binning clustering clusters <- cmp.cluster(apset, cutoff = c(0.7, 0.8, 0.9)) ## The result should look similar to this clusters[1:4,] ids CLSZ_0.7 CLID_0.7 CLSZ_0.8 CLID_0.8 CLSZ_0.9 CLID_0.9 48 CMP48 2 48 2 48 2 48 49 CMP49 2 48 2 48 2 48 54 CMP54 2 54 2 54 2 54 55 CMP55 2 54 2 54 2 54 Please also make sure you are running a recent version of R and ChemmineR. The current version of ChemmineR is 2.11.8, which you can check by running the command sessionInfo() in your R session. Thomas On Thu, Nov 29, 2012 at 05:46:02AM +0000, danishuddin --- wrote: > Dear ChemmieR author > > I am a researcher working on chemoinformatics, recently i come to known about the ChemmineR package, > i found it very helpful for my study, firstly i would give you thanks form my side and also from the > side of resercher working on the chemoinforamtic field. Really, this is very helpful. Well, i am new > to this, thus i am facing some problems. I am using the file of 5000 compounds which are in sdf format > for clustering. while using this i got the error. The image of this error i have attach with this mail. > I request you to help me in coming up with this error. I will wait for you response. Please help me. > > > Thanking you > > > Regards > Danishuddin
Clustering ChemmineR Clustering ChemmineR • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 792 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6