search for: se1

...frame(fitted(m1),d1$R),data.frame(fitted(m2),d2$R)) #Doesn't work... Is it correct to use t-test in this situation and how to compute df for it? My solution is: s1<-coefficients(summary(lm(fitted(m1)~d1$R)))[2,1] #Slopes of regr.line s2<-coefficients(summary(lm(fitted(m2)~d2$R)))[2,1] se1<-coefficients(summary(lm(fitted(m1)~d1$R)))[2,2] #SE of slopes se2<-coefficients(summary(lm(fitted(m2)~d1$R)))[2,2] df1<-df.residual(lm(fitted(m1)~d1$R)) #D. of f. df2<-df.residual(lm(fitted(m2)~d1$R)) kk<-function(se1,se2,df1,df2){(se1^2+se2^2)^2/(se1^4/(df1-1)+se2^4/(df2-1))} #D....

I'm experimenting (using for the first time) with using a cisco3620 to connect to the PSTN via a channelised E1 interface, with * handling all of the SIP calls. If anyone has any installation tips / help / documentation I would be most appreciative :) However, my first question is this: when I am in the setup, I see the following: Current interface summary Controller Timeslots

Maybe embarassing, but I can't quickly find out how to do this: Mix expression() text and regular text __in a variable__ for plotting: I can easily do plot(1,type="n",main=expression(paste("Plot of ", x ^ alpha, "; ", alpha==2.3) However, if alpha is really just in a variable, I don't know what to do alpha <- 2.3 ch.a <- format(alpha) and now

Maybe embarassing, but I can't quickly find out how to do this: Mix expression() text and regular text __in a variable__ for plotting: I can easily do plot(1,type="n",main=expression(paste("Plot of ", x ^ alpha, "; ", alpha==2.3) However, if alpha is really just in a variable, I don't know what to do alpha <- 2.3 ch.a <- format(alpha) and now

Guys, I need to make asterisk connect to another device over a pri using QSIG signalling (master/slave). Is this feasible and has it been done before? Thank you! al. _________________________________________________________________ Get ready for school! Find articles, homework help and more in the Back to School Guide! http://special.msn.com/network/04backtoschool.armx

...20 with a proper T1/PRI card installed with asterisk running on the same LAN. Since I have lit up the line, I can dial out and make calls to regular lands lines. However when a call comes back in it rings the destination phone once and disconnects. Here is an error from my router 15:40:45: ISDN Se1/0:23 SERROR: L3_GetUser_NLCB: EVENT 0X45 No NLCB 2 15:40:45: ISDN Se1/0:23 **ERROR**: Ux_BadMsg: Invalid Message for call state 9, call id 0x253, call ref 0x83DF, event 0x62 15:40:45: ISDN Se1/0:23 SERROR: CCPRI_Go: call id 0x254 event 0x57 No ccb Source->HOST 15:40:45: ISDN **ERROR**: Module...

...if I am suppose to hear announcement or fastbusy. Below is captured ISDN debug from Cisco. It's sending ALERT (progress indicator 8), but I am not sure why * is not relaying that signal to the ATA. I hope someone is out there to help me on this. Thanks in advance. Jan 15 22:53:24: ISDN Se1/0:23: TX -> SETUP pd = 8 callref = 0x0078 Jan 15 22:53:24: Sending Complete Jan 15 22:53:24: Bearer Capability i = 0x8090A2 Jan 15 22:53:24: Channel ID i = 0xA98397 Jan 15 22:53:24: Calling Party Number i = 0x0081, '6882332', Plan:Unknown, Type:Unknown J...

Hi, I found some info messages below in /var/spool/mail/root. To reduce the size of that file, I'm wondering if there's a way to prevent those unimportant messages from entering /var/spool/mail/root. Thanks. 2010-05-10T10:32:27-05:00 <authpriv.info> se1 su: pam_unix(su-l:session): session closed for user root 2010-05-10T10:01:28-05:00 <local1.info> se1 Keepalived_vrrp: VRRP_Instance(VI_1) Received lower prio advert, forcing new election 2010-05-10T10:02:02-05:00 <syslog.info> se1 syslog-ng[12531]: Log statistics; dropped='tcp(........

...the value of "r" for the whole AUC's . > seROC<-function(AUC,na,nn){ > a<-AUC > q1<-a/(2-a) > q2<-(2*a^2)/(1+a) > se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) > se > } > > cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ > se1<-seROC(AUC1,na1,nn1) > se2<-seROC(AUC2,na2,nn2) > > sed<-sqrt(se1^2+se2^2-2*r*se1*se2) > zad<-(AUC1-AUC2)/sed > p<-dnorm(zad) > a<-list(zad,p) > a Could somebody kindly suggest me how to calculate the value of "r" or some ways to calculate the stati...

...or your patience and help. EDIT: So I tried the following, and I'm not sure if I'm doing it properly or making the right conclusions. Furthermore, there has to be a more efficient way to perform this task in R. For a given conversion rate (p) and number of trials (n): p1 = 0.1 n1 = 1000 se1 = sqrt( p1 * (1-p1) / n1 ) se1 se1 * 1.96 (p1 + 1.96*se1) * 100 (p1 - 1.96*se1) * 100 p2 = 0.11 n2 = 1300 se2 = sqrt( p2 * (1-p2) / n2 ) se2 se2 * 1.96 (p2 + 1.96*se2) * 100 (p2 - 1.96*se2) * 100 (8.1, 11.8) # headline 1 (9.2, 12.7) # headline 2 # these confidence intervals for the two headlines...

...g program I found on R-help archives : From: Bernardo Rangel Tura Date: Thu 16 Dec 2004 - 07:30:37 EST seROC<-function(AUC,na,nn){ a<-AUC q1<-a/(2-a) q2<-(2*a^2)/(1+a) se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) se } cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ se1<-seROC(AUC1,na1,nn1) se2<-seROC(AUC2,na2,nn2) sed<-sqrt(se1^2+se2^2-2*r*se1*se2) zad<-(AUC1-AUC2)/sed p<-dnorm(zad) a<-list(zad,p) a } The author of this script says: "The first function (seROC) calculate the standard error of ROC curve, the second function (cROC) compare...

...g program I found on R-help archives : From: Bernardo Rangel Tura Date: Thu 16 Dec 2004 - 07:30:37 EST seROC<-function(AUC,na,nn){ a<-AUC q1<-a/(2-a) q2<-(2*a^2)/(1+a) se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) se } cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ se1<-seROC(AUC1,na1,nn1) se2<-seROC(AUC2,na2,nn2) sed<-sqrt(se1^2+se2^2-2*r*se1*se2) zad<-(AUC1-AUC2)/sed p<-dnorm(zad) a<-list(zad,p) a } The author of this script says: "The first function (seROC) calculate the standard error of ROC curve, the second function (cROC) compare...

...tine dgesv: system is exactly singular" It appears this is because the data frame is passing character strings rather than the numerical data to R. An example of my original data frame minmaxfunc is as follows: min max T1 1.500000e+01 3.999954e+01 SE1 0.000000e+00 1.000000e+00 PRE 0.000000e+00 1.000000e+00 WET 0.000000e+00 5.980000e+00 BE1 4.664642e+00 5.866620e+00 Kappa1 5.500000e+03 2.000000e+04 Kappa3 1.000000e+04 2.000000e+04 Then I created a latin hypercube set using (qunif(x[,i], minmaxfunc$min[i], minmaxfunc$max[...

...ngel Tura > Date: Thu 16 Dec 2004 - 07:30:37 EST > > seROC<-function(AUC,na,nn){ > a<-AUC > q1<-a/(2-a) > q2<-(2*a^2)/(1+a) > se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) > se > } > > cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ > se1<-seROC(AUC1,na1,nn1) > se2<-seROC(AUC2,na2,nn2) > > sed<-sqrt(se1^2+se2^2-2*r*se1*se2) > zad<-(AUC1-AUC2)/sed > p<-dnorm(zad) > a<-list(zad,p) > a > } > -- Waverley @ Palo Alto

Dear Rusers, The error for the following was that it cannot allocate the vector of length 1920165909. a <- expand.grid(se1=0:100/100, sp1=0:100/100, se2=0:100/100, sp2=0:100/100, DR=0:100/100) How to solve it? Maybe setwd(dir) can, i am not very sure about it. Any ideas about it? [[alternative HTML version deleted]]

...676.55377 -9545.785 2422.5340 -7183.686 7 -631.91921 10997.432 -3310.2905 -5348.513 Any help would be much appreciated. Thanks, Daniel Crouch Daniel Crouch Research Student Department of Medical & Molecular Genetics King's College London 8th Floor, Tower Wing Guy's Hospital London SE1 9RT United Kingdom

Dear R-helper, I would like to compare the AUC of two logistic regression models (same population). Is it possible with R ? Thank you Roman Rouzier [[alternative HTML version deleted]]

...n produced rather than, say, a plot/table etc shown on a web page. I've had limited success finding examples on this. Many thanks. David Carslaw Environmental Research Group MRC-HPA Centre for Environment and Health King's College London Franklin Wilkins Building Stamford Street London SE1 9NH david.carslaw at kcl.ac.uk -- View this message in context: http://r.789695.n4.nabble.com/Examples-of-web-based-Sweave-use-tp3425324p3425324.html Sent from the R help mailing list archive at Nabble.com.

I just inherited a Cisco 2621 with a VWIC-1MFT-T1 card in it. Can I make this thing into MGCP gateway or even a SIP gateway for asterisk? Seems like it should bee useful for something! I'm perfectly happy to do my homework, but also don't feel thee need to reinvent the wheel! So, links with relevant info would be appreciated. If there is a config for a 2621 being used as a gateway

...member of any queue. Oct 25 10:27:25 NOTICE[3777] chan_zap.c: PRI got event: HDLC Abort (6) on Primary D-channel of span 3 Oct 25 10:27:25 DEBUG[3777] chan_zap.c: Got event HDLC Abort (6) on D-channel for span 3 *The ISDN Messages on the IAD when the call is flashed* *Mar 3 15:47:58 EST: ISDN Se1/0:23 Q931: RX <- DISCONNECT pd = 8 callref = 0x0057 Cause i = 0x8190 - Normal call clearing *Mar 3 15:47:58 EST: %ISDN-6-DISCONNECT: Interface Serial1/0:0 disconnected from 6782363001 , call lasted 31 seconds *Mar 3 15:47:58 EST: ISDN Se1/0:23 Q931: TX -> RELEASE pd = 8 callref =...