search for: 100,5

...example of what I mean: barnames<-c("I agree completely", "I agree", "I partly agree", "I do not agree", "I really hate this stupid question, don't you?") answers<-data.frame(question=paste("Q",1:5,sep=""),S1=sample(1:100,5),S2=sample(1:100,5),S3=sample(1:100,5),S4=sample(1:100,5),S5=sample(1:100,5)) Width<-50 Cex<-1.5 par(mfrow=c(nrow(answers)+1,1),mar=c(0,1,1,1)) plot.new() plot.window(xlim=c(0,1),ylim=c(0,1)) barnames.plot<-do.call("c",lapply(barnames,function(x)paste(strwrap(x,Width),collapse...

...sar el paquete QCC y colocar mas de gráfico de control en una sola ventana usando función "par" con "qcc" para presentar varios gráficos en la misma ventana. He probado pero no funciona, solo muestra un gráfico miren el código # set unequal sample sizes n <- c(rep(50,5), rep(100,5), rep(25, 5)) # generate randomly the number of successes x <- rbinom(length(n), n, 0.2) par(mfrow=c(1,2)) # plot the control chart with variable limits qcc(x, type="p", size=n) # plot the standardized control chart qcc(x, type="np", size=n) Gracias si alguien me ayuda Ate...

Hello, suppose I have a data frame: > mat id age 1 NA NA 2 NA NA 3 NA NA 4 NA NA 5 NA NA Then I attach the data frame: > attach(mat) I assign some new values... > id <<- sample(100,5) > age <<- rnorm(5,mean=30) Then I want to create a new data frame from the variables id and age which still are attached to position 2 of the R environment... > new_mat <- data.frame(ls(pos=2)) # I want to rescue ALL variables that were created by attach(mat) > new_mat ls.p...

Dear R users, I want to group numerical values in classes with different size and count the values for each classes. My data is in this forma: d 15 12,5 30,4 20,5 80,4 100,5 8,2 40,5 33 21 11 And I want the group them in classes with 4 (5,....etc) cm size like this: class d 16 16 32 24 84 104 12 44 36 24 12 and final to count values for each class: class d n 12 2 16 2 24 2 32 1 36 1 44 1 84 1 104 1 Total 11 Than...

All, I have been beating my head against a wall trying to figure out how I would implement a separate moderator code and participant code for the same conference using meetme, the deal is I dont want the participants to be able to join until the moderator is in the conference. Is it possible to do this using the apps as they are , or is their a way to use an Agi script, is that the only way?

Hi, I've recently upgraded from Asterisk 1.4.22 to 1.4.29, and I'm now having a problem with Originate and chan_local. I'm using the following Manager API action to originate a call: Action: originate Priority: 1 Context: trunk Callerid: 100 Channel: Local/100 at callback/n Exten: 123456789 Variable: USERFIELD=127.0.0.1|USEREXT=123456789 WaitTime: 30 This is intended to first call extension 100 in the callback context, and then when that is answered, call 123456789 in the trunk context. I have the following in the callback context:...

Hi! I will calculate sum??s in the following way: E.g.: a <- rpois(100,20) b <- rpois(100,5) x <- data.frame(cbind(a,b)) # the sum??s should be calculated based on a. attach(x) sort.nace <- unique(sort(x[,1])) sum1 <- matrix(ncol=1, nrow=length(sort.nace)) # I calculate the sum of all values of b, which have the same category in a. Eg.: sum1[1,] &l...

...as de gráfico de control > en una sola ventana usando función "par" con "qcc" para presentar varios > gráficos en la misma ventana. > He probado pero no funciona, solo muestra un gráfico miren el código > > # set unequal sample sizes > n <- c(rep(50,5), rep(100,5), rep(25, 5)) > # generate randomly the number of successes > x <- rbinom(length(n), n, 0.2) > par(mfrow=c(1,2)) > # plot the control chart with variable limits > qcc(x, type="p", size=n) > # plot the standardized control chart > qcc(x, type="np", size=...

Hello! I would like to add names to a table. The table presents median values, numbers of values (n) and percentage(%) for a given agegroup. Thus I would like to add a label above these three variables for each category. But if I try to use colnames I get a message telling me that the number of labels must be equal to the number of columns in the table. I hope someone could understand what I mean!

Hi I have two small question, if you can help me ;=) Problems with Answer+Music my extension: [Cal-In] exten => _811XXXX20,1,Goto(C-Internal,100,1) exten => _811XXXX21,1,Goto(C-Internal,200,1) [C-Phibee] exten => 100,1,Ringing exten => 100,2,Wait,1 exten => 100,3,Answer exten => 100,4,Dial(SIP/201&SIP/200,30) exten => 100,5,Hangup exten => 200,1,Ringing exten => 200,2,Wait,1 exten => 200,3,Answer ex...

...trying to get a inverse matrix of (X^T)X , but I keep getting the error message like: no b argument and no default value for sprintf(gettext(fmt, domain = domain), ...) . -------------------------------------------------------------------------------------------- # my code X<-Matrix(rep(1,500),100,5) X[lower.tri(X)]<-1-10^-7 XtX<- t(X)%*% X XtXu<-lu(XtX) InverseXtX<-solve(XtXu) --------------------------------------------------------------------------------------------- The cause suggested for this error is that the lower.tri function returns a matrix of logicals. If so, is try...

Hi, people. Executing the following command: hist(rpois(100,5), labels=TRUE) yields a graphic in which some labels are truncated (on an X11 device). The truncated labels are those over the highest bars. The hist function should ideally manage enough room for the labels, automatically. (Specifying ylim solves my problem, but yet, hist could be frienlier.)...

...reedom were used to obtain the P value (essentially I don't know how to report the result). The following example (while errors are not overdispersed) otherwise generates a comparable analysis of deviance table to my analysis. Any help would be much appreciated. Jonathan counts <- c(rpois(100,5),rpois(100,20)) sites <- rep(100,200) fac1 <- factor(c(rep("A",100),rep("B",100))) fac2 <- factor(c(rep("C",50),rep("D",100),rep("C",50))) model1 <- glm(counts ~ fac1 * fac2,family=quasipoisson, offset=log(sites)) model2 <- glm(coun...

...than 14. Example code is below. I have gon as far as to select markers and calculate the distances, A while loop with set with max(distances)>14 should work, but I'm not sure where to set this. Any help would be appreciated. Rodrigo #some markers, chromosomes and locations loc<- rep(1:100, 5) marker<- paste("A", seq(1:500), sep="") chr<- rep(1:5, rep(100,5)) mm<- data.frame(marker, chr, loc) selectmarkers<- function(n=10){ tapply(mm$marker, mm$chr, function(m){ # while (max(distances) > 14) { mars<- sort(sample(m, n)) distances&l...

...th 2 TE110P cards. 1 card is plugged to telco line, another card is plugged with a Hicom PBX. I want to send some call to VoIP phones and all other to my PBX. I don't known how to make my dialplan : ===========Extensions.conf========== [incoming_call] exten => 090200000,1,Goto(callcenter,100,1) exten => 022956353,1,Goto(callcenter,100,1) exten => 022956388,1,Goto(callcenter,100,1) exten => 022956355,1,Goto(callcenter,101,1) exten => s,1,Goto(go_to_pbx,200,1) [callcenter] exten => 100,1,Answer exten => 100,2,SetMusicOnHold(default) exten => 100,3,DigitTimeout,5 ext...

...ncluded in every regression. Then i loop over the columns of D to use all unique combinations of covariates in that matrix and save the results for variable "x" in each run. This is the code : (due to the random numbers it will produce 10 similar betas,t.values etc...) M<-matrix(rnorm(100),100,3) D<-matrix(rnorm(100),100,5) y<-matrix(rnorm(100),100,1) x<-matrix(rnorm(100),100,1) beta<-NULL t.value<-NULL sd<-NULL i<-1 while(i<=ncol(D)){   j<-i+1   while(j<=ncol(D)){     k<-j+1     while(k<=ncol(D)){       Res<-lm(y~1+x+M+ D[,i] + D[,j] + D[,k])...

...[local] ignorepat => 9 include => default include => parkedcalls [default] exten => 001,1,Dial(SIP/001,5) exten => 001,n,VoiceMail(001 at default) exten => 001,n,Playback(vm-goodbye) exten => 001,n,Hangup exten => 901,1,Dial(IAX2/901,5) exten => 901,n,Hangup exten => 100,1,Dial(SIP/100,5) exten => 100,n,VoiceMail(100 at default) exten => 100,n,Playback(vm-goodbye) exten => 100,n,Hangup exten => _0XXXXXXXXXX,1,Dial(Zap/g0/${EXTEN}) exten => _ZXXXXX,1,Dial(Zap/g0/${EXTEN}) [from-pstn] ; BT pass in 6 digits only? exten => 032685,1,Answer() exten =&...

...;m having trouble matching the by-hand version with the R code. I have fitted a Cox model - let's use the ovarian data as an example: library(survival) data(ovarian) ova_mod <- coxph(Surv(futime,fustat)~age+rx,data=ovarian) If I want to make predict survival for a new set of individuals at 100 days then that is trivial using predict.coxph e.g.: newdata <- data.frame(futime=rep(100,5),fustat=rep(1,5),age=c(45,50,55,60,65),rx=c(1,2,1,2,1)) preds <- predict(ova_mod,newdata,type="expected") survpreds <- exp(-1*preds) survpreds [1] 0.9967635 0.9969739 0.9859543 0.9868629 0....

...ns where double zeros do not occur: --- R-1.2.2/src/appl/distance.c Sun Oct 15 18:13:25 2000 +++ R-work/src/appl/distance.c Mon Mar 5 10:16:53 2001 @@ -93,5 +93,5 @@ double R_canberra(double *x, int nr, int nc, int i1, int i2) { - double dist; + double dist, sum; int count, j; @@ -100,5 +100,7 @@ for(j=0 ; j<nc ; j++) { if(R_FINITE(x[i1]) && R_FINITE(x[i2])) { - dist += fabs(x[i1] - x[i2])/fabs(x[i1] + x[i2]); + sum = fabs(x[i1] + x[i2]); + if (sum > 0.0) + dist += fabs(x[i1] - x[i2])/sum; count++; } Best wishes, Jari Oksanen -- Jari...

...ns where double zeros do not occur: --- R-1.2.2/src/appl/distance.c Sun Oct 15 18:13:25 2000 +++ R-work/src/appl/distance.c Mon Mar 5 10:16:53 2001 @@ -93,5 +93,5 @@ double R_canberra(double *x, int nr, int nc, int i1, int i2) { - double dist; + double dist, sum; int count, j; @@ -100,5 +100,7 @@ for(j=0 ; j<nc ; j++) { if(R_FINITE(x[i1]) && R_FINITE(x[i2])) { - dist += fabs(x[i1] - x[i2])/fabs(x[i1] + x[i2]); + sum = fabs(x[i1] + x[i2]); + if (sum > 0.0) + dist += fabs(x[i1] - x[i2])/sum; count++; } Best wishes, Jari Oksanen -- Jari...