search for: isotop

Happy Friday fellow R users. I need some help - I am trying to make a graph using ggplot 2 of some lead isotope ratios. Normally, the isotope mass number appears as a superscript before the chemical symbol. However, I cannot figure out how to do this ggplot2's axis labels. The closest I have come is using "." In front of the numbers I am superscripting e.g: ggplot(.... +xlab(expression(.^20...

I really like to use R for all my graphs, and as I work with stable isotopes I want to have a proper chemical notation in my plots I have looked at ?plotmath, but didn't find the answer and also searched the R website. ------------------------------------------------------------------------ -- plot(1:10,xlab=expression(^{14}*C)) # I want to have a superscript with...

I collected eggs laid by Springtails everyday over 28 days after swich to isotopically enriched diet. The eggs were pooled at day 7, 14, and 28 (+ day 0 = initial value) and analyzed for isotopes. After the diet switch the isotopic values of the adults and eggs change towards those of the new diet. Here are the d13C values (y) of the eggs: x y 0 -22.2 0 -22.2 0 -22...

Dear R-List, I would like to write nicely the names of some isotopes on a plot. The code bellow works fine. plot(1:10,1:10) text(c(2,4,8),c(2,4,8),labels=c(expression(italic(phantom(0)^{78}*Ge)), expression(italic(phantom(0)^{137}*Cs)), expression(italic(phantom(0)^{129*m}*Te))), cex=3 ) B...

Dear R-List, I would like to write nicely the names of some isotopes on a plot. The code bellow works fine. plot(1:10,1:10) text(c(2,4,8),c(2,4,8),labels=c(expression(italic(phantom(0)^{78}*Ge)), expression(italic(phantom(0)^{137}*Cs)), expression(italic(phantom(0)^{129*m}*Te))), cex=3 ) B...

...al" y=f(x). I only have the b0,b1's and R-squareds. Should i weigh the separate equations by their R-squared? thanks Remko ^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~' Remko Duursma, Ph.D. student Forest Biometrics Lab / Idaho Stable Isotope Lab University of Idaho, Moscow, ID, U.S.A.

...erent scales. Is this possible? I have searched through all the documentation I could find on the R homepage without luck and I like R very much and would rather not resort to using something like SigmaPlot for this seemingly simple thing. E.g. plotting the concentration of Oxygen and it's isotopic composition in the same graph. Thanks, Henrik ---------------------------------------------------------------------------- -------- Henrik Andersson Netherlands Institute of Ecology - Centre for Estuarine and Marine Ecology P.O. Box 140 4400 AC Yerseke Phone: +31 113 577473 h.andersson@nioo.kna...

Hi there, I am trying to run an MCMC on stable isotope data from certain organisms to determine their dietary habits in the package SIAR. I have prepared my data according to Inger, R., Jackson, A., Parnell, A., Bearhop, S. : SIAR V4 (Stable Isotope Analysis in R) an Ecologist's Guide (also better known as 'SIAR for dummies'), which is as...

...e much better! Anyway, I'm trying to use R to run a gamma distribution model to estimate mean transit times of water moving through a hydrological catchment. My input are 3 .txt format files as follow: Precipitation_18o: 2 columns - column 1 is a date (Excel number format) and column 2 is an isotope ratio (i.e. -8.12) Runoff_18o: same as above Daily_Precip: 2 columns - column 1 is the same date format but column 2 is a weekly bulk precipitation value (i.e. 10mm) When running the script, I keep getting the following error message: Error in data.frame(cQ.ou[ind.mea], cQ[cal.cQ:nrow(cQ), 2]) :...

...)$x > eqscplot(aa.ld,type="n",xlab="LD1", ylab="LD2",las=1) > text(aa.ld,c(rep('f',13),rep('b',10),rep('p',10))) > aa.mean<-lda(aa.ld,AA$group)$means > points(aa.mean,pch=3) Best, Thomas Larsen Leibniz-Laboratory for Stable Isotope Research Max-Eyth-Str. 11-13, 24118 Kiel, Germany Work: +49-431-880-3896 -- View this message in context: http://www.nabble.com/Decision-boundaries-for-lda-function--tp24823213p24823213.html Sent from the R help mailing list archive at Nabble.com.

Dear R-List, I want to use the label function (from Hmisc library) to allow for the names of my isotopes. library(Hmisc) library(lattice) library(grid) num <- c("78","137","129m") nom <- c("Ge","Cs","Te") df <- data.frame(GE78=seq(nom),CS137=seq(nom),TE129m=seq(nom)) if I use this function to create the labels : lab <- functi...

Help! model1 <- siarmcmcdirichletv4(data, source, tef, concdep=0, 500000, 50000) Error in matrix(1, ncol = (numsources + numiso) * numgroups, nrow = (siardata$iterations - : invalid 'ncol' value (too large or NA) In addition: Warning message: In Ops.ordered((numsources + numiso), numgroups) : '*' is not meaningful for ordered factors Previous message, but more problems:

...atex" as well as "pdftex" both run in the (Win) command window. I use R1.81 but also tried in R1.71. Any ideas? Remko ^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~' Remko Duursma, Ph.D. student Forest Biometrics Lab / Idaho Stable Isotope Lab University of Idaho, Moscow, ID, U.S.A.

Hi all... I?m making a chart dealing with frequencies of isotopes of various elements. For instance, I'd like the following text to appear on a chart with the "35" and "37" as superscripts: Based upon: 35Cl: 75% 37Cl: 25% I am having problems properly parsing the superscript that preceeds the "Cl", since there is no characte...

...subset the x and y data appropriately. In the toy problem below I want to plot each inverted V (?model?) as LINES, with a single POINT (?data?) in the center of each inverted V, the same color as the inverted V. The code runs, but I can?t seem to mix lines and points. In my real problem (stable isotopes with ellipses superimposed on data) there will be additional panels but I am creating only one panel here, for simplicity. #generate test "model results" x<-1:9 y<-rep(c(1,3,1),3) model<-c(rep("a",3),rep("b",3),rep("c",3)) modelresults<-data.fr...

Hi, I have network data in the form of a couple of edgelists containing weights in the format "x,y,weight" whereby x represents row header and y represents column header. All edgelists are based on links among 634 nodes and I need to convert them into a 634*634 weighted matrix. I searched for online help using possible keywords I know, but I could not find a clue how to do this in

Hello, I'm using R code that includes a residual permutation that was written as a supplement to the paper: Turner et al. 2010. A general hypothesis-testing framework for stable isotopes ratios in ecological studies. Ecology 91:2227-2233. The supplemental code is available at: http://www.esapubs.org/archive/ecol/E091/157/suppl-1.htm When I execute the function, no warnings are given (seems fine). But the resulting output matrix is simply a matrix of repeating rows (i.e., it app...

Hello All... I gave a task to my students that involved using mclust to look for clusters in some bivariate data of isotopes vs various mining locations. They discovered something I didn?t expect; the data (called tur) is appended below. p <- qplot(x = dD, y = dCu65, data = tur, color = mine) print(p) # simple bivariate plot of the data; looks fine mod1 <- Mclust(tur[,2:3]) mod1$G mod2 <- Mclust(tur[,3:2])...

...ed to naAct in this way. Has anyone encountered this before? What to do? (i can, ofcourse, refit my model with excluding the NAs myself - but thats not so elegant!) Thanks, Remko Duursma ##################################### Remko A. Duursma, Ph.D. student Forest Biometrics Lab / Idaho Stable Isotopes Laboratory Department of Forest Resources University of Idaho - College of Natural Resources Email: remkoduursma at yahoo.com Tel: (208)-885-5165 _____________________________________________________________ Get 25MB, POP3, Spam Filtering with LYCOS MAIL PLUS for $19.95/year. http://login.mail....

...tobit model be modified to do this? Does anyone have experience with this type of dataset? Help is much appreciated, Remko Duursma ^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~'^'~,_,~' Remko Duursma, Ph.D. student Forest Biometrics Lab / Idaho Stable Isotope Lab University of Idaho, Moscow, ID, U.S.A.