mManipulating from particular column

Gene        T0          T1      T2    T3    T4   

'gene1'     'u'     'u'  'd'    'd'  'u'        
'gene2'     'u'    'u'  'u'    'u'  'u'    
'gene3'     'd'    'd'  'u'    'd'  'd'    
'gene4'     'd'    'u'  'd'    'u'  'd'   

Taking T0 as reference for gene1 i need to display values which have 'u' for gene1(T0,T1,T4) have 'u',and corresponding column must be displayed for all genes. for ex gene1

'gene1'     'u'     'u'  'u'        
'gene2'     'u'    'u'  'u'    
'gene3'     'd'    'u'  'd'    
'gene4'     'd'    'u'  'd' 

for gene3 value is d(T0) so d must be taken

'gene3'     'd'    'd'      'd'  'd'    
'gene4'     'd'    'u'      'u'  'd'   

I want to display for all gene i have totally 1000 genes ,please help

Answer 1

[Miro]

Something like this:

[n,~] = size(Data);
sgene = 'u'; %searched letter
referencepos = 1;   %position of the Row in the Cell matrix of reference gene
s=Data{referencepos,:} == sgene;
Output = {};
Datatmp2 = {};
for i = 1 : n
 Datatmp = cell2mat(Data(i,2:end));
 Datatmp = Datatmp(s);
 for k = 1 : length(Datatmp)
   Datatmp2 = [Datatmp2 Datatmp(k)];
 end 
 Datatmp2 = {Data(i,1) Datatmp2}
 Output = [Output;Datatmp2 ]
end

[Pat]

Miro i am getting error

Error using == Too many input arguments.

Error in Untitled (line 11) s=Data{referencepos,:} == sgene;

[Walter Roberson]

s = [Data{referencepos,:}] == sgene;

[Pat]

Miro can you please tell which variable to look for output ,because in Output variable in your above code i did not get the result as i desired.can u please tell

[Miro]

ok so the s should be like this

s=cell2mat(Data(referencepos,:)) == sgene;

[Miro]

so if your Datacell is like this:

'gene1'     'u'     'u'  'd'    'd'  'u'        
'gene2'     'u'    'u'  'u'    'u'  'u'    
'gene3'     'd'    'd'  'u'    'd'  'd'    
'gene4'     'd'    'u'  'd'    'u'  'd'

and you want gene1 as reference you set

referencepos = 1; (for gene2 its referencepos = 2 etc.)

with the searched letter 'u', the output will be

'gene1'     'u'     'u'  'u'        
'gene2'     'u'    'u'  'u'    
'gene3'     'd'    'd'  'd'    
'gene4'     'd'    'u'  'd'

you should use Breakpoints in the script to make debugging. For example set a breakpoint within the loop and see how it changes per iteration. Good Luck

[Jan]

@Miro: Use STRCMP to compare cell strings with a string. CELL2MAT and == is not reliable, when the strings have more then 1 character.

[Miro]

what do you mean with reliable? You mean that it sometimes works and sometimes not? Usually i used to use strcmp too, but just rescently i discovered, that booleans also work, which usually are the fastest method to go for.