Displaying 20 results from an estimated 700 matches similar to: "Re: Polychoric correlations"
2004 Dec 09
1
Re: Tetrachoric and polychoric correlations, Polycor package
A bit late, but you might like to look at
http://www.qimr.edu.au/davidD/polyr.R
Regarding the original posters queries:
You can analyse polychoric correlations as if they were Pearson
correlations using standard software (eg sem), and this usually doesn't do
too badly, or go to AWLS (Browne) in LISREL etc, or ML analysis
of the full multidimensional contingency table using programs such as
2011 May 11
3
Vermunt's LEM in R
I don't know of any R package that can match all the functionality of LEM
eg fitting equality constraints to model parameters a la LISREL.
WRT dumping tables, I would have thought that as.data.frame.table does
pretty much what you want, [not tested]
newtab <- as.data.frame(table(a,b,c))
cat("dim\n")
for(i in seq(1, ncol(newtab)-1) {
cat(nlevels(newtab[,1]," ")
}
2007 Aug 30
0
R-help Digest, Vol 54, Issue 30
Ron Crump wrote:
> Hi,
>
> I have a dataframe that contains pedigree information;
> that is individual, sire and dam identities as separate
> columns. It also has date of birth.
>
> These identifiers are not numeric, or not sequential.
>
> Obviously, an identifier can appear in one or two columns,
> depending on whether it was a parent or not. These should
> be
2000 Jul 17
2
sshd -i problem under 2.1.1p4 as well
I have encountered the same problem with sshd -i (under Mandrake linux 6.1)
as that described by Ben L Perkins, this time with 2.1.1p4: ...
Last login: Mon Jul 17 12:04:50 2000 from orpheus.qimr.edu.au
-bash: ?oe90: command not found
-bash: glorious: command not found
-bash: ?]r90: command not found
-bash: marshall4: command not found
Connection to orpheus.qimr.edu.au closed.
log:
Jul 17
2013 Oct 09
1
Version of L-BFGS-B used in optim etc
Hi.
I just noticed the paper by Morales and Nocedal
Remark on "Algorithm 778: L-BFGS-B: Fortran Subroutines for Large-Scale
Bound Constrained Optimization". TOMS 2011; 38(1): 7
http://www.ece.northwestern.edu/~morales/PSfiles/acm-remark.pdf
which describes a couple of improvements (speed and accuracy) to the
original Netlib code which AFAICT is that still used by optim()
via f2c.
2006 Feb 22
1
var-covar matrices comparison
> Date: Mon, 20 Feb 2006 16:43:55 -0600
> From: Aldi Kraja <aldi at wustl.edu>
>
> Hi,
> Using package gclus in R, I have created some graphs that show the
> trends within subgroups of data and correlations among 9 variables (v1-v9).
> Being interested for more details on these data I have produced also the
> var-covar matrices.
> Question: From a pair of two
2008 May 21
1
table(factor(x), exclude=NULL) (PR#11494)
Hi. I don't know if this a bug or just annoying to me:
> x <- c(1,2,3,NA)
> table(x, exclude=NULL)
x
1 2 3 <NA>
1 1 1 1
> table(factor(x), exclude=NULL)
1 2 3
1 1 1
I don't think many people use factor(x, exclude=NULL): it is not the
default handling of character data by read.table().
Cheers, David Duffy.
--
| David Duffy (MBBS PhD)
2005 Jun 29
2
How to convert "c:\a\b" to "c:/a/b"
I couldn't resist adding a more literal answer
unback <- function(x) {
chars <- unlist(strsplit(deparse(x),""))
chars <- chars[-c(1,length(chars))]
paste(gsub("\\\\","/",chars),collapse="")
}
unback("\n")
| David Duffy (MBBS PhD) ,-_|\
| email: davidD at qimr.edu.au ph:
2005 Jun 29
2
How to convert "c:\a\b" to "c:/a/b"
I couldn't resist adding a more literal answer
unback <- function(x) {
chars <- unlist(strsplit(deparse(x),""))
chars <- chars[-c(1,length(chars))]
paste(gsub("\\\\","/",chars),collapse="")
}
unback("\n")
| David Duffy (MBBS PhD) ,-_|\
| email: davidD at qimr.edu.au ph:
2009 Jan 16
0
R-help Digest, Vol 71, Issue 16
On Fri, 16 Jan 2009, r-help-request at r-project.org wrote:
> Date: Thu, 15 Jan 2009 13:29:03 +0100
> From: Pablo G Goicoechea <pgoikoetxea at neiker.net>
> Subject: Re: [R] How to create a chromosome location map by locus ID
> To: Sake <tlep.nav.ekas at hccnet.nl>
> Cc: r-help at r-project.org
> Message-ID: <496F2C0F.3040304 at neiker.net>
> Content-Type:
2003 Nov 11
0
Re: R-help Digest, Vol 9, Issue 11
Tony Rossini wrote:
>
> For various reasons, I spent part of my time today looking at sample
> size and power calculation tools (don't ask, don't tell...). This
> seems to be one area that R is incredibly weak in (well, nearly all
> stat packages, except perhaps specialized tools and SAS); sure, there
> are a number of functions in various packages:
>
>
2004 Mar 10
1
Re: R-help Digest, Vol 13, Issue 9
C?dric Finet wrote:
>
> I thank you for your answer but I do not understand yet why the Fisher?s exact
> test does not work. And why is a "negative key".
>
> C?dric Finet
>
Running the original TOMS643 fortran code (R uses an f2c translation of
this) says:
FEXACT ERROR: 30
Stack length exceeded in f3xact. This problem should not occur.
The integer hash key is
2004 Nov 09
0
Re: QTL interval mapping in outbred populations
On Mon, 8 Nov 2004 ssim at lic.co.nz wrote:
>
> Is there an add-on package in R for QTL interval mapping for outbred
> population, eg. Haley-Knott regression method ?
>
> Stella
I believe Karl Broman's R/QTL (http://biosun01.biostat.jhsph.edu/~kbroman/software/)
and Brian Yandell's bim (http://www.stat.wisc.edu/~yandell/qtl/software/),
Richard Mott's happyR
2011 Jun 30
1
Looking for Filliben (correlation test)
Hi R-users,
I need to run Filliben correlation test, but I'm not able to find this
function in any package.
Perhaps it isn't in a package, but sometimes anybody have program this
function or create a routine to estimate it.
I found it in google and I hope what I found it's in this page,
http://genepi.qimr.edu.au/staff/davidD/R/filliben.R
but i can't open it.
Anybody can help me?
2005 Sep 27
1
Simulate phi-coefficient (correlation between dichotomous vars)
Newsgroup members,
I appreciate the help on this topic.
David Duffy provided a solution (below) that was quite helpful, and came
close to what I needed. It did a great job creating two vectors of
dichotomous variables with a known correlation (what I referred to as a
phi-coefficient).
My situation is a bit more complicated and I'm not sure it is easily
solved. The problem is that I must
2006 Mar 09
1
import from LISREL output of parameter estimates
I am using R and LISREL for simulation studies. R generates the data
that is analyzed with LISREL.
In LISREL I use "PV" in the LISREL output statement to request estimated
variances. LISREL writes these in a file that looks like this:
1 0 0
0.100331D+01 0.144845D+01 0.141009D+01 0.214423D+01 0.214129D+01
0.194464D+01
0.191531D+01 0.198328D+01 0.100683D+00-0.236392D-01
2013 Dec 17
1
Polychoric Principal Component Analysis (pPCA)
I have data set with binary responses. I would like to
conduct polychoric principal component analysis (pPCA). I know there are several packages used in PCA but I could not find one that directly estimate pPCA and graph the individuals and variables maps. I will appreciate any help that expand these reproducible scripts.
#How to conduct polychoric principal component analysis pPCA using
#either
2008 Sep 01
1
Polychoric and tetrachoric correlation
Hi there,
Am I correct to believe that tetrachoric correlation is a special case
of polychoric correlation when there are only two levels to the ordered
factor? Thus it should be okay to use hetcor from the polycor package
to build a matrix of correlations for binary variables?
If this is true, how can one estimate 95% confidence intervals for the
correlations? My guess would be
mat =
2013 Nov 28
0
polychoric correlation with multiple imputations, a Strate and a Weight
Hi there,
I'm generally more a Stata user than a R user, but I need to computed
something, and I am not able to do it with Stata 13. So, here I am!
I have a database that has multiple imputations (imputations are already
done) with a complex sample design (Strate and Weight).
Is it possible, in R, to run polychoric correlation with multiple
imputation, a Strate and a Weight ?
For the
2004 Nov 28
2
Tetrachoric and polychoric ceofficients (for sem) - any tips?
About two years ago there was a thread about this which suggested that at
that time nobody had these coefficients ready to go.
(a) has anyone in the meanwhile programmed them?
(b) I think I can see how to do the tetrachoric one with mvtnorm on similar
lines to an example on the help page so will try that if nobody else
already has
(c) looking at the polychoric one makes me realise yet again