similar to: how to plot y-axis on the right of x-axis

Dear r-help, I am trying to plot several scatter plots with marginal histograms on one page. Ideally, a page is equally divided into 4 figure regions. Within each figure region, a scatter plot with marginal histograms will be plotted. I followed Dr. Paul Murrell's code released online to successfully plot the scatter plot with marginal histograms. The code applies "layout()" to

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

Hi all, I have a strange problem and rigth now I can't figure out a solution. Trying to calculate an ANOVA with one between subject factor (group) and one within (hemisphere). My dependent variable is source localization (data). My N = 25. My data.frame looks like this: > ML.dist.stack subj group hemisphere data 1 1 tin left 0.7460840 2 2 tin left

Why do I see this error? > library(stats) > require(stats) [1] TRUE > > tHSD <- TukeyHSD(aov) Error in TukeyHSD(aov) : no applicable method for "TukeyHSD" In case it helps: > aov Call: aov(formula = roi ~ (Cue * Hemisphere) + Error(Subject/(Cue * Hemisphere)), data = roiDataframe) Grand Mean: 8.195069 Stratum 1: Subject Terms: Residuals Sum

Dear Experts, I am new to R and unfortunately cannot start with a simply statistical analysis: I manually determined the volume of the right and left hippocampus in a group of meditators and in a group of controls. My data-sheet looks as follows: observation subject group age gender hemisphere volume 1 am04 m 25 f left 3.637 2 am04 m 25 f right 3.713 3 ao08 m 47 m left 3.715 4 ao08 m 47

I'm trying to read in a CSV, with lines looking like: HEADER, Latitude DecDeg, Latitude Hemisphere, Longitude DecDeg, Longitude Hemisphere, Speed knots, Bearing Degrees, fixQualityGga, noOfSatGga, altGga, heightGga, selectionGsa, fixGsa, pdopGsa, hdopGsa, vdopGsa, noOfSatGsv, Time, *checkSum $GPS, 3747.0224, N, 12223.4522, W, 0.36, 348.21, 1, 3, 4.01, 175.5, -25.2, A, 2, 4.14,

Hi, we need some help to define planned comparisons. I've based my understanding of the problem on reading Tabachnick and Fidell (2006), ie: http://www.ablongman.com/catalog/academic/product/0,1144,0205459382,00.html I don't understand how to specify planned comparisons in R. I've not found explanations for this in MASS or elsewhere. There is only discussion of the contrast

Hi I hope this is going to the right place. I am trying to write a program which uses KernSmooth library to estimate mutual information between two time series at various different lags. At the moment it’s producing negative values, which is supposed to be impossible (something is fishy). I am summing across one row of the matrix to get p(value is in bin x) and summing across the columns to get

Full_Name: Arni Magnusson Version: 1.6.2 OS: Windows XP Submission from: (NULL) (210.48.49.68) Hi there. I'm experiencing unexpected behaviour from as.POSIXct: > as.POSIXct("1969-12-24") [1] "1969-12-23 23:00:00 New Zealand Standard Time" > as.POSIXlt("1969-12-24") [1] "1969-12-24" > as.POSIXlt("1969-12-24")+1 [1] "1969-12-23

Hi, I'm new to R (and math ;) Would somebody please be so kind as to direct me in plotting a 3D sphere? I tried something in the lines of: #### y <- x <- seq(-pi, pi, length=pi*10) f <- function(x,y) { z <- sqrt(pi - x^2 - y^2) #z[is.na(z)] <- 0 z } z <- outer(x, y, f) persp(x, y, z, theta = 120, phi = 30) #### I've also tried: .... make.surface.grid(...) ..

First, a disclaimer. I am not affiliatied with the SAMM package. I am only a user of the package, but I have been contacted (off list) by people requesting information about SAMM and so I am posting this information here. SAMM is software for fitting mixed models. Versions are available for both S-Plus and R. More information and downloads of the software (and manual) are available here:

Dear Tom, I see that RPMS for R 3.4.0 are not successfully built on Fedora, nor on RHEL 7. https://koji.fedoraproject.org/koji/packageinfo?packageID=1230 When I build with mock on Fedora 25 I also get a build failure. R is not able to establish the time zone inside the chroot and this triggers a regression tests which is new in R 3.4.0. Specifically, in file tests/reg-tests-1d.R this test

With respect to calculating the epsilon index of sphericity for ANOVA, discussed on pp. 45-47 of: http://www.psych.upenn.edu/~baron/rpsych.pdf It notes that epsilon is not required for a repeated measures design with only k=2 levels, as the minimum value of epsilon (e) is given by: e = 1/(k-1) so for k=2, we have e = 1 (ie, no correction of the F test df; see p. 46). These notes apply to a

Hello, I am a french Ubuntu user since a few weeks. So, please excuse my English and my Technical knowledge. I've got a problem when I play games with wine : my computer crashes and reboots. I tried two games : wow and oblivion. Here is my config : ubuntu 9.04 ? Jaunty Jackalope ? 64 bits + Ext 4 Graphic card GeForce GTX 260 PLUS - 896 Mo DDR3 INTEL Core 2 Quad Q9400 - 2.66 GHz 3 x

Dear all, I have been trying to convert coordinates from longitude/latitude to utm but I got an error. As soon as the longitude coordinate is greater than 90, I get the folloowing error message: "error in pj_transform: latitude or longitude exceeded limits" Here is what I did: SP<-SpatialPoints(cbind(126.59,-14.30),proj4string=CRS("+proj=longlat"))

> On 18 May 2017, at 13:47 , Joris Meys <jorismeys at gmail.com> wrote: > > Correction: Also dlt uses the default timezone, but POSIXlt is not recalculated whereas POSIXct is. Reason for that is the different way values are stored (hours, minutes, seconds as opposed to minutes from origin, as explained in my previous mail) > I would suspect that there is something more subtle

Dear R users, I am currently working in subsetting a zooreg() object using either window or subset. I have a solution but it may be a bit cumbersome when I start working with actual data. Your inputs would be greatly appreciated. Example: I have a zooreg() object that starts in 1997 and ends in 2001. This object contains daily data for the 4 years

> On 18 May 2017, at 14:58 , Martyn Plummer <plummerM at iarc.fr> wrote: > > > >> On 18 May 2017, at 14:51, peter dalgaard <pdalgd at gmail.com> wrote: >> >> >>> On 18 May 2017, at 13:47 , Joris Meys <jorismeys at gmail.com> wrote: >>> >>> Correction: Also dlt uses the default timezone, but POSIXlt is not recalculated

In the autumn and winter it's considered bad form to set power state to D3cold, it might cause the device to freeze to death. This is also the case in the other seasons, so any device in the southern hemisphere is affected too. D3cold is not a valid state in a call to pci_set_power_state, only up to D3hot is allowed. Signed-off-by: Maarten Lankhorst <maarten.lankhorst at canonical.com>

I've posted this to SNOM, but was wondering wheter anyone here has issues with SNOM 190 phones not showing the correct DST adjusted time (using the latest firmware). -- ========================================== Rod Bacon Empowered Communications Ground Floor, 102 York St. South Melbourne Victoria, Australia. 3205 Phone: +613 99401600 ? ?Fax: +613 99401650 FWD: 512237 ? ? ? ? ? ? ? ? ?