No attachment came through the mailserver. Did you follow the
directions in the posting guide regarding acceptable types of
attachments? I also do not see any code that would let us reproduce an
input process.
One method that might work is to offer the results of
dput(coefficients) as text within your message.
--
David
On Jun 4, 2009, at 12:22 PM, Katharina May wrote:
Hi there,
I trying to solve this problem for the whole day not going anywhere,
so I really hope maybe somebody can help
me in this community...
I've got an object coefficient2 which I want to plot in differerent
ways, with colors and labels added to the points,
but somehow there seems to be a problem if a value is NA within the
independent variable, resulting in false labels and false colors for
the points.
plot(coefficient2$intercept ~ coefficient2$average_height,
main="intercepts ::: height", ylab="intercepts", xlab="average height
per site [cm]", xlim=c(20,3020), ylim=c(-2,5), col=coefficient2$color)
highlight(coefficient2$intercept ~ coefficient2$average_height,
lbls=coefficient2$site_no,col="Red", cex = .6)
plot(coefficient2$intercept ~ coefficient2$average_dbh,
main="intercepts ::: dbh", ylab="intercepts", xlab="average dbh per
site [mm]", xlim=c(-10,360), ylim=c(-2,5),col=coefficient2$color )
highlight(coefficient2$intercept ~ coefficient2$average_dbh,
lbls=coefficient2$site_no,col="Red", cex = .6
If I create a temporary object for each plot excluding any NA values
for the x axis variable, somehow all points are displayed and the
labels are correct except for the color (I have e.g. no clue why some
points are red which I do not define at all and no ones are yellow
which I use several times).
#create temp container for all coefficients with average heights for
plotting
coef_avheight <- coefficient2[which(!
is.na(coefficient2$average_height)),]
plot(coef_avheight$intercept ~ coef_avheight$average_height,
main="intercepts ::: height", ylab="intercepts", xlab="average height
per site [cm]", xlim=c(20,3020), ylim=c(-2,5),
col=coef_avheight$color)
highlight(coef_avheight$intercept ~ coef_avheight$average_height,
lbls=coef_avheight$site_no,col="Red", cex = .6)
#create temp container for all coefficients with average dbh for
plotting
coef_avdbh<- coefficient2[which(!is.na(coefficient2$average_dbh)),]
plot(coef_avdbh$intercept ~ coef_avdbh$average_dbh, main="intercepts
::: dbh", ylab="intercepts", xlab="average dbh per site [mm]",
xlim=c(-10,360), ylim=c(-2,5),col=coef_avdbh$color )
highlight(coef_avdbh$intercept ~ coef_avdbh$average_dbh,
lbls=coef_avdbh$site_no,col="Red", cex = .6)
Maybe someone can explain me the color issue and the problem with the
NA values which results in a wrong labeling and to few
points being displayed? I'm new to R as you can guess and my code
isn't really elegant but I really cannot get faults within it...
Attached you can find the referred R object (coefficient2).
highlight requires library(NCStats)...
Thank you very, very much,
Katharina
______________________________________________
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.
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
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.
