search for: rxx

Displaying 14 results from an estimated 14 matches for "rxx".

Did you mean: rx
2005 Jan 13
1
autocorrelation and levinson-durbin
...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 to some LD algo online. am puzzled with the reflection coefficient's computation. from the code, i understand it to be r = [- Rxx(i) - sum(lpc*Rxx)] / [error + 0.003* Rxx(0)] first question is why is negative of Rxx(i) added? from the resources i gather positive Rxx(i) is added. second is what is the motivation behind adding 0.003*Rxx(0) to error? 3) a more general question: i have not been able to find any reason why 10 lp...
2009 Jul 13
0
Partial Correlation
...data <- data.frame(x,y,z)* * * * * * if(corr_type=="pairwise")* * {* * * * print("inside pairwise")* * rxx.z <-P[as.numeric(element1_in_no),as.numeric(element2_in_no)]* * #print("rxx.z")* * #print(rxx.z)* * * * return(rxx.z)* * * * }* *...
2009 Jun 25
2
Error: system is computationally singular: reciprocal condition number
...hod) # is Szz positive definite? zz.ev <- eigen(Szz)$values if(min(zz.ev)[1]<0){ stop("\'Szz\' is not positive definite!\n") } # partial correlation Sxx.z <- Sxx - Sxz %*% solve(Szz) %*% t(Sxz) print(Sxx.z) # this gets printed rxx.z <- cov2cor(Sxx.z)[1,2] #some problem in this function function (V) { print("cov2cor") p <- (d <- dim(V))[1] if (!is.numeric(V) || length(d) != 2L || p != d[2L]) stop("'V' is not a square numeric matrix") Is <- s...
2009 Jun 28
1
ERROR: system is computationally singular: reciprocal condition number = 4.90109e-18
...ata,zdata,m=method) # is Szz positive definite? zz.ev <- eigen(Szz)$values if(min(zz.ev)[1]<0){ stop("\'Szz\' is not positive definite!\n") } # partial correlation Sxx.z <- Sxx - Sxz %*% solve(Szz) %*% t(Sxz) # this gets printed rxx.z <- cov2cor(Sxx.z)[1,2] # probably error in this function cov2cor return(rxx.z) } for some other datasets this pcor.mat function works fine. The solve.default function is this:--------------- solve.default<-function(a, b, tol = ifelse(LINPACK, 1e-07, .Machine$double.eps),LINPACK = F...
2004 May 13
2
please help with estimation of true correlations and reliabilities
...mpared to the true correlation matrix ct of t, however the diagonal entries remain without attenuation. Consequently the correlation matrix of observed variables has different things on- and off-diagonal. I would like to estimate 1) the true correlation matrix ct 2) the measurement reliabilities rxx, i.e. the correlation of a score with itself attenuated by its measurement error (as if we had two measurements of the same score). but I don't have predefined asumptions about structure in the variable set as I guess would be needed for SEM. Is R software available to do this? Best regards...
2015 Dec 21
2
Deutsche Telekom: calls dropped after 15 minutes
...22222 at default-0000003c;2", "2,Outgoing call for 0039015222222 using pbxluca") in new stack == Outgoing call for 0039015222222 using pbxluca -- Executing [0039015222222 at default:14] Dial("Local/0039015222222 at default-0000003c;2", "SIP/pbxluca/0039015222222,,RXx") in new stack == Using SIP RTP CoS mark 5 -- Called SIP/pbxluca/0039015222222 -- SIP/pbxluca-00000126 is ringing -- SIP/pbxluca-00000126 is making progress passing it to Local/0039015222222 at default-0000003c;2 -- Local/0039015222222 at default-0000003c;1 is ringing --...
2001 Feb 16
1
Sub_scribe and a question
...of a 220x220 Toeplitz matrix (loops)_______ (sec): ", = timing, "\n")) remove("b", "j", "k") flush.console() # (5) cumulate <- 0; p <- 0; vt <- 0; vr <- 0; vrt <- 0; rvt <- 0; RV <- 0; j = <- 0; k <- 0; x2 <- 0; R <- 0; Rxx <- 0; Ryy <- 0; Rxy <- 0; Ryx <- 0; Rvmax <- 0 # Calculate the trace of a matrix (sum of its diagonal elements) Trace <- function(y) {sum(c(y)[1 + 0:(min(dim(y)) - 1) * (dim(y)[1] + = 1)], na.rm=3DFALSE)} for (i in 1:runs) { x <- abs(rnorm(22*22)); dim(x) <- c(22, 22) ti...
2013 Jan 29
3
how to suppress the intercept in an lm()-like formula method?
...test version, 0.6-1 from R-Forge install.packages("candisc", repo="http://R-Forge.R-project.org") library(candisc) data(Rohwer) # this bombs: needs intercept removed cc <- cancor.formula(cbind(SAT, PPVT, Raven) ~ n + s + ns + na + ss, data=Rohwer) ## Error in chol.default(Rxx) : ## the leading minor of order 1 is not positive definite #this works as is cc <- cancor.formula(cbind(SAT, PPVT, Raven) ~ -1 + n + s + ns + na + ss, data=Rohwer) cc ## Canonical correlation analysis of: ## 5 X variables: n, s, ns, na, ss ## with 3 Y variables: SAT,...
2004 May 13
0
please help with estimation of true correlations and reli abilities
...mpared to the true correlation matrix ct of t, however the diagonal entries remain without attenuation. Consequently the correlation matrix of observed variables has different things on- and off-diagonal. I would like to estimate 1) the true correlation matrix ct 2) the measurement reliabilities rxx, i.e. the correlation of a score with itself attenuated by its measurement error (as if we had two measurements of the same score). but I don't have predefined asumptions about structure in the variable set as I guess would be needed for SEM. Is R software available to do this? Best regards...
2006 Aug 21
1
Retrieving p-values and z values from lmer output
I can't find a way to retrieve z values and p-values from the output from lmer in the lme4 package. How is this done? Rick B.
2006 Jul 04
1
lmer print outs without T
Hi, I have been having a tedious issue with lmer models with lots of factors and lots of levels. In order to get the basic information at the beginning of the print out I also have to generate these enormous tables as well. Is there a method command to leave off all of the effects and correlations? Or, do I have to go to string commands?
2019 Jul 02
3
How to restore the old network interface name?
Hallo, I need the device eth0 for one tool using centos 7.6. Using this tutorial https://www.certdepot.net/rhel7-restore-old-network-interface-name/ doesn t work. Thanks for a hint. Ralf
2015 Jul 10
2
[HCL] <Fideltronik INIGO> <Viper 1200> supported by <nutdrv_qx>
...u.pl/pub/nutdrv_qx_load.off.txt (4kB) load.on or shutdown.stop (C) - load are switched on after 10s, - stop olny "shutdown.return", - nutdrv_qx output: http://bpa.ksu.pl/pub/nutdrv_qx_load.on.txt (2kB) shutdown.return (S01R0001) - load is switched on immediately, the "Rxx" part does not matter, - "C" command can break the shutdown. - nutdrv_qx output: http://bpa.ksu.pl/pub/nutdrv_qx_shutdown.return.txt (7kB) shutdown.stayoff (S01) - "C" command will not break the shutdown. - nutdrv_qx output: http://bpa.ksu.pl/pub/nutdrv_qx_shu...
2006 Jul 15
3
names() function and lmer()
Hello All, I would like to retrieve some of the results from the lmer(...) function in library lme4. If I run a model, say fm.1 <- lmer(y ~ 1 + (1 | x), data = dog) and try names(fm.1), I get NULL. Is there anyway to retrieve the information? Thanks