search for: hongsheng

...In these two year, R shows the error called 'singular gradient' in the model. Could anyone please help me to fix the error? I don't think there are any problems with my data because I can fit the same model to 2004 data without any problems using SAS. Thank you very much in advance. Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571 Fax:757.683.5293

Attached is a histogram chart with a break on y-axis which my friend made for me using phot shop. Do anyone know how to make such a break using R? Thanks in advance. (See attached file: 2007 age distribution.doc) Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571 Fax:757.683.5293

My stat textbook tells me that using Shapiro-Wilk test for each variable one by one is not equal to a test for multivariate normality as a whole. Does R have a function of testing for multivariate normality? Thanks. Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571 Fax:757.683.5293

Where can I download Bhat package? I need to use plkhci function in this package, which estimates confidence interval using likelihood ratio. Thanks. Hongsheng (Hank) Liao, PhD. Lab Manager Center for Quantitative Fisheries Ecology Old Dominion University 757-683-4571 [[alternative HTML version deleted]]

...When I print out princomp's loading outputs, there is alwasy a section for "SS loading", "Proportional Var" and "Cumulative Var". Anybody can tell what they are for? Or anyone can direct me to some reference to read about? Any help will be highly appricated. Hongsheng [[alternative HTML version deleted]]

...#39;t be converged, such as 2000, the loop stops in that year (2000) even though the rest years (2001 to 2008) can be converged. I would really appreciate if someone helps me on how to force the loop function to finish all 10 rounds even when errors occur due to no convergences for some of years. Hongsheng (Hank) Liao Lab Manager Center for Quantitative Fisheries Ecology Old Dominion University Phone: 757-683-4571 [[alternative HTML version deleted]]

Hi all, I have a huge data set. All I want to do is to change the column names of the it. So if I use read.table, I can read the data in and change colnames and then write back, such as, t <- read.table("a.txt", header = T, sep = "\t") colnames(t) colnames(t) <- c("....) # new column names write.table(t, "a.txt", quote = F, sep = "\t",