Displaying 20 results from an estimated 300 matches similar to: "Partial Correlation"
2009 Jun 25
2
Error: system is computationally singular: reciprocal condition number
I get this error while computing partial correlation.
*Error in solve.default(Szz) :
system is computationally singular: reciprocal condition number =
4.90109e-18*
Why is it?Can anyone give me some idea ,how do i get rid it it?
This is the function i use for calculating partial correlation.
pcor.mat <- function(x,y,z,method="p",na.rm=T){
x <- c(x)
y <- c(y)
2009 Jun 28
1
ERROR: system is computationally singular: reciprocal condition number = 4.90109e-18
Hi All,
This is my R-version information:---
> version
_
platform i486-pc-linux-gnu
arch i486
os linux-gnu
system i486, linux-gnu
status
major 2
minor 7.1
year 2008
month 06
day 23
svn rev 45970
language R
version.string R version 2.7.1 (2008-06-23)
While calculating partial
2005 Jan 13
1
autocorrelation and levinson-durbin
hi,
am trying to understand speex's algo.
have a few questions.
1) autocorrelation:
in the function, _spx_autocorr (for floating point
version), there is a line
ac[0] += 10;
correct me if i am wrong, i suppose the addition of
10 is used to condition the autocorrelation matrix.
wonder how the value of 10 is arrived at?
2) levinson durbin (LD) algo
in the function _spx_lpc,
i referred
2011 Jul 20
0
The C function getQ0 returns a non-positive covariance matrix and causes errors in arima()
Hi,
the function makeARIMA(), designed to construct some state space
representation of an ARIMA model, uses a C function called getQ0,
which can be found at the end of arima.c in R source files (library
stats). getQ0 takes two arguments, phi and theta, and returns the
covariance matrix of the state prediction error at time zero. The
reference for getQ0 (cited by help(arima)) is:
2004 May 24
1
discriminant analysis
Hi,
I have done different discriminant function analysis of multivariat data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also need the predicted values to produce a plot of the analysis, as far as I know.
Here my code:
pcor.lda2<-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr,
2012 Jul 20
1
fitting Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD
Dear friends
i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD
method.
i found these formulas on
http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-model/
this is the mean-reverting process
http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt
and this is the script that i am using.......
ouFit.ML=function(spread) {
n=length(spread)
2009 May 11
0
Partial correlation function required
---------- Forwarded message ----------
From: <r-help-bounces@r-project.org>
Date: Mon, May 11, 2009 at 10:24 PM
Subject: The results of your email commands
To: das.moumita.online@gmail.com
The results of your email command are provided below. Attached is your
original message.
- Results:
Ignoring non-text/plain MIME parts
- Unprocessed:
What is the function for partial
2005 Dec 07
1
KMO sampling adequacy and SPSS -- partial solution
Dear colleagues,
I've been searching for information on the Kaiser-Meyer-Olkin (KMO)
Measure of Sampling Adequacy (MSA). This statistic is generated in
SPSS and is often used to determine if a dataset is "appropriate" for
factor analysis -- it's true utility seems quite low, but it seems to
come up in stats classes a lot. It did in mine, and a glance through
the R-help
2010 Mar 22
1
Dovecot 1.2.11 renames Maildir filenames on first read of foreign POPstore - how to avoid
Greetings,
I'm converting from a custom POPstore ( > 100,000 messages) which is
Maildir-like and uses the filename as the UIDL and doesn't store it
inside the message. I'm wondering whether there is hope for config-file
options to keep devecot from renaming files when it first sees them,
or whether I'm going to need to create a new backend for dovecot based
on Maildir.
2018 Sep 17
3
Cannot access HOME folder after upgrading to 4.8 from 4.6
Hello-
I upgraded Samba from 4.6 to 4.8 on a FreeBSD 11.2 server. After the upgrade, users cannot access the HOME folder share but they can access other shares just fine.
I am using the RID backend on this member server that connects to Windows-based domain controllers. I apologize for the lengthy smb4.conf but here it is:
#======================= Global Settings
2009 Sep 30
1
How to calculate KMO?
Hi All,
How do i calculate KMO for a dataset?
*Dataset:---------------------*
m1 m2 m3 m4 m5 m6 m7 m8
1 2 20 20 2 1 4 14 12
2 9 16 3 5 2 5 5 15
3 18 18 18 13 17 9 2 4
4 7 7 2 12 2 11 11 11
5 7 8 5 19 5 2 20 18
6 7 4 7 4 7 9 3 3
7 5 5 5 12 5 13 13 12
8 6 6 4 3 5 17 17 16
9 12 12 4 2 4 4 14 14
10 5 14
2008 Jun 24
2
persp plot
I have a set of data in the form
x1, y1, z1
x1, y2, z2
...
x1, yN, zN
x2, y1, z(N+1)
x2, y2, z(N+2)
...
x2, yN, z(2N)
...and so on...
xM, yN, val(M*N)
I have been trying to figure out how to get R to use this data in a
persp plot. So far the only thing that I can figure out to do is to
break the data file into three different files. The first file
contains the x-coordinate data:
x1
x2
2004 May 24
2
Manova and specifying the model
Hi,
I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics.
I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor?
If I type:
2008 Dec 08
1
partial correlation
Hej!
I have the following problem:
I would like to do partial correlations on non-parametric data. I checked
"pcor" (Computes the partial correlation between two variables given a set
of other variables) but I do not know how to change to a Spearman Rank
Correlation method [pcor(c("BCDNA","ImProd","A365"),var(PCor))]
Here''s a glimpse of
2002 May 11
2
Bug on Mac version of lm()?
Dear Mac users,
Hi, as you might have probably read the thread of
"[R] Rsquared in summary(lm)" on May 10, it seems that Mac version of
lm() seem to be working incorrectly.
I enclose the script to produce the result both for lm() and manual
calculation for a simple regression. Could you run the script and
report with the version of R, so I don't have to go through every
builds
2009 Nov 11
2
Partial correlations and p-values
I'm trying to write code to calculate partial correlations (along with
p-values). I'm new to R, and I don't know how to do this. I have searched
and come across different functions, but I haven't been able to get any of
them to work (for example, pcor and pcor.test from the ggm package).
In the following example, I am trying to compute the correlation between x
and y, while
2005 Nov 08
1
[PATCH] build warnings in mdf.c
Hi
I just upgraded to http://svn.xiph.org/trunk/speex r10357
and got this build warning:
alfredh@io:$ make -s mdf.o
libspeex/mdf.c: In function 'speex_echo_cancel':
libspeex/mdf.c:321: warning: statement with no effect
libspeex/mdf.c:317: warning: `adapt_rate' might be used uninitialized in this function
Is this intentional? In any case here is a simple fix:
Index:
2009 Jul 10
0
Error when running the examples in plotMultiTS / QRMlib
Dear List,
R ist somtimes a bit frustrating...
I ran the example from
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/QRMlib/man/plotMultiTS.html
and I get an error.
Has someone an explanation for this?
Thanks a lot,
Steve
CODE:
library(R.utils)
library(tseries)
library(wmtsa)
library(waveslim)
library(brainwaver)
library(timeSeries)
library(QRMlib)
> ?plotMultiTS
>
2008 Jul 01
1
[.data.frame speedup
Below is a version of [.data.frame that is faster
for subscripting rows of large data frames; it avoids calling
duplicated(rows)
if there is no need to check for duplicate row names, when:
i is logical
attr(x, "dup.row.names") is not NULL (S+ compatibility)
i is numeric and negative
i is strictly increasing
"[.data.frame" <-
function (x, i, j,
2013 Aug 26
2
Partial correlation test
Dear all,
I'm writing my manuscript to publish after analysis my final data with
ANOVA, ANCOVA, MANCOVA. In a section of my result, I did correlation of my
data (2 categirical factors with 2 levels: Quantity & Quality; 2 dependent
var: Irid.area & Casa.PC1, and 1 co-var: SL). But as some traits (here
Irid.area) are significantly influenced by the covariate (standard length,
SL), I