similar to: Help boxplot to add mean, standard error and/or stadard deviation

I have the following problem. inf p= Int [dnorm(x-d) * pnorm(x) ^ (m-1)] dx -inf Given p and m, I want to find d. For example, p m d .9 2 1.81 Is there a way to solve this problem in R? Ideally I would have a function with arguments p and m, and output d. Thanks very much for any help. Bill -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Hi :) Was just wondering whether someone could help me with adjustments to trellis plots (parallel). I've got two way multivariate data. I want to make parallel plots for one of the factors, and want to color the lines according to the other factor. The first thing I manage, but with the other I'm lost :( Can only change the overall color. This is basically how far I get:

I would like to put confidence intervals on a mantel corellogram they are already calculated in the pmgram object but I am unsure how I get the x value in order to plot them? package(ecodist) X<-1:100 Y<-rnorm(1:100) Z<-rnorm(1:100) XY<-dist(data.frame(X,Y)) YX<-dist(data.frame(Y,X)) my.mgram<-mgram(XY,XZ) plot(my.mgram) print(my.mgram) > print(my.mgram) $mgram

Dear List, I am looking to plot error bars on a line using dispersion. I have values for the upper value and for the lower values, however i am unsure how to plot different values for the upper CI and the lower CI? I have been using dispersion(1:35,sim,simCI,col="red") Where there are 35 points, and simCI relates to a vector containing the lower confidence intervals, however i want

Hello all, I have an optimization routine that is giving me good results, but the results are not in the nice "model" format like "lm". How can I get optim results into a model so that I can use the clever 'fitted', 'residuals', and 'summary' functions? Using optim is the only way that I was able to make a model that 1) sums the betas to 1, 2)

Hello, I have an r function that creates the following dataframe tresults2. Notice that column 1 does not have a column heading. Tresults2: [,1] estparam 18.00000 nullval 20.00000 . . . ciWidth 2.04622 HalfInterval 1.02311 pertinent code: results<-cbind( estparam, nullval, t, pv_left, pv_right, pv_two_t, estse, df, cc, tbox, llim, ulim, ciWidth,

Dear Dr.Goslee and anyone may intrested in matrix manipulate, I am using your ecodist to do mantel and partial mantel test, I have locality data and shape variation data, and the two distance matrixs are given as belowings. When I run the analysis, it is always report that the matrix is not square, but I didn't know what's wrong with my data. Would you please help me on this. I am quite

I would like to perform partial mantel tests on only within group values, with "between group" values assigned to NA. This is possible in package ncf partial.mantel.test, however this sues a different permutation to that used in ecodist.ecodist will not accept data with NA values, returning a "matrix is not square error. is it possible to perform this test in ecodist? many thanks

Helo Greg - I got it with help of your code and I would like to add the value into the boxplot chart. How to do that. y_duration <- c (16.438, 15.321, 12.700, 12.397, 10.795, 9.928, 10.386) On May 09, 2018, at 10:06 AM, Jayaganesh Anbuganapathy <jai6561 at icloud.com> wrote: Actually I would like to get an output on the below snapshot. I have tried various method like points,

Got problems with plotCI (plotrix) I only want to plot the upper part of the error bar in my barplot I had the exact same commands working two months ago Now I wanted to change the legend and when I ran it again plotCI stopped working. To me it seems like there must some bug in R library(plotrix) library (graphics) x = matrix(c( 13.75516276, 3.944604404,

Dear R-users, I'd like to estimate standard errors (for lda) and make a boxplot with the "632plus" and "boot" error estimators included in package ipred (method: errorest). The "boot" estimator returns only a standard deviation but not the whole error data. Thank you in advance, regards, Antoine

What I want to do is create a horizontal bar chart matrix for a set of data that have a common set of variables (e.g., % of As, % of Bs, etc.) listed on the Y-axis and groups (e.g., Classes) on the X-axis. The key here is that the bars for each individual class plot are "centered" rather than left or right-justified. In archaeology plots similar to this are called battleship curves or

Hi, I have a data set like this: Mutant Rep Time OD 02H02 1 0 0.029 02H02 2 0 0.029 02H02 3 0 0.023 02H02 1 8 0.655 02H02 2 8 0.615 02H02 3 8 0.557 02H02 1 12 1.776 02H02 2 12 1.859 02H02 3 12 1.668 02H02 1 16 3.379 02H02 2 16 3.726 02H02 3 16 3.367 306 1 0 0.033 306 2

Hi all, I am have some difficulty with the legend function. I need to add a legend to describe the different line types in a plot. The legend box is small. It did not include sufficient length of each line type to help distinguish the differnt line types. Is there a way to fix this. Thank you Hannah [[alternative HTML version deleted]]

Dear HelpeRs, I have a script where I save several thousands of graphics. These are then used in Latex through Sweave. Unfortunately R crashes while making these plots and Windows pops up some message that I run low on virtual memory. I tried to save the plots using pdf(), windows() and postscript() and also tried to run it with R CMD BATCH myscript.R. But after a while R slows down and crashes

Have made a stacked area plot, but now want to manipulate the x axis: make in even increments (50) in order to suppress the tick marks (forming solid bar under plot). However, the plot functions do not seem to work, and I cannot find documentation for use of xat in stackpoly. This is a time series of data covering 579 years, and has been successfully converted from zoo to matrix. Data plotting

All, There are 3 bugfix updates for CentOS-4 (RC1). 1. mod_perl needed to be recompiled after the RH errata for perl was incorporated, but it was not flagged by the perl update. ( https://bugzilla.caosity.org/show_bug.cgi?id=803 ) (thanks Joshua Hirsh) 2. httpd identified itself as Apache/2.0.52 (Red Hat) Server instead of CentOS ( https://bugzilla.caosity.org/show_bug.cgi?id=806 ) (thanks

Dear all, I am fairly new to package Plotrix and I would like to ask you, if any of you could help me with following. a) I don't know how to set up border line width in the Stackpoly function (seems that lwd from par doesn't work, or at least not in the way I have written) b) I don't know, how to rotate the x axis labels (I used srt=45, but same problem as before). c) Could anyone

I have an issue using mle in versions of 32 bits. I am writing a package which I want to submit to the CRAN. When doing the check, there is an example that has an error running in the 32 bits version. The problem comes from the mle function, using it with a lower constrain. In 64 bits version it works fine but when I put it in the R 32 bits it fails. (same numbers, all equal!) The call is:

I am fittting a linner regression with constrained parameters, saying, all parameters are non-negative and sum up to 1. I have searched historical R-help and found that this can be done by solve.QP from the quadprog package. I need to assess the significance of the coefficient estimates, but there is no standard error of the coefficient estimates in the output. So I can not compute the p-value.