What is structure, and how do I use it on Quarry?

On this page:

• About structure
• Using structure on Quarry

About structure

The program structure can use genotypes of individuals to form clusters that represent populations, to assign individuals to populations, and to describe genotypes of individuals in terms of populations from which they were derived.

The Department of Human Genetics at the University of Chicago distributes structure. For information about structure packages, see the department's Documentation for structure software: Version 2.2 (in PDF format). For sample data sets in structure format, see the department's sample data sets page.

Using structure on Quarry

At Indiana University, structure is available on Quarry. Accounts on Quarry are available to all IU students, faculty, and staff, as well as affiliated researchers. Instructional use is limited to courses approved by the Director for Research Technologies. All users must know the Quarry usage policies.

To use structure on Quarry:

1. Request a Quarry account.
2. Set up softenv and submit jobs.
3. Pass options to structure.

Request a Quarry account

To request a Quarry account, use the Account Management System (AMS); see At IU, if I already have some computing accounts, how do I get others?

For more, see Getting started on Quarry.

Set up softenv and submit jobs

1. To add structure to your softenv, edit your .soft file, adding the line +structure and executing a resoft.
   
   
2. Create a directory that contains all the datafile and paramfiles.
   
   
3. If your job will run for fewer than 20 minutes, call structure in the directory.

   If your job will take longer than 20 minutes, edit a PBS script in the directory:

   $ vi pbs.job #PBS -l ncpus=1 #PBS -l mem=1GB #PBS -l cput=1:00:00 #PBS -e log.err #PBS -o log.out #PBS -m abe cd $PBS_O_WORKDIR structure -m $PBS_O_WORKDIR/mainparams -e $PBS_O_WORKDIR/extraparams -i $PBS_O_WORKDIR/testdata1 -o $PBS_O_WORKDIR/testout
4. Submit pbs.job to run as a batch job: $qsub pbs.job

Output that structure would normally print to the screen is stored in the $PBS_O_WORKDIR/testout and $PBS_O_WORKDIR/log.out files.

Passing options

To run structure from the command line, you will need the following files:

Also, you will have to pass options to structure. If you have the rights to configure in your three files, you can turn an interactive structure program into a batch job by creating a file that contains the options that you need to pass. To do so, create a text file that contains your options, one option per line, or combine all the options on the command line and pipe them to structure.

For example, if the default option is Y, which you need to press and return to keep the structure program running, your option file would have only one line, containing the letter Y. If you name the text file optionfile, then the command would be:

Alternatively, you can echo the option pipe it to structure: /bin/echo "Y\n" | structure -m mainparams -e extraparams -i datafile

To pass a series of options, you have to use newline (\n) to separate them onto different lines:

Note: IU's involvement in the TeraGrid, and the presentation of this material is based upon work supported by the National Science Foundation under Grants No. 0833618, SCI451237, SCI535258, and SCI504075. Any opinions, findings and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

Also see:

• Bioinformatics support