On 03/22/2012 03:22 PM, ens wrote:
I have a matrix of p-values for for each explanatory variable. Each row is an
area of the response variable and each column is an explanatory variable.
e.g.
PSA pval_DOY pval_PDSIconcurrent pval_PDSIantecedent_annual_average
pval_TMAXanomaly pval_FM100anomaly
1 NC06 0.96747495 0.6092668 0.53353019
0.9301115 0.99801334
2 NC04 0.04699659 0.2759152 0.07024752
0.6046828 0.03224094
3 NC01 0.71437394 0.9979173 0.85296024
0.9977558 0.99833623
4 NC08 0.67315904 0.9970511 0.51756714
0.7809994 0.99626038
5 NC07 0.55221280 0.5784208 0.43975219
0.3669491 0.34898877
6 NC05 0.52089881 0.7191645 0.91972153
0.4487460 0.94922430
I want to create a visual display such that instead of #s for the p-values,
I have a circle sized to represent the size of the p-value.
symbolys() is what I've found to do this, but it isn't clear to me how to
format the input for the function.
Any suggestions?
Hi ens,
You might be interested in the "size_n_color" function (plotrix). I'll
be uploading v3.4 shortly, and you can use that to do it like this:
psa.mat
pval_DOY pval_PDSIconcurrent pval_PDSIaaa pval_TMAXanomaly
NC06 0.96747495 0.6092668 0.53353019 0.9301115
NC04 0.04699659 0.2759152 0.07024752 0.6046828
NC01 0.71437394 0.9979173 0.85296024 0.9977558
NC08 0.67315904 0.9970511 0.51756714 0.7809994
NC07 0.55221280 0.5784208 0.43975219 0.3669491
NC05 0.52089881 0.7191645 0.91972153 0.4487460
pval_FM100anomaly
NC06 0.99801334
NC04 0.03224094
NC01 0.99833623
NC08 0.99626038
NC07 0.34898877
NC05 0.94922430
size_n_color(size=psa.mat,col="black",xaxlab=rep("",5),
yaxlab=rownames(psa.mat))
staxlab(1,at=seq(2,10,by=2),colnames(psa.mat))
Jim
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
