a bug of function plot (PR#4405)

From: xzhao@stat.wvu.edu [mailto:xzhao@stat.wvu.edu] 

Full_Name: Ximing Zhao
Version: R 1.7.1
OS: MacOs X
Submission from: (NULL) (64.136.27.51)

When I used function plot in R 1.7.1, I found two bugs.
My whole function is this:
*****************
drawexp<-function(X,x,c,b){
exp<-function28(X,x,c,b)
if(max(exp)>=max(dexp(x,rate=1))){
plot(x,exp, type="l",main= paste("Exponential ,n=",length(X),"
c=",round(c,2),"b=",b),sub="X=exp(rate=1)")
lines(x,dexp(x,rate=1),col="red",type="l") 
lines(density(X,width="SJ-dpi",n=1000,from=0.1,to=4),lty=2,col="blue")
}
else{
plot(x,dexp(x,rate=1),col="red",type="l",main= 
paste("Exponential ,n=",length(X)," 
c=",round(c,2),"b=",b),sub="X=exp(rate=1)")
lines(x,exp, type="l")
lines(density(X,width="SJ-dpi",n=1000,from=0.1,to=4),lty=2,col="blue")
}
text(3,0.35,"-- theoretical",col="red",adj=c(0,0))
text(3,0.3,"-- moment type", adj=c(0,0))
text(3,0.25, "-- Kernel type", col="blue", adj=c(0,0))
}

function28<-function(X,x,c,b){
est<-function(X,x,alfa){
fx<-(1/gamma(alfa))*((alfa-1)/x)*(alfa*X/x)^(alfa-1)*exp(-alfa*X/x)
fx
}
n<-length(X)
result<-colSums(sapply(x,function(x){est(X,x,alfa=(n^c)*log(n)
^b)}))/length(X)
result
}

X1<-rexp(100)

x1<-seq(0.1,4,length=100)

Run the function like this:

drawexp(X1,x1,0,0.5)

***************

For 
"lines(density(X,width="SJ-dpi",n=1000,from=0.1,to=4),lty=2,co
l="blue")",
only lty=1 works, lty=2 and others don't work. I tried the 
function in 1.6.2 before, it was alright. So I am positive 
it's a bug. Please check.