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");