Displaying 20 results from an estimated 1000 matches similar to: "weights vs. offset (negative binomial regression)"
2003 Jul 08
2
Can anybody help me on this?
Hi there:
I have this configuration:
|-----[Server 2]
|
[Internet]--------[Router]----------[Switch]------------ [Server 1]
|
|-----[PC1]
|
|-----[PC2]
|
|-----[PC3]
Server 1 has IP 216.251.XXX.XX1
Server 2 has IP 216.251.XXX.XX2
PC1 has IP 216.251.XXX.XX3
PC2 has IP 192.168.XXX.1
PC3 has IP 192.168.XXX.2
How do I configure shorewall in SERVER 2 to block to/from the Internet Port
22
2013 Jan 10
2
transparency in segments()
Dear all,
i would like to plot each value from my datasets as segment with defined
transparency
However, I didnt find out how to set the transparency.
definition by "col=" in par() or segments() doesnt seem to work
any suggestions?
Thanks in advance.
Kind regards,
Robert Pazur
example code:
xx2 <-read.table("http://www.scandinavia.sk/data/R/0_05.csv", sep=";",
2007 May 11
1
model seleciton by leave-one-out cross-validation
Hi, all
When I am using mle.cv(wle), I find a interesting problem: I can't do
leave-one-out cross-validation with mle.cv(wle). I will illustrate the
problem as following:
> xx=matrix(rnorm(20*3),ncol=3)
> bb=c(1,2,0)
> yy=xx%*%bb+rnorm(20,0,0.001)+0
> summary(mle.cv(yy~xx,split=nrow(xx)-1,monte.carlo=2*nrow(xx),verbose=T),
num.max=1)[[1]]
mle.cv: dimension of the split subsample
2009 Mar 28
1
Error in R??
Can someone explain why I am getting the following error: in the r code
below?
Error in solve.default(diag(2) + ((1/currvar) * (XX1 %*% t(XX1)))) :
system is computationally singular: reciprocal condition number = 0
In addition: There were 50 or more warnings (use warnings() to see the first
50)
The R code is part of a bigger program.
##sample from full conditional
2011 Nov 17
1
Vectorizing for weighted distance
Hi All,
I am trying to convert the following piece of matlab code to R:
XX1 = sum(w(:,ones(1,N1)).*X1.*X1,1); #square the elements of X1,
weight it and repeat this vector N1 times
XX2 = sum(w(:,ones(1,N2)).*X2.*X2,1); #square the elements of X2,
weigh and repeat this vector N2 times
X1X2 = (w(:,ones(1,N1)).*X1)'*X2; #get the weighted
'covariance'
2008 Oct 01
0
xpred.rpart() in library(mvpart)
R-users
E-mail: r-help@r-project.org
Hi! R-users.
http://finzi.psych.upenn.edu/R/library/mvpart/html/xpred.rpart.html
says:
data(car.test.frame)
fit <- rpart(Mileage ~ Weight, car.test.frame)
xmat <- xpred.rpart(fit)
xerr <- (xmat - car.test.frame$Mileage)^2
apply(xerr, 2, sum) # cross-validated error estimate
# approx same result as rel. error from printcp(fit)
apply(xerr, 2,
1999 Feb 16
1
Missing tick marks bug on alpha solved
On some systems (alpha), tick marks don't appear on plots. The easiest
way to see the problem is something like:
> plot(0:1,axes=FALSE)
> axis(1,1:2)
The problem is in X11_Line(...) from .../src/unix/devX11.c, which is
so short I've included the whole function below:
static void X11_Line(double x1, double y1, double x2, double y2,
int coords, DevDesc *dd)
{
2010 Aug 17
3
predict.lm, matrix in formula and newdata
Dear all,
I am stumped at what should be a painfully easy task: predicting from an lm object. A toy example would be this:
XX <- matrix(runif(8),ncol=2)
yy <- runif(4)
model <- lm(yy~XX)
XX.pred <- data.frame(matrix(runif(6),ncol=2))
colnames(XX.pred) <- c("XX1","XX2")
predict(model,newdata=XX.pred)
I would have expected the last line to give me the
2007 Oct 14
1
iPhone 1.1.1 problems.
I'm having problems with the iPhone client in the 1.1.1 version of the
iPhone software. Things with other clients work fine, and the iPhone
worked fine before it was upgraded to 1.1.1. Other phones on 1.1.1
failed. I'm still working on testing against the client with the
iPhone 1.0.2 software on it.
The problem is that the client simply doesn't connect to dovecot. I
can see the
2006 Feb 02
0
Sip - no peer or user found on incoming call
Hi list,
I try to connect to a GW which have one domain eg sip.mydomain.com and
have few IPs related to this domain. I register * to this domain with
host=sip.mydomain.com and type=user. So DNS will decide on which IP of
my domain I will register (or redirection on the GW side).
If an incoming call arrive, I would guess that, as type=user, it will
not try to match the IP from INVITE as I
2009 Oct 05
1
gsub - replace multiple occurences with different strings
Hi,
I search a way to replace multiple occurrences of a string with
different strings
depending on the place where it occurs.
I tried the following;
x <- c("xx y e d xx e t f xx e f xx")
x <- gsub("xx", c("x1", "x2", "x3", "x4"), x)
what I want to get is;
x =
x1 y y e d x2 e t f x3 e f x4
but what I get is;
x =
x1 y y e
2007 Dec 02
2
Optimised qmf_synth and iir_mem16
Hi all,
I've taken preglows ARM versions of qmf_synth and iir_mem16 from
rockboxes speex codec, and tweaked them a bit further for some more
speed.
I attach them here so you can review and take on any changes you
want.
Please let me know if you have questions etc.
Thanks,
Robin
--
Robin Watts, Email: <mailto:Robin.Watts@wss.co.uk>
Warm Silence Software, WWW:
2008 Dec 16
8
sliding window over a large vector
Hi all,
I have a very large binary vector, I wish to calculate the number of
1's over sliding windows.
this is my very slow function
slide<-function(seq,window){
n<-length(seq)-window
tot<-c()
tot[1]<-sum(seq[1:window])
for (i in 2:n) {
tot[i]<- tot[i-1]-seq[i-1]+seq[i]
}
return(tot)
}
this works well for for reasonably sized vectors. Does
2008 Mar 04
1
Difficulty labeling plot axis when "mar" and "oma" are set
Dear all,
Though labeling the x and y axis in the plot command seems to be straight
forward, I can not get it to work if I do the following:
## Creating example data
edata <- c(1,2,1,2)
edata <- matrix(edata, 2, 2, byrow = T)
colnames(edata) <- c("a", "b")
edata <- data.frame(edata)
## plot data
par(mfrow<- c(1,2))
plot(edata$a, edata$b, xlab =
2007 Jan 30
2
dev IFB, few questions
I''ve made some tests...
eth2 is my internal interface, LAN is connected here.
Before I had IMQ device in AB mode...
PREROUTING [A]fter NAT, POSTROUTING [B]efore NAT.
I want the same situation on ifb.
I do this in this way:
---
# incoming traffic here from LAN is before NAT
tc qdisc add dev eth2 handle ffff: ingress
# outcoming traffic here from WAN is after NAT
tc qdisc add dev eth2
2002 Jun 13
2
Samba versus Dreamweaver
Hi!
I hope someone can help me with this one:
We're using Macromedia Dreamweaver 3 to publish local files from a NT4
Server to our internet server running Redhat 7.3 / Samba 2.2.3a. Life could
be so sweet but we're facing a strange problem: Users can't _overwrite_
files edited by other users. Everytime someone tries to overwrite such a
file the message "An error occurred -
2010 Jun 15
4
Unable to pickup an extension, tryi
Hi!
> How to do this ??
> To proceed with your answer on PICKUPMARK, where do I put this ???
Look at the example for Asterisk 1.4 on this page:
http://www.voip-info.org/wiki/view/Asterisk+cmd+Pickup
Philipp
2007 Nov 01
2
ploting a comparison of two scores, including the labels in the plot
Hello r-help!
I have data with two kind of ratings on status of 100 occupations. The
first kind of rating is on the percieved "objective" status that these
occupations have in society at large, and the second kind or rating is
on the status that the respondents think that these occuption *should*
have.
The ratings were originally integer values in the rage 1-9, but in the
current data,
2008 May 19
1
Select certain elements from dataframe
Hello,
I have a specific problem, I have a large dataframe, and after clustering I
want to select certain colums, the elements of a subcluster.
My dataframe looks like this :
> colnames(data)
[1] "101KF4319097339" "102KF4319101170" "103KF4319047549"
"104KF4319046389"
[5] "105KF4319013260" "106KF4319025582"
2011 Jan 30
1
Extract subsets of different and unknown lengths from huge dataset
Dear prospective reader,
I apologize for posting my problem but I've just no idea how to go on by
processing this huge (over 70 MB) dataset. Thank you in advance for any help
or comment! I do appreciate it!
My textfile contains 1 column of interest (numbers/values only). The overall
issue is to extract 'events', starting points of which are defined by at
least 24 preceding values