[R] Density plots

[David L Carlson]

How about providing the data for at least one subject/8 occasions using
dput(dataframe)? This line at the bottom of your message: "provide
commented, minimal, self-contained, reproducible code" is important. Nothing
you sent allows us to reproduce what you are doing and suggest ways to
improve it. Looking at your code it seems interesting that the 8 WSEQ values
are represented as binary. It may be easier to automate if that were
converted to a factor (but perhaps it already is a factor, there is no way
to tell from what you have provided).

Question 2 refers to a function that does not seem to exist -
sm.densityplot.compare - plus you copy the error message but not your
function call. How can we possibly tell you what you did wrong? Is this
really function sm.density.compare() in package sm? 

----------------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77843-4352

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Ayyappa Chaturvedula
> Sent: Friday, August 03, 2012 10:37 AM
> To: r-help at r-project.org
> Subject: [R] Density plots
> 
> Dear group,
> 
> I need help on two problems:
> 
> 1.  I am trying to plot density plots for each individual in 8
> occasions.
> I can do this by subject wiht the code below:
> par(mfrow=c(4,2))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==0]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==1]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==10]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==11]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==100]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==101]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==110]))
> plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WS
> EQ==111]))
> 
> I have total of 100 subjects and I want to automate this and create for
> every subject.
> 
> 2.  I want to also plot each subject in the same plot for all 8
> scenarios
> instead of the way I am doing above.  Could you please help me with the
> coding?
> 
> I tried sm.densityplot.compare but it is giving me an error:
> missing data are removed
> missing data are removed
> Error in if (opt$nbins > 0) { : missing value where TRUE/FALSE needed
> In addition: Warning message:
> In cbind(X, group) :
>   number of rows of result is not a multiple of vector length (arg 2)
> 
> I appreciate your help.
> 
> Regards,
> Ayyappa
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.