Displaying 20 results from an estimated 11000 matches similar to: "Pearson's residuals in logistic regression (PR#2539)"
2013 Mar 28
2
hierarchical clustering with pearson's coefficient
Hello,
I want to use pearson's correlation as distance between observations and
then use any centroid based linkage distance (ex. Ward's distance)
When linkage distances are formed as the Lance-Williams recursive
formulation, they just require the initial distance between observations.
See here: http://en.wikipedia.org/wiki/Ward%27s_method
It is said that you have to use euclidean
2012 Jun 01
1
Violation of sample independence in Pearson's product-moment correlation
Hi all:
There was a concern raised by reviewers of a manuscript of mine over the
proper execution of a Pearson's correlation. In brief, this was undertaken
in order to determine the relationship between the extent of wheel running
(y axis) and ethanol intake (x axis) across three, separate 10 day periods
in 7 animals.
In the paper, the correlational plots for each 10 day-period had 70 data
2009 Jun 03
1
Validity of Pearson's Chi-Square for Large Tables
Is Pearson's Chi-Square test for contingency tables asymptotically unbiased
for large tables (large degrees of freedom) regardless of the expected
values in each cell? The rule of thumb is that Pearson's Chi-square should
not be used when large numbers of cells have expected values < 5. However,
I compared the results on 4x4 contingency tables for R's chisq.test using
chi-square
2003 Apr 11
1
Pearson's Chi-squared Test
How i can perform a Pearson's Chi-squared Test in this data set:
| Outcome
-----------------+-----------+----------------------------------+
Treatment | Sex | None |Some | Marked | Total
-----------------+------------+--------+--------+-------------+
Active | Female | 6 | 5 | 16 | 27
2001 Dec 19
1
Pearson residuals in quasi family
Hi all,
This is a very silly question or something escapes me:
Let obj a simple gam poisson model. Let
>obj<-gam(....,family=poisson)
>obj1<-update(obj, family=quasi(link="log", var="mu"))
>From summary.glm(obj1) the dispersion parameter is estimated 1.165; In fact
it is:
> (predict(obj1, se.fit=T)$se.fit[1:5]/predict(obj, se.fit=T)$se.fit[1:5])^2
4
2001 Nov 16
2
pearson residuals in glm for binomial response (PR#1175)
R version 1.3.0
OS: SunOS 5.7, but I think the same problem occurs with Windows
An incorrect formula seems to be used to calculate the pearson residuals
for a generalized linear model with a binomial response. Here is a
simple program which gives (a) the pearson residuals calculated directly,
(b) the pearson residuals from glm, and (c) the deviance residuals from
glm. The first and last
2007 Oct 17
1
Documentation for Pearson Residuals
Greetings,
I have been using lm to perform weighted linear regressions and resid to
extract the residuals. I happened upon some class notes on the internet that
described how one can specify type="pearson" in resid to extract the
weighted residuals. Where is this option documented? And if you know that,
what is the best way to find such documentation if you don't know that the
2013 Feb 14
0
pearson's correlation and cross-correlation issue
Hello,
I want to compute the pearson's correlation, but even for signals that are
shifted.
For example having two signals like:
1 1 2 1 1
and
1 2 1 1 1
the correlation is very low.. but if we shift them in the right we get much
better correlation.
I know that cross-correlation is used to find the best offset (where
correlation will be bigger).
Is there any metric that can do this job all
2002 Jan 13
0
weighted regression: Osius & Rojek's test for logistic regression models
Dear all,
I am trying to implement goodness-of-fit tests for logistic regression
models described in : Hosmer, D.W., Lemeshow, S., 2000. Applied logistic
regression. New-York, John Wiley & Sons, Inc., 373 p. (pp. 152-154)
Namely I would like to reproduce Osius & Rojek's test (Osius, G., Rojek,
D., 1992. Normal goodness-of-fit tests for multinomial models with large
degrees of
2011 Mar 16
1
Standardized Pearson residuals (and score tests)
Hi Peter and others,
If it helps, I wrote a small function glm.scoretest() for the statmod
package on CRAN to compute score tests from glm fits. The score test for
adding a covariate, or any set of covariates, can be extracted very neatly
from the standard glm output, although you probably already know that.
Regards
Gordon
---------------------------------------------
Professor Gordon K
2008 Nov 27
2
1-Pearson's R Distance
Hi again List,
Well this time I’m writing for a friend (really J). He needs to create a
distance matrix based on an abundance matrix using the 1-Pearson’s R index.
Well I told him to look at the proxy package, but there is only Pearson
Index. He needs it to perform a clustering. Well, as soon as he told me
there proxy only had the Pearson index I thought: “He could just do
something like
2007 Nov 03
1
Pearson residuals
Dear Sirs
What is the best aproximation to the standardized normal distribution:
necessidade = c("sem necessidade","com necessidade")
tipo =c("CE-1", "CE-2", "CE-3")
dados=c(20,34,44,69,9,3)
Tabela =cbind(expand.grid(list(Necessidade=necessidade, Tipo=tipo)),
count=dados)
Tabela.array=tapply(Tabela$count, Tabela[,1:2], sum)
ni =
2010 May 05
1
rcorr p-values for pearson's correlation coefficients
Hi! All,
To find co-expressed genes from a expression matrix of dimension (9275
X 569), I used rcorr function from library(Hmisc) to calculate pearson
correlation coefficient (PCC) and their corresponding p-values. From
the correlation matrix (9275 X 9275) and pvalue matrix (9275 X 9275)
obtained using rcorr function, I wanted to select those pairs whose
PCC's are above 0.8 cut-off and then
2012 Feb 08
0
glm.fit and pearson's correlation coefficient
I did a linear correlation of data using glm.fit and stored the output in the
object "f":
f <- glm.fit(x, y, w)
I am intereseted in estimating the quality of the correlation. I am used to
do it using pearson correlation coefficient "r" or "r^2". Can I extract this
coefficient from the output of glm.fit?
Is there another number in the output of glm.fit that
2008 Apr 05
2
pearson's correlation
Hello,
I used the function cor to calculate the pearson correlation coefficient between variables. However, the resulting values do not correspond to the outcome of my excel-calculations, for which I used the formula Cor(x,y)=Cov(x,y)/(SD(x)*SD(y))
So my question is: How does the function "cor" compute the pearson correlation coefficient?
Thank you in advance,
Ake Nauta
2009 Mar 21
1
Goodness of fit for negative binomial model
Dear r list,
I am using glm.nb in the MASS package to fit negative binomial models to data on manta ray abundance, and AICctab in the bbmle package to compare model IC. However, I need to test for the goodness of fit of the full model, and have not been able to find a Pearson's Chi Squared statistic in any of the output. Am I missing it somewhere? Is there a way to run the test using
2010 Jan 05
1
bootstrapping a matrix and calculating Pearson's correlation coefficient
Hi All,
I have got matrix 'data' of dimension 22000x600. I want to make 50
independent samples of dimension 22000x300 from the original matrix 'data'.
And then want to calculate pearsons CC for each of the obtained 50 matrices.
It seems it is possible to do this using 'boot' function from library boot
but I am not able to figure out how? I am really stuck. Please help!
2012 May 29
2
a question about "by" and "ddply"
Hi all,
I have a data set (df, n=10 for the sake of simplicity here) where I have two continuous variables (age and weight) and I also have a grouping variable (group, with two levels). I want to run correlations for each group separately (kind of similar to "split file" in SPSS). I've been experimenting with different functions, and I was able to do this correctly using ddply
2010 Mar 11
2
logistic model diagnostics residuals.lrm {design}, residuals()
I am interested in a model diagnostic for logistic regression which is normally distributed (much like the residuals in linear regression with are ~ N(0,variance unknown).
My understanding is that most (all?) of the residuals returned by residuals.lrm {design} either don't have a well defined distribution or are distributed as Chi-Square.
Have I overlooked a residual measure or would it be
2002 Oct 21
2
More Logistic Regression Tools?
I've been using R to do logistic regresssion, and that's working
well, but there are two things I haven't figured out how to do.
(1) Is there some pre-existing function that will let you compute the
odds ratios and confidence intervals for them for a specific fit. I know
how to do this manually or even write a function that I can call with
the coefficients and se, but