search for: km2

Hello, This is Elaine. I am drawing a plot with x-axis label with "km square" as the unit. Now I want to print km square in the form of km2 and output "2" as the uppercase. Please kindly help suggest command to show the uppercase. Thank you. Elaine [[alternative HTML version deleted]]

...el script que he utilizado en el análisis. https://www.dropbox.com/s/d96itird8ms42yx/dataframe.Rdata?dl=0 <https://www.dropbox.com/s/d96itird8ms42yx/dataframe.Rdata?dl=0> sapply(levels(df0$bmi0),function (x){ #####SURVIVAL CURVE dfx=filter(df0,bmi0==x) surv2=Surv(dfx$bmi,dfx$MetS) km2=survfit(surv2~dfx$apoe4)##start.time=20,type='kaplan') plot(km2,lty=2:1,xlim=c(20,41),xlab='BMI at onset',main=x,mark.time = F) legend('bottomleft',c('E4','no-E4'),lty=2:1) cox=list(coxph(surv2~relevel(dfx$apoe4,ref='no-E4'))) }) sapply(levels(...

...of nine colors? 3. I have one group who's y range is much greater than all the others. I would like to split the plot somehow so that the bottom part shows ylim=c(0,200) and the top shows ylim=c(450,550). Is this possible? What I have so far is: library(lattice) xyplot(m.dp.area$Area.km2 ~ m.dp.area$DataPoint, m.dp.area, groups = m.dp.area$Manta, main = "Cummulative area of 100% MCP", xlab = "Data Point", ylab = "MCP Area", ylim = c(0,150), scales = list(tck = c(1, 0)), #Removes tics on top and r-axis pch...

...exy) <- c ("xCord", "yCord") ID <- mysample$ID datltraj <- as.ltraj(mysamplexy, mysampletime, id=ID) Ddat <- BRB.D(datltraj, Tmax=21600, Lmin=36) BRBdat <- BRB(datltraj, D= Ddat,type=c("UD"),Tmax=21600,Lmin=36, hmin=100) kernel.area(BRBdat, unout=c("km2")) # unfortunately my data are not a very good example. Sorry about that. But I guess you know what I mean. # Because I wish to run the functions three times I do want a data frame (or matrix) in the end, which includes the results (kernel areas): That means 12 columns and the three lines...

...uot; ) datltraj <- as.ltraj( mysamplexy, mysample$mysampletime, id=mysample$ID ) Ddat <- BRB.D( datltraj, Tmax=21600, Lmin=36 ) BRBdat <- BRB( datltraj, D = Ddat, type = "UD", Tmax = 21600, Lmin = 36, hmin = 100 ) result[ , i ] <- kernel.area( BRBdat, unout = "km2" ) } #> Warning in kernel.area(BRBdat, unout = "km2"): The grid is too small to allow the estimation of home-range #> for the following value of percent: 70,75,80,85,90,95. You should rerun kernelUD with a larger extent parameter #> Warning in kernel.area(BRBdat, unout = &qu...

Dear HelpeRs, I would like to do sg similar to: plot(c(0,100), c(0,100), xaxs="i", yaxs="i", type="n") txt <- paste("a =", a, "%\n b =", b, " km2", sep="") text(95, 95, txt, adj=c(1,1)) just with the "km2" formatted with 2 in a superscript. I thought txt <- substitute("a = "*a*" %\n b ="*b*" km"*phantom()^2, list(a=a,b=b)) would do it. However it does not work. Is it possible to fo...

Hello R list, I need to compute empirical variograms using data from a large geographic area (~10^6 km2). Although I could not find a specific reference, I assume that both geoR and gstat calculate distances among data points assuming points are on a flat surface (using the Pythagorean Theorem). Because the location of my data is large and located near the pole, assuming that latitude and longi...

...ccuracy measure. In the same way the 2 observations replicate1_subreplicate2 of model A and B (1 for A and 1 for B) are also linked. Is there a way to introduce such pairing in a mixed effect model? 2) several continuous covariates, attributes of species (number of points for modelization, size in km2 of the range) may influence the measures of accuracy and I may be interested in investigating those effects. How could I include a covariate in such model? How should I strucutre it given that there are species covariates (high order of nesting)? I hope that my questions are relevant! Thank you ve...

...", "yCord") > ID <- mysample$ID > datltraj <- as.ltraj(mysamplexy, mysampletime, id=ID) > Ddat <- BRB.D(datltraj, Tmax=21600, Lmin=36) > BRBdat <- BRB(datltraj, D= Ddat,type=c("UD"),Tmax=21600,Lmin=36, hmin=100) > kernel.area(BRBdat, unout=c("km2")) > > # unfortunately my data are not a very good example. Sorry about that. But > I guess you know what I mean. > > # Because I wish to run the functions three times I do want a data frame > (or matrix) in the end, which includes the results (kernel areas): That > means...

...1.000 0.8575587 1.000 0.9406087 0.778 0.5861066 0.600 0.6204616 1.000 0.8585725 0.000 0.2949169 0.100 0.1291246 0.444 0.7627612 observed = proportion of positive responses to crop damage questionnaires in a 25 km2 grid predicted = values produced by the final glm(binomial) model on the same dataset as used to develop the model Thanks a lot in advance for any suggestion Ahimsa Ahimsa Campos Arceiz The University Museum, The University of Tokyo Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033 phone +81-(0)3-5841-28...

...t;)] d_ltraj <- as.ltraj(dxy, d$myframestime, id = d$AnimalID) dLiker <- liker (d_ltraj, sig2 = 18, rangesig1= c(1, 100) ) dBB <- kernelbb(d_ltraj, sig1=dLiker$Ernie$sig1, sig2=18) data.frame(AnimalID = unique(d$AnimalID), day = unique(d$day), h = kernel.area(dBB, unout=c("km2"))) } ddply(myframe3, .(day), hrfun) # i.e. I calculated the size of the homerange for Animal "Ernie" at every single day. (The several numbers per day # represent the several Homerange levels). # What I want to do now is to calculate the homerange for Animal "Ernie" fo...

Dear exports,I have created a hypsometric curve (area-elevation curve) for my watershed by using simple command hypsometric(X,main="Hypsometric Curve", xlab="Relative Area above Elevation, (a/A)", ylab="Relative Elevation, (h/H)", col="blue")It plots the hypsometric curve in "RGraphics window", My question is how can I export values which is used

Hello, I tried to construct my very first loop today and completly failed :-( Maybe someone can help me? I have a dataframe somewhat like this one: myframe <- data.frame (Timestamp=c("24.09.2012 09:00", "24.09.2012 10:00", "24.09.2012 11:00", "25.09.2012 09:00", "25.09.2012 10:00", "25.09.2012