search for: zb

This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c("ZB:X","ZC:X") # and test.terms = ~ ZB:X + ZC:X in 'regTermTest' differ? I thought that both wo...

Hi all, I have been previously been using scatterplot3d package to create some graphs but unfortunately it does not allow me to rotate the plot on all three axis. The cloud() function in the lattice package does allow me to do so. When I was using scatterplot3d I was using a script (Shown Below) to calculate the mean, quartiles and range limits for all three axis and I was representing that on the

...t of an nuw add with an nuw > add of zexts changes the behavior of a well-defined program. > Correct me if I am wrong but we are talking about transforming: %maybe_poison = add nuw i32 %a, %b %x = zext i32 %maybe_poison to i64 %y = lshr i64 %x 32 To: %za = zext i32 %a to i64 %zb = zext i32 %b to i64 %x = add nuw i64 %za, %zb %y = lshr i64 %x 32 ? If so, this seems fine in the model given by the RFC. In the before case: %maybe_poison doesn't infect the upper bits of %x with its poison which allows %y to have the well defined result 0. In the after case: %za an...

Hi all, I want to combine two data sets (ZA and ZB to get ZAB). The common variable between the two data sets is ID. Data ZA ID F M 1 0 0 2 0 0 3 1 2 4 1 0 5 3 2 6 5 4 Data ZB ID v1 v2 v3 3 2.5 3.4 302 4 8.6 2.9 317 5 9.7 4.0 325 6 7.5 1.9 296 Output (ZAB) ID F M v1 v2 v3 1 0 0 -9 -9 -9 2 0 0 -9 -9 -9 3 1...

Dear all As a new user of R, can someone please help me with the following I have created a programme to analyse laboratory data and one of the graphs is a bar plot of 'Z' scores. On the bar plot I am using the following line to plot some results barplot (zb[,c("ZBW")], ylim = c(-6,6), names.arg=zb[,c("LabNo")],xlab="Lab Code Number", cex.names = .5 , axis.lty=1 , main = " Total Coliform ") What I would like to do is move the xlab parameter to sit on top or on each bar and not on the x-axis at y=0. Any help wou...

...this is out of my statistical trianing as biologist. Here is my brief description of the problem. I have a survey that VESSELs are selected at random with the probability of p(j). Then the tows within the jth VESSEL are sampled at random with probability of p(i|j). I write my model as y = XB + Zb + e where XB is fixed part, Zb is for random effect (VESSEL) and e is within-vessel error. I feel that I should weight the Zb part by p(j) and the e-part by p(i,j)=p(j)*p(i|j). Is this a correct weighting? How can I implement the weightings in nlme (or lme)? I think that p(i,j) can be specified...

Is it possible to install and run Netscape whitout having access to a native windows. -- Best regard, Svend Bjerrum Jensen http://www.sbj.esenet.dk

I am interested to know how to make for clouds: - aspect ratio = 1 - labels attached to points - vertical lines from the points to the x/y base plane I tried: t = c( 'A', 'B', 'C', 'D' ) x = c( 100, 0, 200, 100 ) y = c( 0, 100, 0, 100 ) z = c( 80, 0, 20, 40 ) q = data.frame( x, y, z ) rownames( q ) = t print(cloud( z ~ x *

...ass = "factor"), x = c(0.3, 0.2, 0.1, 0, 0, 0, 0.2, 0.3)), .Names = c("z", "x"), row.names = c(NA, -8L), class = "data.frame") I would like know if it's possible to use model.matrix() to construct the following design matrix in some sensible way: za zb 1 1 0 2 1 0 3 1 0 4 0 0 5 0 0 6 0 0 7 0 1 8 0 1 In other words, I want column 1 to be (z=="a" & x>0) and column 2 to be (z=="b" & x>0). I can construct this matrix using sweep(model.matrix(~z-1,dd),1,dd$x>0,"*") and then stick it...

I'm new to R, after many years of using Matlab. I've found the R function 'polygon' to be nearly equivalent to the Matlab function 'patch'. For example, the R commands: plot(c(0, 5), c(0, 4), type = 'n', asp = 1, ann = FALSE) x <- c(1, 2, 2, 1.5, 1) z <- c(1, 1, 2, 1.7, 2) polygon(x, z, col = 'green') produce a plot with a small green shape exactly

Hi R Users, I have produced a simulated scatter plot of y versus x tightly clustered around the 45 degree line through the origin with the following code: x <- seq(1,100) y <- x+rnorm(100,0,10) plot(x,y,col="blue") abline(0,1) Is there some way to generate a 3-dimensional analogue of this? Can I get a similar simulated scatter plot of points in 3 dimensions where the points

...ans$ne + loess.ne<-loess(ans$ne ~ pc, span=l.span) + lines(pc,loess.ne$fit,lty=1,col=4) + old.ne<-loess.ne$fit + ###readline() + } + return(ne=ans$ne, Ep=ans$Ep) +                } + } > > ## needed for method 1 > nef2<-function(rc,nc,re,ne,alpha,power){ + za<-qnorm(1-alpha) + zb<-qnorm(power) + xe<-asin(sqrt((re+0.375)/(ne+0.75))) + xc<-asin(sqrt((rc+0.375)/(nc+0.75))) + ans<- 1/(4*(xc-xe)^2/(za+zb)^2-1/(nc+0.5)) - 0.5 + return(ans) + } > ## needed for method 2 > nef<-function(rc,nc,re,ne,alpha,power){ + za<-qnorm(1-alpha) + zb<-qnorm(power) + xe...

Hi folks, I was wondering how to run a mixed models approach to analyze a linear regression with a user-defined covariance structure. I have my model y = xa +zb +e and b ~ N (0, C*sigma_square). (and a is a fixed effects) I would like to provide R the C (variance-covariance) matrix I can easily provide an example, but at this point I am first trying to know what is the best package the allows an unstructured covariance matrix. I was trying the functio...

...f. This option determines whether Shorewall accepts TCP packets which are not part of an established connection and that are not ''SYN'' packets (SYN flag on and ACK flag off). 2. The need for the ''multi'' option to communicate between zones za and zb on the same interface is removed in the case where the chain ''za2zb'' and/or ''zb2za'' exists. ''za2zb'' will exist if: a. There is a policy for za to zb. b. There is at least one rule for za to zb. 3. The /etc/shorewall/blacklist file...

Hello, I'm attempting to set up a daily Robocopy task on a Windows Server 2016 to duplicate a large directory tree -- including data, timestamps and ACLs -- to a Samba 4.5.12 server. This is my Robocopy command: Robocopy.exe C:\source \\foo\bar /XA:SH /XJ /MIR /NODCOPY /COPY:DTSO /A+:A /ZB I noticed that for a significant number of files, Robocopy will systematically identify them as modified even though no aspect of them will have changed in any way. After much trial and error, I tracked the cause down to the archive bit. If the archive bit is set the same way on the source and de...

...FreeBSD. Kris (CCed) observed strange situation. In function arc_read() he had a panic on assertion that we try to unlock a lock which is not beeing held: rw_enter(&pbuf->b_hdr->b_datalock, RW_READER); err = arc_read_nolock(pio, spa, bp, done, private, priority, flags, arc_flags, zb); rw_exit(&pbuf->b_hdr->b_datalock); <--- THIS ONE The only possiblity was that b_hdr for pbuf was changed somewhere. We diagnozed this further and the b_hdr field is changed in arc_release() function, here: buf->b_hdr = nhdr; Backtrace for this change is the following: arc_r...

Hi, I'm a beginner on VPN's World. I think TINC is a solution for me. I want to connect thru Internet 2 different locations, both networks are Windows based. I downloaded TINC Installer for Windows and followed Installation Guide step by step, everything is ok but when I tried to run TINC Deamon a message it's displayed said: "The services has started and finished immediately

...the data to be placed on a regular grid and you need to specify limits, xn <- which(xf>=your_lon_low_limit & xf<=your_lon_high_limit) yn <- which(yf>=your_lat_low_limit & yf<=your_lat_high_limit ) x <- xf[xn] y <- yf[yn] z <- zf[xn,yn] za <- apply(z,1,rev) zb <- t(za) yb <- -rev(y) xb <- x nx <- length(xb) dx <- (xb[nx]-xb[1])/(nx-1) ny <- length(yb) dy <- (yb[ny]-yb[1])/(ny-1) mf <- 5 xo <- seq(xb[1],xb[nx],by=dx/mf) yo <- seq(yb[1],yb[ny],by=dy/mf) nxo <- length(xo) nyo <- length(yo) zo1 <- matrix(NA,nrow=nxo,nc...

...ist. Here is my brief description of the problem. >> >>I have a survey that VESSELs are selected at random with the probability >>of p(j). Then the tows within the jth VESSEL are sampled at random with >>probability of p(i|j). I write my model as >> >>y = XB + Zb + e >>where XB is fixed part, Zb is for random effect (VESSEL) and e is >>within-vessel error. >> >>I feel that I should weight the Zb part by p(j) and the e-part by >>p(i,j)=p(j)*p(i|j). Is this a correct weighting? >> >>How can I implement the weightings...

...e to assume a maximally fragmented system where each +; extent contains only one sector. +; + alignz 2 MaxInitDataSize equ 96 << 10 MaxLMA equ TEXT_START+SECTOR_SIZE+MaxInitDataSize -SectorPtrs times MaxInitDataSize >> SECTOR_SHIFT dd 0 +SectorPtrs zb 10*(MaxInitDataSize >> SECTOR_SHIFT) SectorPtrsEnd equ $ ; ---------------------------------------------------------------------------- As you can see from the subject, I'm overflowing Sector 1 by 5 bytes. Could you explain why you changed 0x3FC to 0x3F6? I admit to not knowing an...