similar to: Response surface plot

Colin, If your propose is to create a ternary plot with points and vectors, I think easier do this with graphics based plots instead of trellis based plots. Although, with little work you can do with trellis too. I gave you a reproducible code to put an arrow in a ternary plot. I use the function locator() to extract coordinates. The code is the following

I'm in e-mail contact with a geochemist who maintains a well regarded geochemistry web site. He's drawing diagrams either with a Turbo Pascal program running in a DOS window or with Excel. I'm trying to persuade him that R would be a better choice. Something he's particularly keen on is ternary diagrams. I think he is talking about ternary phase diagrams, and if so, it looks as

I'm having trouble using the read.mtp function in the foreign package to import datasets from minitab (.mtw) format. Specifically, each file I try to import fails to load any data beyond a header row stating the version of Minitab that saved the dataset. I get this error: incomplete final line found by readtableHeader on 'income.mtw' The dataset appears to be complete (ie, using

Dear list: I encountered some problems using the function point.in.polygon() of the sp package, when trying to determine whether some points lye inside, outside, on the border or on a vertice of a polygon. I have a list of point I know should lye right on the border of a polygon, but some of them are not classified as such by point.in.polygon() (see the example code below). To make a long story

Hi I am performing an analysis of variance with two factors, each with two levels. I have differing numbers of observations in each of the four combinations, but all four combinations *are* present (2 of the factor combinations have 3 observations, 1 has 4 and 1 has 5) I have used both anova(aov(...)) and anova(lm(...)) in R and it gave the same result - as expected. I then plugged this into

Hello I am trying to add a ternary plot as a corner inset graph to a larger main ternary plot. I have successfully used add.scatter in the past for different kinds of plots but It doesn't seem to work for this particular function. It overlays the old plot rather than plotting as an inset. Here is a simple version of what I'm trying. Note that if I change the inset plot to be an ordinary

Dear All, Would like to ask the inconsistency in the autocorrelation from R with SPSS/Minitab. I have tried a dataset x with 20 data (1-20) and ask R to give the autocorrelation of different lags using the command < acf(x, lag.max=100, type = "correlation"), However while SPSS and Minitab give the same answers (0.85 for lag1), R gives 0.3688 which is much smaller. Obviously, the

Dear R list, I have been trying to do a linear model, extracting the effect of a covariate.... and the results do not match, when I do it with other programs (e.g. minitab).... so it is obvious that I was doing something wrong. Whan I do it with minitab, I have this results: (sector is a factor and depth is the covariate): Source DF Seq SS Adj SS Adj MS F P

Hi, I ran a small data set from a factorial experiment through R, Minitab and JMP... the result from R is significantly different from what Minitab or JMP give... The data set is at the following link: http://www.personal.psu.edu/nug107/Uploads/2x3_16repsANOVA.txt The first 5 columns are the factors and the next three are responses. In particular, for the response beta11MSE, two of the

Hi R-help. I am trying to run a linear mixed model with nested factors with either lme or lmer and I am having no luck obtaining the same results as Minitab. Here is Minitab's code: MTB > GLM 'count' = site year replicate(site year) site*year; SUBC> Random 'year' 'replicate'; Can you tell me how to code this in R? The settings are typeII, Tukey,

Hello, Given a probability function: p(x) = 12 / (25(x+1)) , x=0, 1, 2, 3 we generate the following values: C1 C2 0 0.48 1 0.24 2 0.16 3 0.12 Now, I'm supposed to create 50 random values using this table. In MiniTab, I simply selected Calc -> Random Data -> Discrete, and selected the columns, and it created 50 random values in a new column.[1] How do I do the

A strange thing happened. After installing the most recent RC version of 2.7 (2008-04-20 r45403) I compiled it as usual with no errors. After installing it I ran "make fullcheck". I got as far as testing recommended packages and I got the following: ... -------- Testing package foreign -------- Running examples in 'foreign-Ex.R' ... Comparing `foreign-Ex.Rout' to

I recently downloaded R for Windows, running on Win XP. I'm trying to create a perspective plot but not having any luck after reading the R manual and several examples found on the Internet. I have a 100 x 100 matrix of Z data as a tab-delimited text file exported from Minitab. I read this in to R using read.delim; this seemed to go ok. I created X and Y using seq() to get 100 divisions

Does anyone else miss email-friendly character graphics such as the following example, produced using Minitab? Histogram of C6 N = 478 N* = 21 Each * represents 2 observation(s) Midpoint Count -12 16 ******** -11 53 *************************** -10 63 ******************************** -9 83

Hi, I am trying to create a triangular plot to show the 'composition' of a set of items with three variables (historically the percent sand, silt and clay in soil). So far I have tried the 'soil texture triangle plot' in the package plotrix and the 'ternary or triangular plots' in the package cwhtool (cwhmisc). Both have strengths and weaknesses, but neither has a

I want to specify a two-factor model in lme, which should be easy? Here's what I have: factor 1 - treatment FIXED (two levels) factor 2 - genotype RANDOM (160 genotypes in total) I need a model that tells me whether the treatment, genotype and interaction terms are significant. I have been reading 'Mixed effects models in S' but in all examples the random factor is not in the main

I need to plot multiple sets of data, on one Ternary plot (from the VCD package). Since they all need to be formatted differently, I can't just coerce them into one matrix, I need three separate calls. Is there a way to do this? Thanks -- View this message in context: http://www.nabble.com/Plot-multiple-datasets-on-a-VCD-ternary-graph-tp18412092p18412092.html Sent from the R help mailing

hi folks, Bit of a newbie, but I've spent a fair bit of time looking for an answer on this, with no joy. Can anyone help me? Dataset: A single column of values in a csv file (eg. 52, 53, 54, 85, etc) Goal: In Minitab, you have what they call a dot plot. It's a histogram, where a single dot represents a set of identical values (eg. 57, 57, 57 would be one dot). Multiple dots are

Dear R-Users, I'd like to have some tips for a ternaryplot ("vcd"). I have this dataframe: a<- c (0.1, 0.5, 0.5, 0.6, 0.2, 0, 0, 0.004166667, 0.45) b<- c (0.75,0.5,0,0.1,0.2,0.951612903,0.918103448,0.7875,0.45) c<- c (0.15,0,0.5,0.3,0.6,0.048387097,0.081896552,0.208333333,0.1) d<- c (500,2324.90,2551.44,1244.50, 551.22,-644.20,-377.17,-100, 2493.04)

Hi all, I am drawing a ternary graph. Everything is fine with both ternaryplot (package vcd) and triangle.plot (package ade4), but I want to present scales in neither percents nor from 0 to 1 (this is actually the only option I found in both functions). I want the scales to be in a natural scale (from 0 to k, k being the number of objects). Is it at all possible? (Descriptions of both functions