Displaying 6 results from an estimated 6 matches for "resf".
Did you mean:
ref
2011 Jun 07
1
variable selection in linear regression
...matrix(mydata[,2:p])
y <- as.matrix(mydata[,1])
X <- as.matrix(rep(1,n))
b <- lm(y~1,data=mydata)$coefficients
yhat <- X%*%b
res <- y-yhat
sigma.hat <- sqrt(sum(res^2)/(n-ncol(X)))
cv <- sigma.hat^2*ginv(t(X)%*%X)
se <- sqrt(diag(cv))
pc <- matrix(0,nrow=1,ncol=d)
resF <- matrix(0, nrow=n, ncol=d)
pf <- matrix(0, nrow=1, ncol=d)
for(j in 1:d){
pc[,j] <- cor(x=(x[,j]), y=(mydata[,1]))
resF[,j] <- lsfit(x[,j], y)$residuals
sseF <- t(as.matrix(apply(resF^2, 2, sum)))
resR <- lm(y~1,data=mydata)$residuals
sseR <- sum(resR^2)
dfF <- n-2
dfR...
2013 Feb 01
29
cumulative sum by group and under some criteria
Thank you very much for your reply. Your code work well with this example.
I modified a little to fit my real data, I got an error massage.
Error in split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) :
Group length is 0 but data length > 0
On Thu, Jan 31, 2013 at 12:21 PM, arun kirshna [via R] <
ml-node+s789695n4657196h87@n4.nabble.com> wrote:
> Hi,
> Try this:
>
2013 Feb 28
11
new question
...ta[i]; colnames(x)<-NULL;x}))
indx<-combn(names(resNew),2)
resPval<-do.call(cbind,lapply(seq_len(ncol(indx)),function(i) {x<-as.data.frame(apply(resNew[,indx[,i]],1,t.test.p.value)); colnames(x)<-paste("Pvalue",paste(indx[,i],collapse=""),sep="_");x}))
resF<-cbind(resCounts,resPval)
resF
}
fpv(resCounts)
A.K.
________________________________
From: Vera Costa <veracosta.rt at gmail.com>
To: arun <smartpink111 at yahoo.com>
Sent: Thursday, February 28, 2013 11:30 AM
Subject: new question
Sorry about my question, but I need a ne...
2011 Apr 01
1
Node Recovery locks I/O in two-node OCFS2 cluster (DRBD 8.3.8 / Ubuntu 10.10)
I am running a two-node web cluster on OCFS2 via DRBD Primary/Primary
(v8.3.8) and Pacemaker. Everything seems to be working great, except during
testing of hard-boot scenarios.
Whenever I hard-boot one of the nodes, the other node is successfully fenced
and marked ?Outdated?
* <resource minor="0" cs="WFConnection" ro1="Primary" ro2="Unknown"
2011 Jun 30
2
Saving fExtremes estimates and k-block return level with confidence intervals.
I am estimating a large model by groups. How do you save the results and?returns
the associated quantiles?
For this example I need a data frame
n?? ?xi??????? mu????????beta
1?? 0.1033614? 2.5389580 0.9092611
2? ?0.3401922? 0.5192882 1.5290615
3?? 0.5130798? 0.5668308 1.2105666
I also want to apply gevrlevelPlot() for each "n" or group.
?
#Example
n <- c(1, 1, 1, 1, 1, 1, 2, 2, 2,
2013 Feb 15
10
reading data
Hi,
#working directory data1 #changed name data to data1.? Added some files in each of sub directories a1, a2, etc.
?indx1<- indx[indx!=""]
lapply(indx1,function(x) list.files(x))
#[[1]]
#[1] "a1.txt"??????? "mmmmm11kk.txt"
#[[2]]
#[1] "a2.txt"??????? "mmmmm11kk.txt"
#[[3]]
#[1] "a3.txt"??????? "mmmmm11kk.txt"
#[[4]]
#[1]