This code consistenly segfaults for me. Can someone please take a look
and tell me if the problem is due to something I am doing or is there a
problems with the dwt (idwt) functions in the waveslim library.
Thanks
tom
library(waveslim)
data<-c(936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,936.944,916.994666666667,909.592,941.8,961.933333333333,987.8,969.2,934.866666666667,912.8,913.266666666667,928.2,949.533333333333,933.4,948.533333333333,929.866666666667,935.2,937,934.666666666667,923.733333333333,911.8,937.066666666667,921.133333333333,914.6,941.933333333333,961.933333333333,950.333333333333,925.6,921.066666666667,916.2,927.266666666667,921.266666666667,938.066666666667,927.133333333333,949.866666666667,938.533333333333,934.933333333333,945.8,962.466666666667,967.266666666667,965,935.533333333333,935,927.933333333333,919.733333333333,930.466666666667,948,945.733333333333,972.933333333333,942.066666666667,921.533333333333,921.2,920.4,956.333333333333,945.666666666667,971.4,980,965.6,964.133333333333,960.066666666667,963,932.866666666667,952.466666666667,938.8,924.8,925.466666666667,949.4,931.466666666667,924.066666666667,926.866666666667,931.133333333333,937.266666666667,961.333333333333,952.2,958.133333333333,964,941.733333333333,936.133333333333,936.2,955.866666666667,964.6,962.333333333333,946.866666666667,971.533333333333,952.666666666667,954.933333333333,959.866666666667,969.2,970.266666666667,954.8,974.266666666667,947.133333333333,941.466666666667,949.066666666667,974.333333333333,943,929.333333333333,933.666666666667,926.666666666667,939.933333333333,936.6,952.4,954.6,947.466666666667,945.533333333333,948.733333333333,950.066666666667,941.666666666667,926.466666666667,934.333333333333,948.933333333333,946.733333333333,933.266666666667,913.266666666667,924.466666666667,937.866666666667,980.6,975.6,954.266666666667,943.266666666667,954,948.133333333333,947.066666666667,915,924.933333333333,938.2,942.666666666667,941.266666666667,956.066666666667,964.066666666667,958.666666666667,947.066666666667,940.4,926.133333333333,931.066666666667,940.933333333333,956.466666666667,958.2,941.866666666667,941.4,935.066666666667,943.933333333333,939,948.666666666667,964.133333333333,996.533333333333,957.466666666667,945.733333333333,922.533333333333,919.2,931.733333333333,952.2,921.066666666667,933.466666666667,948.666666666667,950.8,954,958.333333333333,942.6,936.266666666667,936.8,938.266666666667,943.666666666667,936.533333333333,920.133333333333,902.266666666667,933,948.733333333333,924.866666666667,941.866666666667,943.8,961.533333333333,968.333333333333,984.466666666667,949.4,933.266666666667,931.266666666667,934.533333333333,937.066666666667,959.533333333333,942.4,953.266666666667,975.866666666667,983.6,966.2,953.866666666667,938.066666666667,936.666666666667,931.733333333333,926.866666666667,965.8,970.133333333333,944.6,956.8,937.266666666667,933.066666666667,947.333333333333,935.733333333333,944.933333333333,946,943.6,969.733333333333,965.8,979.2,973.533333333333,970.666666666667,960.533333333333,939.266666666667,929.8,919.266666666667,925.333333333333,913.333333333333,924.266666666667,927.666666666667,936.6,919.066666666667,919.333333333333,936.2,954,938.066666666667,947.266666666667,945.066666666667,943.8,941.733333333333,942.133333333333,935.8,940.066666666667,942.466666666667,950.6,932.6,923.866666666667,920.2,916.6,939.866666666667,922.2,941.066666666667,927.466666666667,916.6,914.4,924.533333333333,958.533333333333,960.333333333333,932.733333333333,937,971.733333333333,981.866666666667,966.666666666667,949.333333333333,930.733333333333,918.933333333333,917.533333333333,948.133333333333,941.066666666667,952.8,966.533333333333,971.333333333333,962.866666666667,948.333333333333,970.066666666667,952.333333333333,937.2,951.933333333333,928.133333333333,937,948.066666666667,944.066666666667,963.733333333333,960.533333333333,931.266666666667,958.933333333333,971.4,969.2,963.133333333333,963.933333333333,965.133333333333,963,962.066666666667,938.6,945.733333333333,938.933333333333,936.8,954.2,961.133333333333,973.466666666667,989.933333333333,959.4,946.866666666667,959.733333333333,953.666666666667,953.8,945.466666666667,954.533333333333,957.6,951.533333333333,954.6,959.733333333333,943.466666666667,938.333333333333,954,966.466666666667,944.733333333333,941.666666666667,936.733333333333,938.933333333333,930.2,947.933333333333,927.333333333333,923.866666666667,908.533333333333,942.466666666667,950.066666666667,931.266666666667,935.2,941.933333333333,941.4,944.866666666667,956.8,935.266666666667,945.333333333333,963.666666666667,968.8,962.4,942.2,953.666666666667,932.8,920,921.2,933.933333333333,918.733333333333,933.8,945.933333333333,964.2,968.733333333333,961.666666666667,962.2,978.933333333333,973.866666666667,953.6,957.133333333333,940.333333333333,946.733333333333,950.733333333333,938.2,942.933333333333,951.866666666667,934.4,939.266666666667,936,916.533333333333,921.133333333333,913.2,903.466666666667,920.866666666667,920.466666666667,919.2,958.133333333333,954.866666666667,965.8,960.066666666667,943,970.4,976.8,967.933333333333,960.333333333333,934.6,953.466666666667,949.266666666667,934.4,936.666666666667,949.933333333333,931.933333333333,941.666666666667,957.266666666667,939.933333333333,925.733333333333,899.066666666667,932.733333333333,984.266666666667,969.866666666667,962.733333333333,945.333333333333,939.133333333333,932.8,933.666666666667,943.8,941.066666666667,937.2,942.8,938.333333333333,940.733333333333,938.466666666667,930.066666666667,941,937.066666666667,941.333333333333,954,937.733333333333,927.666666666667,918.933333333333,918.066666666667,938.666666666667,924.4,931,931.733333333333,922.866666666667,921.6,933.4,929.866666666667,926.066666666667,949.266666666667,967.866666666667,949.933333333333,947.666666666667,937.8,938.733333333333,926.733333333333,922,937.133333333333,940.6,911.666666666667,924.533333333333,952.066666666667,932.8,921.733333333333,917.466666666667,921.466666666667,955.4,967.933333333333,943.133333333333,896.666666666667,899.6,897.666666666667,917.733333333333,918.133333333333,919.8,923.666666666667,942.733333333333,942.266666666667,935,941.266666666667,951,960.6,949.4,964.8,947.466666666667,957.866666666667,952.466666666667,940.266666666667,927.133333333333,948.2,971,978.733333333333,954.333333333333,939.066666666667,936.866666666667,942.066666666667,941.733333333333,946.666666666667,947.866666666667,927.066666666667,953.266666666667,946,936.066666666667,925.266666666667,957.333333333333,965.333333333333,958.333333333333,942.6,927.866666666667,951.8,963.733333333333,982.8,994.266666666667,982.333333333333,966.2,957.333333333333,959.4,991,986.066666666667,973.733333333333,996.4,1012.2,1023,1049.2,1066.46666666667,1083.2,1085.73333333333,1113.8,1167.6,1241.8,1352.53333333333,1461.4,1597.8,1713,1764.33333333333,1776.86666666667,1744.26666666667,1696.06666666667,1642.73333333333,1609.8,1616.4,1613.2,1589.06666666667,1577.33333333333,1568.2,1553,1528.66666666667,1521.2,1512.8,1485.66666666667,1427.06666666667,1344.06666666667,1306.2,1279.46666666667,1264.6,1260.46666666667,1263.4,1233.6,1203.4,1187,1149.4,1158.46666666667,1150.66666666667,1154.8,1137.4,1109,1114.8,1108.06666666667,1118.4,1117.33333333333,1120.66666666667,1110,1131.66666666667,1149.6,1169.53333333333,1166.73333333333,1163,1130.6,1106.2,1113.66666666667,1119.06666666667,1101.06666666667,1118.06666666667,1137.93333333333,1193.6,1221.13333333333,1264.93333333333,1279.6,1250,1211.06666666667,1179.53333333333,1162.8,1176.13333333333,1172.66666666667,1167.93333333333,1186.73333333333,1193.33333333333,1194.8,1184.06666666667,1167.26666666667,1165.2,1165.86666666667,1183.13333333333,1208.86666666667,1201.53333333333,1196.66666666667,1183.8,1167.26666666667,1156,1121.06666666667,1115.06666666667,1105.4,1087.86666666667,1132.13333333333,1152.8,1152.8,1121.4,1095.86666666667,1107,1113.53333333333,1118.13333333333,1120.73333333333,1117.4,1113.86666666667,1094.46666666667,1102.4,1122.4,1132.6,1157.8,1155.33333333333,1124,1129.13333333333,1144.13333333333,1127.13333333333,1129.73333333333,1175.93333333333,1201.66666666667,1199.4,1196.2,1176.66666666667,1150.33333333333,1154,1170.53333333333,1193,1180.13333333333,1147.46666666667,1094.53333333333,1096.26666666667,1116.6,1147.33333333333,1130.8,1111.6,1094.73333333333,1092.33333333333,1104,1081,1085.53333333333,1098.33333333333,1089,1054.2,1042.4,1073.73333333333,1091.53333333333,1119.4,1081.06666666667,1080.66666666667,1050.66666666667,1073.53333333333,1082.4,1119.53333333333,1132.86666666667,1143.8,1197.13333333333,1206.93333333333,1206.46666666667,1202.6,1206.26666666667,1213.46666666667,1221.06666666667,1248.4,1277.2,1319.66666666667,1378.4,1447.06666666667,1532.33333333333,1636.13333333333,1727.6,1773.26666666667,1772.8,1755.66666666667,1738.46666666667,1709.4,1684.4,1673,1741.4,1806.13333333333,1857.2,1891.53333333333,1902.86666666667,1894.6,1860.13333333333,1809.6,1762.93333333333,1705.2,1656.93333333333,1652.46666666667,1671.46666666667,1720.6,1748.4,1753.06666666667,1750.2,1737.73333333333,1701.4,1636.66666666667,1549.2,1492.13333333333,1454.66666666667,1395.26666666667,1347.6,1320,1307.06666666667,1296.13333333333,1267.46666666667,1241.66666666667,1212.33333333333,1179.66666666667,1159.06666666667,1160.73333333333,1173.13333333333,1184.93333333333,1197.26666666667,1168.6,1135,1129.13333333333,1095.33333333333,1100.53333333333,1101.86666666667,1099.46666666667,1076.46666666667,1070.93333333333,1076.6,1051,1056.06666666667,1074.2,1072.6,1047.6,1066.6,1095.46666666667,1121.6,1135.93333333333,1152.53333333333,1151.86666666667,1152.33333333333,1173.06666666667,1199.33333333333,1219.53333333333,1262.06666666667,1265.93333333333,1252.13333333333,1206.26666666667,1170.13333333333,1123.33333333333,1141.86666666667,1108.4,1080.86666666667,1090.26666666667,1077.66666666667,1056.8,1068.6,1085.86666666667,1081,1059.06666666667,1072.26666666667,1050.8,1069.66666666667,1055,1046,1004.06666666667,984.333333333333,991.066666666667,993.133333333333,999.066666666667,998.4,991,975.533333333333,954.533333333333,932.933333333333,919.533333333333,956.866666666667,961.933333333333,961.8,957,951,956.733333333333,985.866666666667,962,961.666666666667,967.533333333333,978.866666666667,990,989,972.266666666667,978.8,956.6,963.8,975.4,976,962.8,958.333333333333,968.4,942.066666666667,931.133333333333,949.4,936.733333333333,970.466666666667,998.933333333333,971.333333333333,961.733333333333,954.666666666667,952.933333333333,951.2,963.2,975.266666666667,967.666666666667,940,959.066666666667,944.466666666667,968.133333333333,965.066666666667,956,961.533333333333,965.8,952,934.533333333333,945.066666666667,952.466666666667,953.133333333333,956.266666666667,952.866666666667,967.2,955.2,951.133333333333,955.533333333333,959.933333333333,949.266666666667,952.133333333333,948.533333333333,959.733333333333,950.4,960.866666666667,944.533333333333,940.2,945.066666666667,951.133333333333,958.666666666667,943.4,936.466666666667,937.6,941.333333333333,950.533333333333,967.066666666667,962.466666666667,975.2,977.6,959.733333333333,928.466666666667,949.533333333333,975.8,976.133333333333,959.933333333333,948.8,986.133333333333,976.333333333333,945.866666666667,947.266666666667,928.8,925.266666666667,956.466666666667,965.6,959.8,972.4,969.866666666667,948.533333333333,943.133333333333,962.133333333333,968.066666666667,945.8,931.6,940.066666666667,952.066666666667,922.666666666667,926.6,945.466666666667,952.933333333333,964.933333333333,986.533333333333,986.866666666667,977.933333333333,960.466666666667,944.933333333333,925.466666666667,947.666666666667,956.2,957.4,973.466666666667,977.066666666667,967.866666666667,973.4,968.333333333333,973.866666666667,963.266666666667,940.4,943.4,955.6,948.333333333333,932.066666666667,943.933333333333,946.266666666667,946.866666666667,933.733333333333,934.6,958.866666666667,950,940.266666666667,937.2,927.333333333333,931.2,942.4,961.866666666667,938.6,947.6,939.6,949.933333333333,948.666666666667,957.533333333333,952.733333333333,932.933333333333,931.2,938.4,921.066666666667,922.6,923.666666666667,933.2,953.066666666667,970.933333333333,970,973.733333333333,938.4,920.4,939.333333333333,945.133333333333,949.866666666667,963.866666666667,973.8,968.466666666667,950.933333333333,937.133333333333,934.8,924.733333333333,949.066666666667,930.466666666667,920,919.866666666667,925.666666666667,923.266666666667,929.2,947.2,935.4,952.133333333333,965,953.733333333333,940.4,955.133333333333,935.2,936.8,934.4,935.666666666667,935.8,943.6,954.666666666667,927.866666666667,924.733333333333,928.2,919.4,937,933.066666666667,927.333333333333,920.8,922.8,934.6,954.666666666667,976.933333333333,950.6,928.133333333333,930.8,914.4,926.866666666667,912.666666666667,926.466666666667,955.933333333333,957.933333333333,937.866666666667,928,922,946.066666666667,960.333333333333,971.4,938.933333333333,914.533333333333,930.133333333333,919.6,931.866666666667,936.266666666667,919.266666666667,923.466666666667,945.066666666667,955.733333333333,939.666666666667,935.533333333333,932.933333333333,935.266666666667,939.8,940.066666666667,941.666666666667,938.8,969.066666666667,972.2,932.933333333333,923.666666666667,931.266666666667,933.8,923.533333333333,938.333333333333,936.133333333333,913.866666666667,921.066666666667,952.4,951.133333333333,932.666666666667,924.6,959.4,947.333333333333,924.266666666667,906.666666666667,922.533333333333,950,985.266666666667,983.533333333333,964.2,945.6,945.6,947.733333333333,958.8)
setZeros<-function(data,factors=list()){
for(factor in factors){
sFac<-paste('d',factor,sep='')
data[sFac]<-rep(0,length(data[sFac]))
}
return(data)
}
data.dwt<-dwt(data[[2]],n.levels=8)
opar<-par(mfrow=c(4,2),mar=c(2,2,2,2))
mlist<-c(1:8)
for(iFac in 1:8){
#flist<-mlist[mlist!=iFac]
ndata<-setZeros(data.dwt,iFac)
plot(idwt(ndata),type='l')
}