search for: 10,11

...ug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### In version 0.99.0 for Windows NT I get the following: > .Random.seed [1] 0 21012 4461 28299 > rbinom(10,10,.3) [1] 4 5 5 1 5 2 4 2 6 2 > rbinom(10,9,.3) [1] 0 0 0 0 0 0 0 0 0 0 > rbinom(10,10,.3) [1] 3 3 2 2 3 1 3 3 3 2 > rbinom(10,11,.3) [1] 0 0 0 0 0 0 0 0 0 0 I tried similar inputs in version 0.65.1: > .Random.seed [1] 0 16049 70 30225 > rb...

...s only show the data in 'y'. I'm wondering how to show the data in 'x' as well. I also want to add a legend to show that blue points corresponds to 'x' and yellow points correspond to 'y'. Could somebody let me know what the correct commands should be? x=rbind(c(10,11),c(10,11)) y=cbind(-1:0,-1:0) plot(y,col='yellow') points(x,col='blue')

Hi, I have a little problem with dc10, dc11. I use three quad dc cards, so far from dc0 up to dc8 with no problems. All (dc0 to dc11) are displayed correctly with pciconf and with ifconfig. The trouble is with dc10 and dc11 that they don't send any data out and also don't react to arp requests etc. - at least using tcpdump won...

...he missing data issue as follows: tmtA <- c(20,26, 16,29,22,NA) tmtB <- c(12,24,17,21,NA,17) require(lme4) dv <- c(20,12,26,24,16,17,29,21,22,17) subject <- rep(c("s1","s2","s3","s4","s5","s6"),each=2) subject <- subject[-c(10,11)] myfactor <- rep(c("f1","f2"), 6) myfactor <- myfactor[-c(10,11)] mydata <- data.frame(dv, subject, myfactor) am2 <- lmer(dv ~ myfactor + (1|subject)), data = mydata) summary(am2) anova(am2) subject <- subject[-c(10,11)] Any help would be greatly appreciated...

...om Gelman, Carlin, Stern, and Rubin's Bayesian Data Analysis solutions, problem 3.5. I was trying to modify this code in post.a2 when I ran into this error. post.a1 <- function(mu,sd,y){ ldens <- 0 for (i in 1:length(y)) ldens <- ldens + log(dnorm(y[i],mu,sd)) ldens} y <- c(10,10,12,11,9) mugrid <- c(10,11) sdgrid <- c(1,1.2) logdens1 <- outer (mugrid, sdgrid, post.a1, y) #*** no error messages *** post.a2 <- function(mu,sd,y) { ldens <- sum(log(dnorm(y,mu,sd))) ldens } y <- c(10,10,12,11,9) mugrid <- c(10,11) sdgrid <- c(1,1.2) logdens2 &lt...

Hello, I need to extract sequences from an ordered vector. For example, if a<-c(1,2,3,6,10,11,13) I need to get the followings 4 vectors (1,2,3),(6),(10,11),(13) Thank You -Bruno ------------------------------------------------------ Leggi GRATIS le tue mail con il telefonino i-mode? di Wind http://i-mode.wind.it/

...le alpha testing 2.8. The resulting function for splinefun gives incorrect deriv values when x is= less than the smallest x-value used to create the function (at least in on= e circumstance), but does the correct thing for x greater than the largest = x-value. Here is an example: > x <- 1:10 > y <- sin(x) > > splfun <- splinefun(x,y, method=3D'natural') > > # these should be linear (and are) > splfun( seq(0,1, 0.1) ) [1] 0.5682923 0.5956102 0.6229280 0.6502459 0.6775638 0.7048816 0.7321995 [8] 0.7595174 0.7868352 0.8141531 0.8414710 > > # these s...

...liance, probably pointing to a non-existent directory, and everything goes downhill from there. Avoid this by making sure that /tmp and /var/tmp are real directories. --- appliance/init | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/appliance/init b/appliance/init index 3c5ef1b..3973f18 100755 --- a/appliance/init +++ b/appliance/init @@ -10,6 +10,11 @@ RUNLEVEL=S PREVLEVEL=N export RUNLEVEL PREVLEVEL +# Make sure /tmp /var/tmp are real directories, not symlinks. +rm -f /tmp /var/tmp +mkdir /tmp /var/tmp +chmod 1777 /tmp /var/tmp + # Make sure to find all the libraries, also tho...

...nd calls so that I can take appropriate action in the dial plan. "dahdi show channel 1" shows "Fax Handled: no". Does that mean that I don't have it configured correctly? [channels] ; Span 1: WCTE2/0/1 "WCTE23X (PCI) Card 0 Span 1" (MASTER) ESF/B8ZS RED group=10,11 context=from-pstn switchtype = national signalling = pri_cpe faxdetect=incoming faxdetect_timeout=0 faxbuffers => 12,half channel => 1-23 ; Span 2: WCTE2/0/2 "WCTE23X (PCI) Card 0 Span 2" ESF/B8ZS RED group=10,12 context=from-pstn switchtype = national signalling = pri_cpe chann...

i read a book from WOOD, there's an example which is talking about the pollutant. library(gamair) library(mgcv) y<-gam(death~s(time,bs="cr",k=200)+s(pm10median,bs="cr")+s(so2median,bs="cr")+s(o3median,bs="cr")+s(tmpd,bs="cr"),data=chicago,family=Possion) lag.sum<-function(a,10,11) {n<-length(a) b<-rep(0,n-11) for(i in 0:(11-10)) b<-b+a[(i+1):(n-11+i)] b} death<-chicago$death[4:5114] time<-...

Hi, I'm experiencing difficulties to get a B410P running with Asterisk 10.3.1 and DAHDI 2.6.1. Am I supposed to use DAHDI for this card and ISDN BRI for my country (Belgium)? thx, BC

On 08/05/2010 06:46 AM, Wenbin Zhang wrote: > Hi all, > I'm using postDominatorTree to do some program analysis. My code works > well for small tests, but when I run it on real applications, the > following error occurs: > /Inorder PostDominator Tree: DFSNumbers invalid: 0 slow queries. >...

Hi Andrew, I don't understand why 2 systems running the exact same version of Samba have different behaviour. Is this an option I can disable? regards, John On 19/04/16 11:29, Andrew Bartlett wrote: > On Tue, 2016-04-19 at 10:29 +1000, John Gardeniers wrote: >> I'm setting up a test domain in order to try out Sudoers LDAP and >> have >> run into a problem that has my puzzled. On our production domain I >> can >> run a query such as: >> >> ldapsearch -LLL -p389 -h DC -u me a...

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m<-matrix(0,10,11) m[1,1]<--1 m[1,2]<-1 m[2,1]<--1 m[2,3]<-1 m[3,1]<--1 m[3,4]<-1 m[4,1]<--1 m[4,5]<-1 m[5,1]<--1 m[5,6]<-1 m[6,1]<--1 m[6,7]<-1 m[7,1]<--1 m[7,8]<-1 m[8,1]<--1 m[8,9]<-1...

...know how to use the text() function for regular y~x plots but in this particular case, I'm kinda lost... Here is my code: # creating labels for x-axis labelling lb <- c() for( n in names(o) ) lb <- c(lb, paste("\n",n,"\n",o[[n]], sep="")) boxplot(x[,-c(10,11)], varwidth=TRUE, col=c("#FF0000","#EE0000","#DD0000","#00FF00","#0000FF","#0000DD","#0000BB","#000099"), col.lab="dark red", col.axis="blue", main=&q...

Hi, I try to use write.xls from dataframes2xls inside a function. The write.xls work normally in console, but inside a function it dont find the object. Look this example: > library(dataframes2xls) > test <- function(x){ + a <- data.frame(A=c(1,2),B=c(10,11)) + write.xls(a,file="a.xls") + } > test() Erro em get(s[i]) : objeto 'a' não encontrado Why the write.xls dont find the object a inside a function? Thanks Ronaldo -- 9ª lei - Se você ficaria infeliz ao perder seus dados, faça um back-up (cópia) permanente d...

** Symptom: plot(1:1000,1:1000,xlim=c(10,11),type="l") makes a plot with lots of (nearly) horizontal lines, not just the intended one from (10,10) to (11,11). So do lots of similar commands where the xlim or ylim is much shorter than the range of the data. ** R version: _ platform...

hi. I have a list like x <- list(1:10,11:20,21:30) It's a sort of a 3 x 10 matrix in list form. I would like to reduce the dimension of this list. it would be something like list(1:3, 11:13, 21,23) I tried x[,1:3] does not work of course. Neither lapply(x, [1:3]) works... Any suggestions? ishida [[alternative HTML versi...

...n(), svd() and qr() can all do it. Andy > From: Jonathan Baron > > On 06/17/04 19:04, Fred wrote: > >Dear R-listers: > > > >I am trying to test an algorithm on a set of linearly > independent vectors > >{x1,x2,...,xn}. > > Well, here's an idea, for 10 vectors of length 10, > as columns of a matrix m1. The 11th seems to be needed. > > m1 <- matrix(rnorm(110),10,11) > for (i in 2:11) { > m1[,i] <- resid(lm(m1[,i] ~ m1[, 1:(i-1)])) > } > > Test it with cor(m1[,-11]) > > I'm sure there are better ways. &g...

It took me some time to find this bug in my code. Is this a feature of R? Am I doing something wrong? > a=data.frame(x=1:10,y=1:10) > b=data.frame(x=11:20,y=11:20) > z=data.frame(1:10,11:20) > a$z=z > b$z=z > rbind(a,b) Error in `row.names<-.data.frame`(`*tmp*`, value = c("1", "2", "3", "4", : duplicate 'row.names' are not allowed In additio...