R help - Sep 2011 - Dealing with NA's in a data matrix

hello... I am trying to get this code to work, but as I get to the predict
command, it displays an error due to the length of the data sets from the
removal of the NA's.   Here is the data, and the code that I am using so
far, if you run it, you'll see the error pop up.... please help me to get
around this problem. Thanks in advance.

Location Dist. Size
low1 .5 10.5
low2 .5 23
low3 .5 NA
low4 .5 NA
mid1 3 15
mid2 3 11.5
mid3 3 15
mid4 3 NA
high1 6 12.5
high2 6 20
high3 6 21
high4 6 22
wall1 10 13
wall2 10 12.5
wall3 10 13
wall4 10 12

d<-read.table("dilomaaethiops.txt",header=T)
d
str(d)
range(d$Dist.)
s<-d$Size
range(s)
s[s<.1]<-NA
range(s,na.rm=T)
summary(d)
length(s)
length(d$Dist.)
plot(s~d$Dist.,xlab="Distance from Shoreline (m)",ylab="D.
aethiops Size
(mm)", pch=7)
a1<-glm(s~d$Dist.,family="quasipoisson",data=d)
scatter.smooth(d$Dist.,s,ylab="D. aethiops Length
(mm)",xlab="Distance from
Shoreline (m)",col="red")
summary(a1)
a2<-aov(s~d$Dist.,data=d)
summary(a2)
shapiro.test(a2$residuals)
qqnorm(a2$residuals)
qqline(a2$residuals)
fligner.test(s~d$Dist.,data=d)
kruskal.test(s~d$Dist.,data=d)
df<-df.residual(a1)
t<-qt(.975,df)
k<-data.frame(s=seq(10.5,max(s,na.rm=T),1))
p<-predict(a1,k,se=T)
est<-exp(p$fit)
low<-exp(p$fit-t*p$se.fit)
hight<-exp(p$fit-t*p$se.fit)
plot(d$Dist.,s,ylim=c(min(s),max(s)),xlab="Distance from Shorline
(m)",
ylab="Length of D. aethiops (mm)", pch=7)
lines(k$s,est,lty=1,col="red")
lines(k$s,low,lty=2,col="blue")
lines(k$s,high,lty=2,col="blue")

	[[alternative HTML version deleted]]

Please provide the data as self-contained code, as requested in the posting
guide, so that helpers can directly paste it into R. Alternatively, you can
provide the dilomaaethiops.txt.

Best,
Daniel


Watching19 wrote:
> 
> hello... I am trying to get this code to work, but as I get to the predict
> command, it displays an error due to the length of the data sets from the
> removal of the NA's.   Here is the data, and the code that I am using
so
> far, if you run it, you'll see the error pop up.... please help me to
get
> around this problem. Thanks in advance.
> 
> Location Dist. Size
> low1 .5 10.5
> low2 .5 23
> low3 .5 NA
> low4 .5 NA
> mid1 3 15
> mid2 3 11.5
> mid3 3 15
> mid4 3 NA
> high1 6 12.5
> high2 6 20
> high3 6 21
> high4 6 22
> wall1 10 13
> wall2 10 12.5
> wall3 10 13
> wall4 10 12
> 
> d<-read.table("dilomaaethiops.txt",header=T)
> d
> str(d)
> range(d$Dist.)
> s<-d$Size
> range(s)
> s[s<.1]<-NA
> range(s,na.rm=T)
> summary(d)
> length(s)
> length(d$Dist.)
> plot(s~d$Dist.,xlab="Distance from Shoreline (m)",ylab="D.
aethiops Size
> (mm)", pch=7)
> a1<-glm(s~d$Dist.,family="quasipoisson",data=d)
> scatter.smooth(d$Dist.,s,ylab="D. aethiops Length
(mm)",xlab="Distance
> from
> Shoreline (m)",col="red")
> summary(a1)
> a2<-aov(s~d$Dist.,data=d)
> summary(a2)
> shapiro.test(a2$residuals)
> qqnorm(a2$residuals)
> qqline(a2$residuals)
> fligner.test(s~d$Dist.,data=d)
> kruskal.test(s~d$Dist.,data=d)
> df<-df.residual(a1)
> t<-qt(.975,df)
> k<-data.frame(s=seq(10.5,max(s,na.rm=T),1))
> p<-predict(a1,k,se=T)
> est<-exp(p$fit)
> low<-exp(p$fit-t*p$se.fit)
> hight<-exp(p$fit-t*p$se.fit)
> plot(d$Dist.,s,ylim=c(min(s),max(s)),xlab="Distance from Shorline
(m)",
> ylab="Length of D. aethiops (mm)", pch=7)
> lines(k$s,est,lty=1,col="red")
> lines(k$s,low,lty=2,col="blue")
> lines(k$s,high,lty=2,col="blue")
> 
> 	[[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.
> 
--
View this message in context:
http://r.789695.n4.nabble.com/Dealing-with-NA-s-in-a-data-matrix-tp3791462p3792737.html
Sent from the R help mailing list archive at Nabble.com.