Could you provide us with data to test the code? use dput (and limit the size!!!!!)
eg: dput(vowel_features) dput(mask_features) Without this information, it's impossible to say what's going wrong. It looks like you're doing something wrong in the selection. What should vowel_features[15] return? Did you check it's actually what you want? Did you use str(G) to check the type? Cheers Joris On Thu, May 27, 2010 at 5:28 PM, cobbler_squad <la.f...@gmail.com> wrote: > > Joris, > > You are a life saver. Based on two sample files above, I think lda should > go > something like this: > > vowel_features <- read.table(file = "mappings_for_vowels.txt") > mask_features <- data.frame(as.matrix(read.table(file = > "3dmaskdump_ICA_37_Combined.txt"))) > G <- vowel_features[15] > > cvc_lda <- lda(G~ vowel_features[15], data=mask_features, > na.action="na.omit", CV=TRUE) > > ERROR: Error in model.frame.default(formula = G ~ vowel_features[15], data > = > mask_features, : > invalid type (list) for variable 'G' > > I am clearly doing something wrong declaring G (how should I declare > grouping in R when I need to use one column from vowel_feature file)? Sorry > for stupid questions and thank you for being so helpful! > > --------- > again, sample files that I am working with: > > mappings_for_vowels.txt: > > V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 > V21 V22 V23 V24 V25 V26 > 1 E 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 1 0 > 0 0 0 0 0 0 > 2 o 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0 > 1 0 1 0 0 0 > 3 I 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 0 > 0 0 0 0 0 0 > 4 ^ 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 1 0 > 0 0 0 0 0 0 > 5 @ 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 1 > 0 0 0 0 0 0 > > and the mask_features file is: > > V42 V43 V44 V45 V46 > V47 V48 V49 > [1,] 2.890891625 2.881188521 2.887784444 -2.882606612 -2.888877341 > 2.879834384 2.886483229 2.883815864 > [2,] 2.763404707 2.756198683 2.761863881 -2.756827983 -2.762268531 > 2.754305072 2.760017050 2.758399799 > [3,] 0.556614506 0.556377530 0.556247414 -0.556300910 -0.556098321 > 0.557495060 0.557383073 0.556867424 > [4,] 0.367065248 0.366962036 0.366870087 -0.366794442 -0.366644148 > 0.366613343 0.366537320 0.366953464 > [5,] 0.423692393 0.421835623 0.421741829 -0.421897460 -0.421659824 > 0.421567705 0.421465738 0.422407838 > > -- > View this message in context: > http://r.789695.n4.nabble.com/Linear-Discriminant-Analysis-in-R-tp2231922p2233333.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help@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. > -- Joris Meys Statistical Consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control Coupure Links 653 B-9000 Gent tel : +32 9 264 59 87 joris.m...@ugent.be ------------------------------- Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php [[alternative HTML version deleted]] ______________________________________________ R-help@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.