similar to: ERROR: system is computationally singular: reciprocal condition number = 4.90109e-18

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)

Why do we get Partial correlation values greater than 1? I have used the default function pcor.mat :-- I have manipulated the default pcor.mat function a bit so ignore tha variables corr_type,element1_in_no,element2_in_no,P.Please ignore the ?pairwise? section and have a look at athe ?listwise ? part i.e else part. *pcor.mat <-

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,

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:

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

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

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

Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados,

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

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

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

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

Hi All, Can someone tell me if solve function shown below for my version of R is proper or not? I am using R 2.7.2 .Wherever i have used this function ,i got results which were different from the expected results as computed using SPSS. Description of this function says:-- Solve a System of EquationsDescription This generic function solves the equation a %*% x = b for x, where b can be either

Hi, I'm having problems while working with large data sets with R 1.5.1 in windows 2000. Given a integer matrix size of 30 columns and 15000 rows my function should return a boolean matrix size of about 5000 rows and 15000 columns. First of all I tried to run this function on computer with 256 MB of RAM. I increased memory limit of R with memory.limit() up to 512 MB. I was inspecting

Hello everyone My name is Vangelis and I want to ask a question about partial correlation. I have used the command "pcor.shrink" to evaluate the partial correlations of a data.frame but the problem is that in the output results I cannot see whether these correlations are significant or not. Is there any command which can show me if these correlations are significant at 95% level or

I've noticed unstable behavior of nlminb on some Linux systems. The problem can be reproduced by compiling R-3.5.2 using gcc-8.2 and running the following snippet: f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) opt <- nlminb(rep(0, 10), f, lower=-1, upper=3) xhat <- rep(1, 10) abs( opt$objective - f(xhat) ) < 1e-4 ## Must be TRUE The example works perfectly when

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:

Hello, after thorough searching of the R help files as well as S+-help, I'm coming to the list: Is there a possibility to compute partial correlation coefficients between multiple variables (correlation between two paired samples with the "effects of all other variables partialled out")? All I seem to find are the standard Pearson correlation coefficients (with cor()) and no clue

r-help: I have an R package suggestion. After spending several hours the other day installing about a dozen packages, I had an idea. In xemacs, there is a "sumo" package which allows me to install a large bundle of xemacs packages at one time (about a 120 modes including ESS). I think R should have a similar bundle. It would be so much easier than hunting/downloading/installing.

Hello everyone I tried to install the library GeneNT in order to use the command pcor.confint because I want to construct confidence intervals for partial correlations but among other demanding the specific library needs the library "Graph" which I don't have it and I cannot find it at this site. Is there any other site that I can download this library? Thanks Kind regards