similar to: Obtaining x-values from ECDF

Dear R users, Quick question on the use of ecdf for producing a cumulative distribution function. After having produced a cumulative distribution using ecdf, I would like to extract the exact x-value for a given percentile. For example I can ascertain the 1st and 3rd quartile using the summary function: cawa.cdp <- ecdf(cawaocc$LEFF80) plot(cawa.cdp, ylab="Fn(x)")

Hello everyone, Quick question...is there a way of specifying a y-intercept value within a lm statement. For example, if I wanted to specify the regression to pass through the origin I would enter lm(y~0+x). But can I specify an actual term such as 1,2,3,4, etc. as an intercept value? Thank you! -- Benjamin Zuckerberg Doctoral Candidate State University of New York College of

Hello all! I have a problem to extract values greater that for example 1820. I try this code: x[x[,1]>1820,]->x1 Please help me! Thank you! The data structure is: structure(c(2.576, 1.728, 3.434, 2.187, 1.928, 1.886, 1.2425, 1.23, 1.075, 1.1785, 1.186, 1.165, 1.732, 1.517, 1.4095, 1.074, 1.618, 1.677, 1.845, 1.594, 1.6655, 1.1605, 1.425, 1.099, 1.007, 1.1795, 1.3855, 1.4065, 1.138, 1.514,

I am working with a RCB 2x2x3 ANCOVA, and I have noticed a difference in the calculation of sum of squares in a Type III calculation. Anova output is a follows: > Anova(aov(MSOIL~Forest+Burn*Thin*Moisture+ROCK,data=env3l),type=3) Anova Table (Type III tests) Response: MSOIL Sum Sq Df F value Pr(>F) (Intercept) 22.3682 1 53.2141 3.499e-07 *** Forest

Deare R users I am new to R. What I want to do is explained below;- I have table called States.Prob which is given below Prob of States Changes State1 State2 State3 State4 A Pa1 Pa2 Pa3 Pa4 B Pb1 Pb2 Pb3 Pb4 C Pc1 Pc2 Pc3 Pc4 D Pd1 Pd2 Pd3 Pd4 and I have a dataframe called

Hello, This is an urgent request. I want to convert a list of 3 elements into a matrix and I am not sure how to do it. The list looks like this: List of 3 $ : num [1:15364, 1] 0.133 0.622 0.588 1.024 0.583 ... ..- attr(*, "dimnames")=List of 2 .. ..$ : chr [1:15364] "6420681" "3610072" "2260458" "60689" ... .. ..$ : NULL $ :

Dear Users This is my dataset called mydata4. I want to sort the dataframe on the first column PxMid which is basically a column with dates. I've tried mydata4<-mydata4[order(mydata4$PxMid),] but it doesnt work. Could it be because these are dates? Please help I'm really stuck !! Thank you for your time. Regards Meenu PxMid EU0006MIndex.x DMSW1Curncy.x DMSW2Curncy.x DMSW3Curncy.x 1

Hi I have a question regarding an output of a binomial lmer-model. The model is as follows: lmer(y~diet * day * female + (day|female),family=binomial) The corresponding output is: Generalized linear mixed model fit by the Laplace approximation Formula: y ~ diet * day * female + (day | female) AIC BIC logLik deviance 1084 1136 -531.1 1062 Random effects: Groups Name Variance

Hi - I'm just starting out with R and have come up against a problem with what should be a simple operation. I'm plotting a range of safety standards for 7 different chemicals. I don't want the standard box and whisker plot - I just a dot for each value. Whilst using "plot.default" I get the error "invalid xlim value" - which is confusing since using the same script

Having spent the last few weeks trying to decipher R, I feel I may finally be getting somewhere, but i'M still in need of some advice and all my tutors seem to be on holiday! Basically a bit of background, I have data collected on a population of Lizards which includes age,sex, and body condition. I collected data myself this year and I have data previously collected from 1999, 2002 and

Hi there, I've been using the SMATR package to do standardized major axis (SMA) regression on allometric data, and I've been able to to pairwise comparisons of slope & elevation between multiple study sites. Now I'm trying to do the same thing using major axis (MA) regression, but I'm getting errors (SMATR should be able to do MA as well as SMA regression for everything,

Hi, I migrate Xen4.0 to XEN4.1, and found XEN4.1 performance is not as good as 4.0.2. xm shutdown takes longer time, control thru libvirt take longer time, and since I''m testing PCI/VGA passthru, the "createdevice: pci" also take around 5~8 seconds comparing to 1 second in XEN4.0.2. below is my configuration, do you have any suggestions to make the performance better?

I am trying to use R lm() with quantitative and qualitative predictors, but am getting different results than those that I get in SAS. In the R ANOVA table documentation I see that "Type-II tests corresponds to the tests produced by SAS for analysis-of-variance models, where all of the predictors are factors, but not more generally (i.e., when there are quantitative predictors)." Is

Hi: I am using a subset of the below dataset to predict PRED_SUIT for the whole dataset but I am having trouble with 'newdata'. The model was created with 153 records and want to predict for 208 records. wolf2 <- structure(list(gridcell = c(367L, 444L, 533L, 587L, 598L, 609L, 620L, 629L, 641L, 651L, 662L, 674L, 684L, 695L, 738L, 748L, 804L, 805L, 872L, 919L, 929L, 938L, 950L, 958L,

I am encountering an error when I attempt to reference a glm model within a function. The function uses the segmented.glm command (package = segmented). Within the segmented.glm command one specifies an object, in this case a logistic regression model, and specifies a starting threshold term (psi). I believe this is an environment problem, but I do not have a solution. Any assistance