similar to: Throwing expressions around

Hi all, I have a problem to estimate Std. Error and t-value by ?polr? in library Mass. They result from the summary of a polr object. I can obtain them working in the R environment with the following statements: temp <- polr(formula = formula1, data = data1) coeff <- summary(temp), but when the above statements are enclosed in a function, summary reports the following error:

Hi all, I'm plotting to get the intersection value of three curves. Defining the x-axis as dsm, the following code works; dsm = c(800,600,NA,525,NA,450,400,NA,NA,NA,0) s3 = seq(0.05,1.05,0.1) plot(dsm,s3,col="blue",las=1,ylab="fraction",xlab="distance (km)") fc <- function(x,a,b){a*exp(-b*x)} fm <- nls(s3~fc(dsm,a,b),start=c(a=1,b=0)) co <- coef(fm)

Hi all, I'm trying to get the value of y when x=203 by using the intersect of three curves. The horizontal curve does not meet with the other two. How can I rectify the code below? Thanks Muhammad ts <- 1:10 dd <- 10:1 ts <- seq(200,209,1) dd <- c(NA,NA,NA,NA,1.87,1.83,1.86,NA,1.95,1.96) plot(ts,dd,ylim=c(1.5,2)) abline(lm(dd~ts),col="blue",lty=2)

Hi I want to draw a rectangle behind a barplot such that it highlights those particular bars from the rest of the plot. I have figured out how to draw a barplot(), and then how to draw a rectangle of the correct shape and size (using rect()), but when I use rect() it draws over the top of the bars, and then when I re-draw the bars, it draws with a white background, thus eliminating my rectangle!

Hello List I'm moving this over from the bioC list as, although the problem I'm working on is biological, the current bottle neck is my poor understanding of R. I wonder if someone would help me with the following function. cumulMetric <- function(deMirPresGenes, deMirs){ ??? #need to match position of each miR in deMirPresGenes with its FC to form a vector of FC in correct order ?

On Mon, Feb 29, 2016 at 6:17 PM, John Chambers <jmc at r-project.org> wrote: > The Wikipedia statement may be a bit misleading. > > S was never open source. Source versions would only have been available with a nondisclosure agreement, and relatively few copies would have been distributed in source. There was a small but valuable "beta test" network, mainly university

After upgrading my small test system from Debian Etch->Lenny via a complete reinstall, I find my g729 hostid has changed. Same machine, same CPU, same NIC! It doesn't seem reasonable that I have to burn my one "no-hassle" re-registration for a simple OS upgrade. The README only says that hostid is based on MAC addresses of all NICs, but that doesn't seem to be true. Does

Dear all, in the coming Winter Semester, I will be a teaching assistant for a course in Survival Analysis. My job will be to do the lab sessions. The software used for these lab sessions will be R. Most of the students have a background in social sciences and the only stats package they used so far is most likely SPSS. So I assume they might be quite surprised the first time they see R

Hi list, I have a little curve fitting problem. I would like to fit a sigmoid curve to my data using the following equation: f(x) = 1/(1 + exp(-(x-c)*b)) (or any other form for that matter) Where x is the distance/location within the dataframe, c is the shift of the curve across the dataframe and b is the steepness of the curve. I've been playing with glm() and glm.fit() but without

Dear All I am trying to convert from the type "expression" to the type "numeric". The following works: > x <- expression(6.2) > as.numeric(as.character(x)) [1] 6.2 However, the following does not work: > x <- expression(62/100) > as.numeric(as.character(x)) [1] NA Warning message: NAs introduced by coercion Any idea about how to deal with the second case?

A colleague is receiving some data from another person. That person reads the data in SAS and it takes 30s and uses 64k RAM. That person then tries to read the data in R and it takes 10 minutes and uses a gigabyte of RAM. Person then goes on to say: It's not that I think SAS is such great software, it's not. But I really hate badly designed software. R is designed by

Hi all, I am trying to run some plots on data, but when loading he CSV data file R is stopping and I am getting an out of memory error. Anyway to tweak this somehow to get it to run? Using WinXP with 4 GB RAM Tnx Bruce

I don't think anyone denies that you *could* make an EXE to do all that. The discussion is on *how easy* it should be to create a single file that contains an initial "main" function plus a set of bundled code (potentially as a package) and which when run will install its package code (which is contained in itself, its not in a repo), install dependencies, and run the main()

On Thu, Jul 26, 2018 at 12:22 AM, G?bor Cs?rdi <csardi.gabor at gmail.com> wrote: > I am trying to control a background R session, connected via a fifo / > named pipe. Is the fifo significant here? If I read the same R code from a file via `<` I get the input echoed (R 3.4.4, Ubuntu). Barry

I'm trying to plot an extrapolated logistic sigmoid curve using glm(..., family=binomial) as follows, but neither the fitted() points or the predict()ed curve are plotting correctly: > year <- c(2003+(6/12), 2004+(2/12), 2004+(10/12), 2005+(4/12)) > percent <- c(0.31, 0.43, 0.47, 0.50) > plot(year, percent, xlim=c(2003, 2007), ylim=c(0, 1)) > lm <- lm(percent ~ year)

If you index a vector with a vector that has NA in it, you get NA back: > x=101:107 > x[c(NA,4,NA)] [1] NA 104 NA > x[c(4,NA)] [1] 104 NA All well and good. ?"[" says, under NAs in indexing: When extracting, a numerical, logical or character ?NA? index picks an unknown element and so returns ?NA? in the corresponding element of a logical, integer,

HI there, I have been trying to use a code posted on R help to be able to calculate area under the curve for complicated data points and there seems to be an issue with the code: no "b" object found. I am not a good R user and can''t find were the problem is. Any help? Thanks!! This is the code ( as a test run I gave it this info because I know the answer: x<-seq(1:50)

A little while ago here we had a short discussion about Task Views - I think ignited by someone saying 'how many times do I have to say "have you read the Optimisation Task View?"?' and I poured some fuel on that fire by saying "Task Views" was a stupid name. Anyway, I did say that Task Views were rather brilliant, but were let down by their hidden position on the R

I can set up an entry widget (thanks to an old post by Barry Rowlingson) that gets a password and exits when the user clicks on the "OK" button. Anyone have any clever ideas for returning/ destroying the window when the user types a carriage return/ENTER in the text window? I've messed around a little with validate, validatecommand, but don't see any obvious way to do it ...

Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have done this]?? BUT, I would like a symbol like // (but rotated 90 degrees) to indicate that the y-axis has been truncated. I'd need to be able to