similar to: How to obtain nls parameter estimates

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

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

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

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

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 am running 3.0.10-1.4E on RHEL4. The machine is a ADS member server. I would like to statically map the ADS group "Domain Admins" to the built in "wheel" group so all members of "Domain Admins" are in the "wheel" group. I have looked at the username map option, but I don't want a group of users mapped to a UID (this would defeat what I am trying to

I have setup a SAMBA + OPENLDAP server following Samba-3 By Example, but I have run into several problems. All the tests described in the Chapter 6 such as pdbedit -Lv, slapcat, and ldapsearch -x -b give the desired results. But, getent can't read passwd or group information from ldap backend. I have done my best to solve the problem, but it just doesn't work. I have compiled the nss_ldap

"My Documents" folder doesn't sync with the server unless users log off from the machine and log back in. Is there a way to sync the "My Documents" folder during the session? Thank you. --------------------------------- Do you Yahoo!? Yahoo! Search presents - Jib Jab's 'Second Term'

I am running a 3-parameter nonlinear fit using the default Gauss-Newton method of nls. initialValues.L = list(b=4,d=0.04,t=180); fit.nls.L = nls( myModel.nlm , fData.L, start = initialValues.L, control = nls.control(warnOnly = TRUE), trace=T ); summary.nls.L = summary(fit.nls.L); I run the same analysis in SAS proc NLIN. proc nlin data=apples outest=a; parms b=4 d=.04 t=180; model Y =

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:

Hello, I''ve got Package resource on my Windows machine that keeps creating itself even though the package has been properly installed. I''ve got the name attribute set to the DisplayName found in the Registry for the package. Every time the agent polls, I see in the Event Logs: /Stage[main]/Mssql_rtwo_sptwo/Package[mssql2008r2sp2_exec]/ensure: created This is