```
# NOT RUN {
data(codcluslen)
temp<-codcluslen[codcluslen$number>0,]
temp$station<-c(paste(temp$region,"-",temp$tow,sep=""))
total<-aggregate(temp$number,list(temp$station),sum)
names(total)<-c("station","total")
temp<-merge(temp,total,by.x="station",by.y="station")
newdata<-data.frame(NULL)
for(i in 1:as.numeric(length(temp[,1]))){
for(j in 1:temp$number[i]){
newdata<-rbind(newdata,temp[i,])
}
}
newdata<-newdata[,-c(5)]
clus.t.test(popchar=newdata$length,cluster=newdata$station,
group=newdata$region,rho=0.72,
alpha=0.05,alternative="two.sided")
# }
```

Run the code above in your browser using DataLab