search for: youngjin

Displaying 12 results from an estimated 12 matches for "youngjin".

2006 Mar 06
2
[Q] BIC as a goodness-of-fit stat
Dear R-List I have a question about how to interpret BIC as a goodness-of-fit statistic. I was trying to use "EMclust" and other "mclust" library and found that BIC was used as a goodness-of-fit statistic. Although I know that smaller BIC indicates a better fit, it is not clear to me how good a fit is by reading a BIC number. Is there a standard way of interpreting a BIC
2006 Oct 19
1
[Q] How to fit data to a straight line
Dear R users I have a question about how to fit data to a straight line. I tried nls to do it, but it didn't work. The reason I want to fit data to a straight line is that I need to compare AIC or BIC values of the two models (a straight line model vs a nonlinear curve model). Fitting data to a nonlinear curve is straightforward, but I could not figure out how to fit the data to a straight
2006 Mar 24
1
Multiple error bar plots
...Yn2 // End of the second set of measurements X1, Y13, error of Y13, X2, Y23, error of Y23, ... Xn, Yn3, error of Yn3 // end of the third set of measurements I need to put one errbar plot that include all three cases with different symbols for each data set. How can I do this? Thanks in advance. Youngjin [[alternative HTML version deleted]]
2006 Apr 07
1
[Q] Format of a plot axis label
Dear R-users Can anyone please tell me how to format a label of a plot? I found that R uses a floating point number (e.g. 2.0) sometimes and an integer number (e.g. 2) other times in lables of a plot even though no floating number is needed to show the lables. How can I get R to use only integer numbers in the plot labels? Thanks in advance. Young-Jin [[alternative HTML version deleted]]
2006 May 30
0
weighted correlation coefficient
Dear R-listers Is there a R-package that I can use to compute a weighted correlation coefficient along with its p-values? I want to be able to assign weights to the data points before computing a Pearson correlation coefficient. Thanks in advance. Young-Jin [[alternative HTML version deleted]]
2010 Aug 07
1
[Q] a dummy variable used with sapply
Hello, While learning how to manipulate data with R, I found one example that I could not understand. In the following example, the function definition of "maxcor" has an argument named "i" and I don't understand why. Could someone explain why the maxcor function definition needs to have this argument? maxcor = function(i, n = 10, m = 5) { mat =
2010 Sep 01
1
[Q] Goodness-of-fit test of a logistic regression model using rms package
Hello, I was looking for a way to evaluate the goodness-of-fit of a logistic regression model. After googling, I found that I could use "resid(fit, 'gof')" method implemented in the rms package. However, since I am not used to the "le Cessie-van Houwelingen normal test statistic," I do not know which statistic from the returned from the "resid(fit,
2010 Sep 11
0
[Q] How to extract cross validation results from e1071's svm model
Dear all, Is it possible to extract cross-validation results from e1071's svm model? For example, the following R code shows the result from the 10 fold cross-validation. model = svm(spam ~ ., data = spam, cross = 10) summary(model) But, I could not figure out how to get to the accuracy values from the cross-validation. I looked at the svm method, but did not find any return values. Any
2006 Mar 28
2
[Q] How to make a multi-line title with expression()
Dear R-lister Could anyone know how to make a multi-line title for a plot with expression()? In my plot, the title should be writeen in two lines (because it is two long for one line) and it should use a mathematical expression. I tried to use "\n", but "\n" is ignored in expression() call: hist(diffChangeRequestHintsBeforeAnswering[,4], br = 50, xlab = "Skill Change in
2006 Apr 12
3
[Q] Bayeisan Network with the "deal" package
Dear R-users I am looking for a help in using the "deal" package. I followed the manual and paper from the author's web site to learn it, as shown below, but I could not figure out how to access the local and posterior probability of the nodes in the constructed network. library(deal) data(ksl) ksl.nw <- network(ksl) ksl.prior <- jointprior(ksl.nw) banlist <-
2007 Mar 09
4
How to create a list that grows automatically
Dear R users I would like to know if there is a way to create a list or an array (or anything) which grows automatically as more elements are put into it. What I want to find is something equivalent to an ArrayList object of Java language. In Java, I can do the following thing: // Java code ArrayList myArray = new ArrayList(); myArray.add("object1"); myArray.add("object2");
2007 Mar 09
4
How to create a list that grows automatically
Dear R users I would like to know if there is a way to create a list or an array (or anything) which grows automatically as more elements are put into it. What I want to find is something equivalent to an ArrayList object of Java language. In Java, I can do the following thing: // Java code ArrayList myArray = new ArrayList(); myArray.add("object1"); myArray.add("object2");