similar to: Help in Customising NLS function to spit out Mean and SD of the new fit!!!

i'm coding in R(first time) for a paper my colleague is publishing. i plotted a histogram for 6000 values. Taking the mean and sd of the midpoints i did a dnorm and got the densities. pl<-dnorm(trimmedvals,mean=midsmean,sd=midsSD) (now in a loop of 5 times) i plotted these experimental values against theoretical values using the nls function the following way.

Hi all, My data genes [,1] [,2] [,3] [,4] [1,] 25 72 23 55 [2,] 34 53 41 33 [3,] 26 43 26 44 [4,] 36 64 64 22 [5,] 47 72 67 34 stu<-t.test(genes[,1:2],genes[,3:4]) > stu$p.value [1] 0.4198002 i get 1 pvalue for the entire col1:col2 Vs col3:col4. I am trying to get 5 p values for the 5 rows i have. I am trying to avoid a for loop coz my

Hi, my vector V<- c(1,0.5,0.06,0.056,0.01,0.04,0.4,0.9,0.82,0.1) if( V > 0.5) { V <- 1 - V } I get a warning saying only the first element will be used in comparing (if V > 0.5). However, my results tell me vis-versa ,it actually compares every element of the vector V with 0.5 and that is waht i want it to do. Using a for loop is expensive and time consuming since my actual

Hi all, I am trying to find sd of my rows in a matrix and i get column sd inspite of extracting rows. I tried to do the sqrt(var(x)) but that did'nt work as well, Here is my data genes 15 24 63 40 25 42 46 35 23 53 37 45 30 37 50 55 40 51 30 48 x<-sd(genes[1:5,]) y<-sqrt(var(genes[1:5,])) I get 4 sds for the 4 columns instead of 5 sds for my 5 rows. Thanks you in advance.

I'm trying to overlay a density plot on a previously plotted histogram. However, i need to use the same axis as of the 1st(histogram) plot to plot the second. My second plot is creating its own axis and causing my plot to extend the entire histogram instead of getting a subplot on a portion of the histogram. I tried 'fig' and 'new' parameters with no luck. Thanks in advance.

Hi all, I think the last time i posted this topic i started on the wrong foot. Thnaks alot to everyone who responded. i'm coding in R(first time) for a paper my colleague is publishing. i plotted a histogram for 6000 values. I am told to plot experimental vs theoretical vlaues from the histogram and do a non linear least square curve minimization and compute the mean and sd of the new x

Hi all, I'm trying to use the nls function to do a nls curve fit. I have a histogram for my data. I need to get athe best curve fit for this. i am plotting experimental vs the theoretical values. my experimental values are got from my dnorm function, mean=mean of my mid values(hist) and sd=sd of my mid values. For the theoretical i'm using mean and sd of all 20000 raw values as the start

How can I receive parameter estimates for a given curve that has been fit with the NLS function. The problem is that I have a "n parameter" curve and I want the optimal fit. The NLS procedure gives me final function values and not the individual parameter estimates that were used to define this "best" fit. What function can I use to get these parameters? Thanks

Good Day All ! I have a 2-D vector of mode numeric and a parallel 1-D vector of mode numeric. Here are my values. 2-D vector 1-D vector 80 75 85 80 80 70 80 90 80 80 60 70 80 70 70 85 75 95 85 85 70 60 90 60 70 My 1-D vector is the average across the rows of my 2-D vector. I process each column of the 2-D vector in a loop. When I trim my 2-D vector column

Hi all, I downloaded all the required r-tools from Dr.Brian Ripley's website. RCMD build dAnal It built the tar file for me. After which i tried RCMD check dAnal and i get the following error. Making package dAnal adding build stamp to DESCRIPTION file installing R files installing data files installing man source files installing indices Error in

Doesn't look like nonparametric fit to me, since nls() is used to fit to a gaussian density, so the result is a gaussian density (with estimated parameter). What I do not understand is why people would do this. This is not the first time I've seen people doing this, on both R-help and S-news (if my memory is still any good). If the objective is to fit a Gaussian distribution to the

Hi All, I get the following error when i perform a t.test. studentt<-apply(tlr, 1, function(x) t.test(x[2:6], x[7:11],var.equal=TRUE)$p.value) # tlr is a table of 11 columns and 22500 rows. I am not able to post tlr due to its size. Error in if (stderr < 10 * .Machine$double.eps * max(abs(mx), abs(my))) stop("data are essentially constant") : missing value where

Hi Group, I performed the following commands to build my package in R 2.0 under Windows XP I got all my tools from Dr. Duncan Mudroch's website. I did a RCMD build dnal and it built a tar file for me. I did a RCMD INSTALL dnal and it installed well. When i do RCMD check dnal i get the following 2 WARNINGS with no Errors. checking package directories..WARNING Subdirectory 'src'

https://bugzilla.mindrot.org/show_bug.cgi?id=3137 Bug ID: 3137 Summary: -f keeps stdin and stderr open Product: Portable OpenSSH Version: 8.2p1 Hardware: All OS: All Status: NEW Severity: normal Priority: P5 Component: ssh Assignee: unassigned-bugs at mindrot.org

I would like to customise my FCKeditor-instance so that when the "add image" button is pressed, it opens my own add_image popup. In this popup i would handle the image upload/image browser/etc. tasks with rails. How can I open the popup and pass the picture information from there back to FCKeditor''s textarea? -- Posted via http://www.ruby-forum.com/.

Hi, Does anyone have any advice on customising firmware of an SPA921 so that it can be locked to a sip provider and display logos on the config pages. Many thanks Dan Journo -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20091019/f6aa2510/attachment.htm

On 01/08/2015 11:18 AM, James Bishop wrote: > I apologise if this is not the appropriate list for the following > issue, but... > > Is it possible to upgrade the Linux kernel on a kickstart CD? > > I've changed the vmlinuz and initrd.img files in the isolinux > directory from the distro's ISO image; so kernel > 2.6.32-504.3.3.el6.i686 now boots from a DVD, and

I'm looking for an image for CentOS v7 to download.... Any suggestion on a suitable website? Sent from my iPhone > On Jan 8, 2015, at 1:30 PM, Greg Bailey <gbailey at lxpro.com> wrote: > >> On 01/08/2015 11:18 AM, James Bishop wrote: >> I apologise if this is not the appropriate list for the following issue, but... >> >> Is it possible to upgrade the

Hi, I want to customise the segments on an xyplot. Below is a simple example of what I'm trying to do... #Example dataset x <- c(-0.25, 0.25, 0.8) y <- c(-0.5, 0, 0.75) gp <- c("A", "I", "C") my.data <- cbind.data.frame(x,y,gp) #setting up the parameters to customise the lines with. ltype <- c(1,2,3) env.col <-

Newb here; I feel like im cheating when i actually write to the list, but i''ve been stuck on this for a while now and cant seem to find the answer in tutorials... heres the simplified version of the problem controller says: def show @exhibition = Exhibition.find(params[:id]) end content says: <% for exhibition in @exhibitions %> <%= link_to