similar to: Re: normal distribution in samples of soil organisms.

Hi Sirs and Madams. My question is more statistical than related with the use of R software and I hope it will not seems so silly and elemental. I'm analyzing a set of data of some soil organism collected in diferent landscapes, soils taxa, and depths. The sample was performed thinking in a factorial structure with four factors: Specie, Landscape, Soil and Depth. Because not all the

Dear all, I have several sets of x-y-z points and I need to estimate the volume that encompass all my points. Recently I got some adivice to show the "convex hull" of my points using geometry package (see code below). But now I need to calculate the volume of my set of points. Any advice are wellcome. Thanks in advance Miltinho Astronauta Brazil ----- require(geometry) library(rgl)

Hi All, I would like to use Kennard-Stone algorithm for splitting a dataset. > mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/binary.csv ")) > library("soil.spec") > ken.sto(mydata,per.n=0.3) Error in ken.sto(mydata, per.n = 0.3) : subscript out of bounds I found that other people run into this problem as well:

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

http://r.789695.n4.nabble.com/file/n3031344/RSV.Rdata RSV.Rdata I want to split my dataset to training set and test set using kennard-stone(KS) algorithm, it is lucky there is R packages soil.spec to implement it. but when I used it to my dataset, it does not work, who can help me, how reasons is it, below, it is my code, and my data in the attachment.

Dear R users, has anybody made an attempt to create the soil texture triangle graph in R? For an example see here: http://www.teachingkate.org/images/soiltria.gif I would like to get the lines in black and texture labels in gray to allow for plotting my texture results on top. Any examples or suggestions are very welcome! Thanks in advance, Sander. --

Hello, Using the dataset below I produce a soil.texture plot (R code for this included at the end). One will notice I have plotted the points based on whether or not they are in a region called 'upstream' or 'downstream'. I'm curious if there is a way to somehow extract counts of the number of points (red points and blue points) falling within each classification (e.g. silty

Hi, I have a dataframe wf existing of a header with different labels and beneath the values of those labels : wf: label1 label2 ... 0,45 0,21 0,10 0,45 .... .... I have a list fl <- c("label2","label3",..) Isn't possible to use the list elements in the list in order to subtract values from the dataframe? like : wf$fl[[1]] When I do in R I get :NULL

Hi people, I have made some progress trying to work out how to solve this problem but I have got a bit stuck - sorry if this turns out to be a simple exercise . . Allelic Differentiation (AD) in genetics measures the number of different alleles between (say) two populations eg: Organisms in Pop 1 have alleles: a, b, c, d, e Organisms in Pop 2 have alleles: b, b, c, d, e Different

Hi, I wonder if you could help me with my barcharts, I am using the barchart function from the lattice library. I have data like this: Complex,Organisms,Percentage 130,0,50 130,1,10 130,2,20 130,3,15 130,4,5 133,0,10 133,1,15 133,2,20 133,3,50 133,4,5 I draw barcharts using this command...

Hello, i'm trying to convert my physical web servers to a virtual guest. What i'm experiencing is a poor disk i/o, compared to the physical counterpart (having strace telling me that each write takes approximately 100 times the time needed on physical). Tested hardware is pretty good (HP Proliant 360p Gen8 with 2xSAS 15k rpm 48 Gb Ram). The hypervisor part is a minimal Centos 6.5 with

Hi, I have a problem about graphics. I would like to plot two graphs: a barplot and curve. Here is the code: > barplot(dpois(0:45,20),xlim=c(0,45),names=0:45) > curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T) Both graphs are drawn in the same figure, however the scale in both graphs dooes not match. For some reason the second plot is shifted to left. I think there is a problem

Dear All, I am newbie to R, and I wanted to plot a barplots with R and in such a way that It will also show me position which I can plot on the bar line. Here is my code that I am using to plot, > chromosome <- c(40.2, 35.6, 36.1, 29.6, 31, 29.6, 31, 29.4, 28.2, 23, 23, 28.2) >barplot (chromosome, col="purple", xlab="Oryza sativa Chromosomes", border = NA, space =

Hi everyone, I want to plot a 3D interpolation of the concentration of aquatic organisms. My goal would be to have the result represented as clouds with a density proportional to the abundance of organisms, so that I could fly (well, swim actually ;) ) through the scene and see the patches here and there. Basically, I want to do something like this: http://www.youtube.com/watch?v=27mo_Y-aU-c but

Would someone be able to help with this question? I'm using the Gcmrec, Survrec, and Design packages to do a power analysis on simulated data. I'm receiving an error after using the Survr function that all data must have a censoring time even after using the gcmrec function: newdata<-addCenTime(olddata). My program is below. I'd greatly appreciate any help!

Hallo, I would like to organise at my pleasure the layout of a trellis plot. Currently I have a 3x3 matrix display and 7 plots. Is it possible to choose which specific panels will stay empty? I tried index.perm to arrange the order. Then there is perm.cond which I could not understand if it can serve for my purpose since I cant really find which kind of data it accepts. Thank you for the help.

Hi Guys, I have this part of a program: library(survival) Gastric <- cbind.data.frame(Gp=c(rep(1,45),rep(0,45)), ### 2nd gp 0 time=c(1,63,105,129,182,216,250,262,301,301,342,354,356,358, 380,383, 383,388,394,408,460,489,499,523,524,535,562,569,675,676, 748,778,786,797,955,968,1000,1245,1271,1420,1551,1694,2363,2754,2950,

Dear list, I am looking for a possibility to present results in a more graphical way by adding an axis. But I have trouble relating my data to the added axis. Imagine the following example: a <- c(10, 20, 30, 40) b <- c(50, 250, 500, 600) ba <- b/a par(las=1, mar=c(5,5,.5,5)) plot(a,b, type="b", pch=22, cex=2, col=4, lwd=2, ylim=c(0,650), xlim=c(0,45)) axis(4,

Greetings I need to merge two data sets. The first data set contains information on individual organisms roosting in clusters,e.g., cluster one has 5 individuals: 3 females, two males, and one juvenile and so on for hundreds of clusters. The second data set contains temperature data on each of the plots from the first data set. The problem: the second data set only has one data point per

Hello, I'm struggling with understanding the output on the ancestral.pars() command from the phangorn package, I'm new to doing phylogenetic analyses using R. I used it on nucleotide data, and it works fine, I'm just not sure how to read the output. The output is phyDat class, and outputs a matrix for each node/leaf in the tree. I figured out that the matrix columns represent the four