fork download
copy 
  1. from scipy.stats import multinomial
  2. import pandas as pd
  3.  
  4. n = 10
  5. p1 = 3/100
  6. p2 = 12/100
  7. p = [p1, p2, None]
  8.  
  9. df = pd.DataFrame(
  10.   { f'SSR{n1}':
  11.       [ f'{multinomial.pmf(x=[n1, n2, n-(n1+n2)], n=n, p=p)*100:.4f}' if n1+n2 <= n else None
  12.         for n2 in range(0, n+1)
  13.       ]
  14.     for n1 in range(0, n+1)
  15.   }, index=[ f'SR{n2}' for n2 in range(0, n+1) ] )
  16.  
  17. print(df.to_csv(sep='\t'))
Success #stdin #stdout 0.78s 81820KB
comments ()
stdin
copy 
Standard input is empty
stdout
copy 
	SSR0	SSR1	SSR2	SSR3	SSR4	SSR5	SSR6	SSR7	SSR8	SSR9	SSR10
SR0	19.6874	6.9485	1.1036	0.1039	0.0064	0.0003	0.0000	0.0000	0.0000	0.0000	0.0000
SR1	27.7940	8.8287	1.2464	0.1026	0.0054	0.0002	0.0000	0.0000	0.0000	0.0000	
SR2	17.6574	4.9856	0.6159	0.0435	0.0019	0.0001	0.0000	0.0000	0.0000		
SR3	6.6475	1.6423	0.1739	0.0102	0.0004	0.0000	0.0000	0.0000			
SR4	1.6423	0.3478	0.0307	0.0014	0.0000	0.0000	0.0000				
SR5	0.2782	0.0491	0.0035	0.0001	0.0000	0.0000					
SR6	0.0327	0.0046	0.0002	0.0000	0.0000						
SR7	0.0026	0.0003	0.0000	0.0000							
SR8	0.0001	0.0000	0.0000								
SR9	0.0000	0.0000									
SR10	0.0000
https://ideone.com/5kd7ED
language:
Python 3 (python  3.12)
created:
2 years ago
visibility:
secret

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!