similar to: Least square minimization (non-linear)

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

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, 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

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, 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

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.

Hi all, I'm trying to customize the nls function code to spit out the sd and mean of the new fit. All it gives me are the fitted values. THis function seems to work fine for my project otherwise. All i need is the SD and MEAN. Can someone help me to customise. I frankly did not understand the nls code at all. I have 10 x and 10 y values from a histogram. I need to use nls to do the best fit

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

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'

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

Mar 29 10:57:02 k8ux dovecot: POP3(stf): mbox sync: Expunged message reappeared in mailbox /var/mail/stf (UID 123 < 60016, seq=2, idx_msgs=0) Mar 29 10:57:02 k8ux dovecot: pop3-login: Login: user=<stf>, method=... Mar 29 10:57:03 k8ux dovecot: POP3(stf): mbox sync: UID inserted in the middle of mailbox /var/mail/stf (60016 > 123, seq=2, idx_msgs=1) Mar 29 10:57:03 k8ux dovecot:

Hi Do you have an idea why am I getting random reboots and (sometimes) this stacktrace? : http://people.eisenbits.com/~stf/load/20151222-kernel-stack-trace.png This should probably go as a bug report, but I have already uninstalled the driver --- sorry. This is a fresh install of Debian 8.3 (latest stable) (Jessie). After switching to VESA everything is stable --- the only problem is that this

Hi, I need some help understanding the output from PAM. When I look at the output it doesn't list the cluster number by the median vlaues on each of the variables (like it does with k-means) Instead I have the following: So I know for instance cluster 1 has a mean for variable1 of 33.33, however when I run PAM i get: variable 1 variable2 29 32 12 97 12

Hi, The colClasses seem to be supressing 'NA' vlaues. How do I fix this? R script and first 5 lines of output is below. File "test2.dat" has blanks that are read as "NA" when I do not use 'colClasses', but as blanks when I use 'colClasses'. temp.df <- read.fwf("test2.dat", width=c(10,1,1,1,1,2,2,3,3,1),

Was having problems getting NT Workstation to see the samba server as the PDC and do domain logins (Win9x does just fine). Then found this comment in Sam's Teach Yourself Samba in 24 hours with regard to configuring a samba PDC: "You need to configure the server to be in user-level security and you must enable password encryption. There is no way around either of these settings."

Hello, I would like to use by to create a series of plots, but I'm not sure how to design the function. Here's what I've got: 1. A data frame of observations of measured value 2 measured vlaues and a pch code. V1 V2 pchCode .0045 123 1 .0034 145 2 .0045 123 1 .0046 167 3 ... So what I want to do is create a single scatterplot of the multiple conditions (specified by

i installed red hat 7.0, got teh codeweavers version of wine up and runnin, and i fire up the setup program for ms word 97, and it runs into this error: setup is unable to open the data file X:\~MSSETUP.t\~msstfof.t\Word97.stf';run setup again from were you ariginally ran it im stumped, what should i be doing here differently...thanks ............................... Daniel W. Schar