[R] Problems using hetcor (polycor)

Mark Difford mark_difford at yahoo.co.uk
Thu Aug 7 13:23:05 CEST 2008


Hi Birgitle,

You need to get this right if someone is going to spend their time helping
you. Your code doesn't work: You have specified more columns in colClasses
than you have in the provided data set.

> TestPart<-read.table("TestPart.txt", header=TRUE,row.names=1,
> na.strings="NA" ,colClasses = Classe81)
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, 
: 
  line 1 did not have 82 elements

> TestPart<-read.table("TestPart.txt", header=TRUE,row.names=1,
> na.strings="NA")
Warning message:
In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings,  :
  number of items read is not a multiple of the number of columns

> length(names(TestPart))
[1] 72

> length(Classe81)
[1] 82

That will never work.

HTH, Mark.


Birgitle wrote:
> 
> Sorry if this post should be long but I tried to give you a piece of my
> data to reproduce my error message using hetcor:
> 
> Fehler in result$rho : $ operator is invalid for atomic vectors
> Zusätzlich: Warning messages:
> 1: In polychor(x, y, ML = ML, std.err = std.err) :
>   1 row with zero marginal removed
> 2: In polychor(x, y, ML = ML, std.err = std.err) :
>   the table has fewer than 2 rows
> 
> 
> Error in result$rho : $ operator is invalid for atomic vectors
> Additional: Warning message:
> 1: In polychor(x, y, ML = ML, std.err = std.err) :
>   1 row with zero marginal removed
> 2: In polychor(x, y, ML = ML, std.err = std.err) :
>   the table has fewer than 2 rows
> 
> Use tab delimited data at the end of the post.
> Copy in Texteditor and save as TestPart.txt 
> 
> Then use the following code
> 
> library(methods)
> setClass("of")
> setAs("character", "of", function(from) as.ordered(from)) 
> 
> Classe81<-cclasses <- c(rep("factor", 64),  rep("numeric",6), rep ("of",
> 12))
> 
> TestPart<-read.table("TestPart.txt", header=TRUE,row.names=1,
> na.strings="NA" ,colClasses = Classe81)
> 
> str(TestPart)
> 
> library(polycor)
> 
> TestPart.hetcor<-hetcor(TestPart, use="complete.obs")
> 
> this will produce the above mentioned error message that I can not
> interprete.
> 
> Would be great if somebody could help me to solve the problem.
> 
> Thanks
> 
> B.
> 
> 
> 
> 
> 	1	2	3	4	7	8	9	10	12	13	14	15	16	17	18	19	21	22	23	25	27	28	29	30	31	33	34
> 35	36	37	38	39	40	41	42	43	44	45	46	47	48	49	50	51	52	53	54	55	56	58	59	60
> 61	62	63	64	65	66	67	68	69	70	71	72	73	74	75	76	77	78	79	80
> AX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	1	0	0	1	1	0	0	0	0	0	1	0	0	0	1	0	0	0	0
> 0	1	0	0	0	0	0	0	0	0	0	1	1	0	0	1	0	1	25	5	9	1	8.5	2.5	3	5	2	2	3	3	1	1	1	2	1
> 2
> BX	1	1	0	0	1	0	0	1	NA	NA	NA	0	0	0	0	1	0	0	1	0	0	1	0	NA	NA	NA	NA	NA	NA	NA
> NA	0	0	0	1	0	NA	NA	NA	NA	NA	NA	NA	0	0	0	0	1	1	0	0	0	1	1	NA	NA	6	1	3.25
> 2.25	5	5	2	2	3	3	1	1	1	1	1	1
> CX	1	1	0	0	1	0	0	1	1	0	0	0	1	0	1	0	0	1	0	1	0	0	0	0	1	1	0	0	0	0	0	0	1	0	0	0
> 0	1	0	0	0	0	0	0	0	0	1	0	0	0	0	1	0	0	15	3.5	6	1	5.5	5.5	5	5	2	2	1	2	1	1	1	1
> 2	2
> DX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	0	1	0	0	1	0	0	1	0	1	0	0	0	1	0	0	0	0	1
> 1	0	0	0	0	0	0	0	0	0	1	0	1	0	1	0	1	0	50	17.5	7.5	2.5	8.5	5	5	5	2	2	2	3	1	1
> 1	1	3	3
> EX	1	0	1	0	1	0	0	1	NA	NA	NA	0	0	0	1	1	0	1	1	0	1	0	0	0	0	0	0	1	0	0	0	0	1	0
> 0	0	0	0	0	0	1	0	0	0	0	0	1	0	1	0	0	0	1	0	NA	NA	14.5	30	13	2.5	3	3	1	1	4	4	1
> 1	1	1	1	1
> FX	1	0	1	0	1	0	0	1	0	0	0	0	0	0	0	1	0	0	1	0	1	0	0	0	1	1	1	0	0	1	1	0	0	0	1	1
> 1	0	0	0	1	0	0	0	0	0	1	0	1	0	1	1	0	0	165	25	11.5	15	12	6.5	5	5	1	1	3	3	1	1
> 1	1	4	5
> GX	1	0	1	0	1	0	0	1	0	0	1	0	0	0	0	1	1	1	0	0	1	0	0	0	0	0	1	0	0	1	1	0	1	0	0	1
> 0	0	0	0	1	0	0	0	1	0	0	0	0	0	1	0	1	0	40	20	14.5	9.5	11	10	3	3	1	1	1	3	1	1	3
> 4	1	3
> HX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	1	0	0	1	0	0	NA	NA	NA	NA	NA	NA	NA	NA
> NA	NA	NA	NA	NA	1	0	0	0	1	0	0	0	0	0	0	0	1	1	1	0	1	0	70	18.5	1	12.5	7.2	7.2
> 5	5	2	2	3	3	3	3	1	1	2	4
> JX	1	0	1	0	1	0	0	1	0	0	0	0	0	1	1	0	0	1	0	1	0	0	0	0	1	0	1	0	0	0	0	0	1	1	0	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	1	0	65	9	11	18	4.5	3.75	1	3	2	3	2	2	1	1	1
> 1	2	2
> KX	1	0	1	0	1	0	0	1	NA	NA	NA	0	0	1	0	0	0	0	1	0	1	0	0	0	0	0	1	1	0	0	1	0	0	1
> 1	0	0	1	0	0	0	0	0	0	0	0	0	1	0	0	0	1	0	0	NA	NA	4.25	4.5	2.1	1.6	3	3	2	2	2	3
> 1	1	1	2	1	1
> LX	1	0	0	1	1	0	0	1	0	1	0	0	0	0	1	0	0	0	1	1	0	0	0	0	0	0	0	1	1	0	0	1	1	0	0	0
> 0	1	0	0	0	0	0	0	1	0	0	1	0	0	0	0	1	0	80	32.5	13	4	8.5	5.25	2	3	3	3	3	3	1	1
> 1	1	3	4
> MX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	1	1	0	1	0	0	1	0	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	1	1	0	0	0	0	0	0	0	0	0	1	45	7.5	14	1	11	5.5	3	4	2	2	3	3	1	1	1	1
> 2	3
> NX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	1	0	0	1	0	0	1	1	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	1	42.5	12.5	17.5	1	11	NA	5	5	2	2	3	3	1	1
> 1	1	1	2
> OX	1	1	0	0	1	0	0	1	NA	NA	NA	0	0	0	0	1	0	0	1	0	1	0	1	0	1	1	0	0	0	0	0	0	0	0
> 1	0	0	0	0	0	1	0	0	0	1	0	0	1	0	0	0	0	0	1	NA	NA	25	1	25	NA	3	4	1	2	3	3	1	1	3
> 3	1	1
> PX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	1	1	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	42.5	11.5	21.5	2	25	NA	4	5	2	2	3	3	1	1
> 1	2	1	2
> QX	1	1	0	0	1	0	0	1	NA	NA	NA	0	0	0	1	1	0	0	1	1	0	0	0	0	1	1	1	0	0	0	0	0	0	0
> 1	0	0	1	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	NA	NA	27.5	1.5	19	NA	3	3	2	2	3	3	1
> 1	1	2	1	1
> RX	1	1	0	0	1	0	0	1	NA	NA	NA	0	0	0	0	1	0	1	0	0	0	1	0	0	0	1	0	0	0	0	0	0	0	0
> 1	0	0	0	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	NA	NA	30	1	24	NA	3	4	2	2	3	3	1	1	3
> 3	1	1
> SX	1	1	0	0	1	0	0	1	0	0	1	0	0	0	1	1	0	1	0	1	0	1	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	1	47.5	20	22.5	2	15	NA	5	5	2	2	3	3	1	1	1
> 1	2	2
> TX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	1	0	0	1	0	1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	1	40	11	15	1.5	11.5	NA	5	5	2	3	3	3	1	1	1
> 1	2	2
> UX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	0	1	1	0	0	0	0	1	1	0	0	0	0	0	0	0	0	1	0
> 0	1	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	65	35	39	2.5	20	NA	3	4	2	2	3	3	1	1	2	2
> 1	2
> VX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	1	1	0	0	1	1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	85	10.5	19.5	1.5	13.5	NA	5	5	2	3	3	3	1
> 1	2	2	2	2
> WX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	1	1	1	0	0	0	0	0	1	1	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	1	155	18	30	2	17.5	NA	5	5	2	3	3	3	1	1	1
> 2	2	3
> XX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	1	90	10.5	25	1	19	NA	5	5	2	2	3	3	1	1	1	1
> 2	3
> YX	1	1	0	0	1	0	0	1	NA	NA	NA	0	0	0	0	1	0	0	1	1	1	1	0	0	0	1	0	1	0	0	0	0	0	0
> 1	0	0	1	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	1	NA	NA	30	2.5	20.25	NA	3	3	2	2	3	3
> 1	1	2	2	1	1
> ZX	1	1	0	0	1	0	0	1	0	0	0	0	0	0	0	1	0	0	1	0	0	1	0	0	0	0	0	1	0	0	0	1	1	0	0	0
> 0	0	0	0	1	0	0	0	0	0	1	0	0	0	0	0	0	1	105	25	37.5	1	37.5	2.25	5	5	2	3	3	3	1
> 1	1	1	2	2
> AY	0	1	1	0	1	1	1	0	1	0	0	0	0	0	0	1	0	1	1	0	0	1	0	1	0	1	0	0	0	0	0	0	0	0	1	1
> 0	1	0	0	0	0	0	0	0	0	0	1	1	0	0	1	0	0	85	6.5	11.5	3	11.5	3.25	5	5	2	2	3	4	1
> 1	1	2	2	4
> BY	0	0	1	0	1	0	1	0	0	1	0	0	0	0	0	1	0	1	0	1	0	1	0	0	1	0	1	0	0	0	0	0	1	1	0	0
> 0	1	0	0	1	0	0	0	0	0	1	1	0	0	0	0	1	0	38	10	6.5	8	3	2.5	4	5	2	2	2	3	1	1	1	2
> 2	3
> CY	0	1	0	0	0	1	1	0	1	0	0	0	1	1	1	1	0	1	0	0	1	0	0	1	1	0	0	0	0	0	1	0	1	1	1	1
> 1	1	0	0	0	0	0	0	0	0	1	0	1	0	0	1	0	0	15	3.5	5	1.5	3.25	3.5	5	5	2	3	2	2	1	1
> 1	1	2	2
> DY	0	1	0	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	0	0	1	0	NA	NA	NA	NA	NA	NA	NA	NA
> NA	NA	NA	NA	NA	1	0	0	0	0	0	0	0	0	1	0	0	1	0	0	1	0	0	70	20	NA	NA	5.5	3.5	5	5
> 3	3	3	3	3	3	1	1	9	9
> EY	0	0	0	1	1	0	0	1	1	0	0	0	0	0	0	1	0	1	0	0	1	0	0	0	1	0	0	0	0	0	0	0	1	0	0	0
> 0	0	0	0	1	0	0	0	0	0	0	0	1	0	0	0	1	0	20	7	8	25	5.5	3.25	3	3	3	3	3	3	1	1	1	1
> 2	2
> FY	0	0	0	1	1	0	0	1	0	1	0	0	0	0	0	1	0	0	1	0	1	0	0	0	1	0	1	0	0	0	0	0	0	1	1	1
> 1	0	0	0	1	0	0	0	0	0	1	0	1	0	0	0	1	0	150	25	13	24	11.5	4.5	5	5	3	3	4	4	1	1
> 1	1	5	6
> GY	0	0	0	1	1	0	0	1	1	0	0	0	0	0	0	1	0	0	1	0	1	0	0	0	1	0	1	0	0	0	0	0	1	0	0	0
> 0	0	0	0	1	0	0	0	1	0	0	0	0	0	0	0	1	0	50	15	12.5	17.5	8.5	5.75	3	3	3	3	2	3	1
> 1	1	1	2	3
> HY	0	1	0	0	0	1	0	1	0	1	0	0	0	0	0	1	0	1	0	0	1	0	0	0	1	0	0	0	0	0	0	0	0	1	0	1
> 1	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	1	0	110	15	9.5	7.5	11	3.35	5	5	3	3	4	4	3	3
> 1	1	4	5
> JY	0	0	0	1	1	0	0	1	1	1	0	0	1	1	0	0	0	1	0	1	0	0	0	0	0	1	0	1	0	0	0	0	1	0	0	0
> 0	1	0	0	0	1	0	0	0	0	0	1	0	0	0	0	1	0	60	11.5	10	30	3.25	3.75	3	3	2	2	1	1	1
> 1	1	1	2	3
> KY	0	0	0	1	1	0	0	1	NA	NA	NA	1	0	1	0	0	1	0	1	0	1	1	0	1	1	0	1	0	0	0	0	0	0	0
> 1	0	0	1	1	0	0	0	0	0	0	0	0	1	0	0	0	1	0	1	NA	NA	6.5	23	2.25	1	1	2	2	2	3	3	1
> 1	2	2	2	2
> LY	0	0	1	0	1	0	0	1	0	1	0	0	0	0	1	0	0	1	0	1	1	0	0	0	1	0	1	0	0	0	0	0	1	0	0	0
> 0	1	0	0	0	0	0	0	0	0	1	1	0	0	0	1	0	0	52.5	20	11	5	6.25	4.25	2	2	2	2	3	3	1	1
> 1	1	2	4
> MY	0	1	1	0	0	1	0	1	0	1	0	0	1	1	0	0	0	0	1	0	1	0	0	0	0	0	1	1	0	0	0	0	1	0	0	0
> 0	0	1	0	0	1	0	0	0	0	1	0	0	0	0	1	0	0	60	15	3	15	1.5	1.25	5	5	2	2	2	2	1	1	1
> 1	6	6
> NY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	1	1	0	1	1	0	1	1	0	0	0	1	1	0	0	1	0	0	1	0	0	0
> 0	0	0	0	1	0	0	0	0	0	1	0	0	0	0	1	0	0	50	14	3.75	25	2.85	2.4	5	5	2	2	1	2	1	1
> 1	1	5	6
> OY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	0	0	1	0	0	0	0	1	0	0	1	0	0	1	0	1	0
> 0	0	0	0	1	0	0	0	0	0	0	1	0	0	0	1	0	0	45	14	4	15	2.75	1.75	5	5	2	2	2	3	1	1	2
> 3	5	5
> PY	0	0	1	1	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	0	0	1	0	0	0	0	1	0	0	1	0	0	1	1	1	0
> 0	0	0	0	1	0	0	0	0	0	0	1	0	0	0	1	0	1	60	18.5	3.75	15	2.75	2.1	5	5	2	2	2	3	1
> 1	1	2	4	6
> QY	0	0	1	0	0	1	0	1	0	1	0	0	0	1	1	1	0	0	1	1	0	1	0	0	1	1	1	0	0	0	0	0	0	1	1	0
> 0	0	0	0	1	0	0	0	0	0	1	1	0	0	0	0	0	1	290	90	8.5	22.5	2.75	2.25	5	5	2	2	2	3
> 1	1	1	2	7	7
> RY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	0	0	1	0	0	0	0	1	0	0	1	0	0	1	1	1	0
> 0	0	0	0	1	0	0	0	0	0	0	1	0	0	0	0	1	0	55	15.5	7.5	14	7	2.25	5	5	2	2	3	3	1	1
> 2	3	2	4
> SY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	1	0	0	0	0	0	0	0	0	0	1	0	0	1	0	0	0
> 0	0	0	0	1	0	0	0	0	0	1	1	0	0	0	0	1	0	55	20	5.5	30	4.5	2.5	5	5	2	2	3	3	1	1	1
> 1	4	5
> TY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	0	1	1	0	0	0	0	1	0	0	1	0	0	1	1	0	0
> 0	0	0	0	1	0	0	0	0	0	1	0	0	0	0	1	0	0	40	15	4.5	25	2.75	1.75	5	5	2	2	3	3	1	1
> 1	1	4	6
> UY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	1	1	0	0	1	1	0	0	0	0	0	1	1	0	0	1	0	0	1	0	1	0
> 0	1	0	1	1	0	0	0	0	0	1	1	0	0	0	0	1	0	200	27.5	6.5	25	2.45	2.25	5	5	2	2	2	3
> 1	1	1	2	7	7
> VY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	1	1	0	0	1	1	0	1	0	0	1	0	1	0	0	1	0	0	1	1	1	0
> 0	0	0	0	1	0	0	0	0	0	0	1	1	0	0	0	1	0	135	17.5	5.5	17.5	2.75	2.25	5	5	2	2	2
> 3	1	1	1	2	5	7
> WY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	0	1	0	0	1	1	0	0	0	0	1	0	1	0	0	1	1	0	1	0	1	1
> 0	0	0	0	1	0	0	0	0	0	0	1	0	0	0	1	0	1	240	35	5.5	14	3.25	2.75	5	5	2	2	2	3	1
> 1	1	2	7	7
> XY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	1	1	0	0	1	1	0	0	0	0	0	1	0	0	0	1	0	0	1	0	0	0
> 0	0	0	1	1	0	0	0	0	0	0	1	0	0	0	1	0	0	175	27.5	4	10.5	2.5	1.875	5	5	2	2	2	2
> 1	1	1	1	6	7
> YY	0	0	1	0	0	1	0	1	0	1	0	0	0	0	1	1	0	0	1	1	0	1	0	0	1	1	1	0	0	0	0	0	0	0	1	0
> 0	0	0	0	1	0	0	0	0	0	1	1	0	0	0	0	0	1	230	42.5	5.6	25	2.5	2.5	5	5	2	2	2	2	1
> 1	1	2	5	7
> ZY	0	1	0	0	0	1	0	1	0	1	0	0	0	0	1	0	0	1	0	0	0	1	0	0	0	0	0	0	0	0	1	0	0	0	0	1
> 1	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	1	115	35	10.5	20	7.5	6.5	5	5	1	1	3	3	1	1
> 1
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 	
> 

-- 
View this message in context: http://www.nabble.com/Problems-using-hetcor-%28polycor%29-tp18867343p18868284.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list