[julia-users] how can I create an incidence matrix

[Gerson Jr.]

Hi
I would like to create an incidence matrix.
I have a file with 3 columns, like:

id  x   y
A 22   2
B  4  21
C 21 360
D 26   2
E 22  58
F  2 347

And I want a matrix like (without col and row names):

  2 4 21 22 26 58 347 360
A 1 0  0  1  0  0   0   0
B 0 1  1  0  0  0   0   0
C 0 0  1  0  0  0   0   1
D 1 0  0  0  1  0   0   0
E 0 0  0  1  0  1   0   0
F 1 0  0  0  0  0   1   0

I have started the code like:

haps = readdlm("File.txt",header=true)         #read the file
hap1_2 = map(Int64,haps[1][:,2:end])          # create a subfile with 2 and 
3 column
ID = (haps[1][:,1])                                          #subfile with 
the 1 column
dic1 = Dict()

for (i in 1:21)
    dic1[ID[i]] = hap1_2[i,:]
end

X=[zeros(21,22)];                         #the original file has 21 rows 
and 22 columns 
hcat(ID,X)



The problem now is that I don't know how to fill the matrix with 1s in the 
specific columns as in the example above.  
I'm also not sure if I'm on the right way.

Any suggestion??

Thanks!