similar to: pearson's correlation

Hello list, I'm trying to figure out how exactly the specification of nested random effects works in the lmer function of lme4. To give a concrete example, consider the rat-liver dataset from the R book (rats.txt from: http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/ ). Crawley suggests to analyze this data in the following way: library(lme4) attach(rats) Treatment <-

Hello, I asked this question a little while ago ( https://stat.ethz.ch/pipermail/r-help/2008-April/158761.html ) but got no response. Can anybody explain to me the difference between /, :, and %in% in the definition of random effects in lmer, such as: (1|A/B), (1|A:B), (1|B %in% A)? My understanding is that (1|A/B) is the same as (1|A) + (1|A:B), but I have not seen this stated explicitly

Colleagues I am working in a Linux OS with R 2.6.2. I need to execute a telnet session to another Linux machine from R, perform some operations, then return to the original computer. When I am in an R session, this is easy to accomplish by typing: system("telnet -l username machinename") I am then asked for my password; once that is entered, I am connected to the remote

Hi All, I am trying to lable the y-axis on my scatterplot with the following: "Soil moisture content (m3m-3)" I am using the following coding for plotting the graph: plot(soilmoisture~gradientlevel, xlab="Levels of droughting gradient", ylab="Soil moisture content (m3m-3)", bty="l", font.main="2", pch=16, las=1, cex.lab="1.13") I have

Hi, I have an unbalanced dataset on which I would like to perform a one-way anova test using R (aov). According to Wannacott and Wannacott (1990) p. 333, one-way anova with unbalanced data is possible with a few modifications in the anova-calculations. The modified anova calculations should take into account different sample sizes and a modified definition of the average. I was wondering if the

Hi, When I create an environment object with new.env() and populate it with values then how can I save it into an .RData file properly, so it can be loaded later on in a new session? Saving an environment object with save() or save.image() results in an error message when loading again: Error: protect(): protection stack overflow Regards, benjamin ======================================

hi, I'd like to do graphs like shown on the right panel on the following link: http://www.clinchem.org/cgi/content/full/53/8/1544/F1 (if link doesn't work see attached file) Any help and ideas highly appreciated, thanks in advance. Thomas von K?nel Human Genetics Uni Berne

Dear all I am new R user, and trying to learn more. I am doing linear regression analysis in R with my data. I am trying to find the way to calculate the slope value (coefficient of x) to degree of slope. Please give me idea on this. Thanking you in anticipation Warm regardMS _________________________________________________________________ [[alternative HTML version

Dear R-user Suppose I have the following data y=c(2,1,5,8,11,3,1,7,50,21,33,7,60) x=data.frame(y) for(i in 4:nrow(x)) x[i,] =var(x[i-3:i-1,]) I'm trying to get a new variable with the variance of the 3 previous values (just an example) and with NA in the three first positions. I know that my for() is wrong but I'm not able to find my error. Any idea? Thanks, Marlene.

Dear list, I encounter an error when I try to use ddply to generate means as follows: fun3<-structure(list(sector = structure(list(gics_sector_name = c("Financials", "Financials", "Materials", "Materials")), .Names = "gics_sector_name", row.names = structure(c("UBSN VX Equity", "LLOY LN Equity", "AI FP Equity",

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

Hello all; sorry for the (maybe) simple question... How can I join an Ubuntu (18.04) Work Station to an samba4 Active Directory (4.11.4) -- Rommel Rodriguez Toirac rommelrt at nauta.cu

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

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

How do I get non-squared correlation coefficient in some more sensible way than sqrt(summary(lm(y~x))$r.squared)? Thanks Matej -- Matej Cepl, matej at ceplovi.cz, Finger: 89EF 4BC6 288A BF43 1BAB 25C3 E09F EF25 D964 84AC 138 Highland Ave. #10, Somerville, Ma 02143, (617) 623-1488 The difference between death and taxes is death doesn't get worse every time Congress meets -- Will

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!

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

Rommel Rodriguez Toiracrommelrt at nauta.cu On 18/11/2020 17:34, Rommel Rodriguez Toirac via samba wrote:>?? > In my network I have a samba 4.11.4 as Active Directory Domain Controller installed in CentOS 7 (gtmad.gtm.onat.gob.cu - 192.168.41.17). I have recently installed samba 4.13.2 in CentOS 8 (gtmad1.gtm.onat.gob.cu - 192.168.41.18) and following the wiki.samba.org guide I have joined

Hi. Pre15 compiles out-of-the-box on AIX 4.3.2 ...almost. No patch included this time, but the following were the gotchas: - The __P() prototyping doesn't work (as discussed earlier) - bsd-daemon.o wasn't linked into libssh.a (though configure seemed to detect the need for it) DCE patch will follow shortly. Regards, Tor-?ke ______________________________________________________ Get

El 6 de abril de 2017 12:37:35 GMT-04:00, Marc Muehlfeld via samba <samba at lists.samba.org> escribió: >Hi Rommel, > >Am 06.04.2017 um 15:47 schrieb Rommel Rodriguez Toirac via samba: >> The problem is that I can not share directory using Windows > > or POSIX ACLs. Trying with Windows ACL: I use the Windows 7 > > RSAT. I use the Computer Management and the option