search for: pdata

Just spoke to Jerry, we'll be doing a 3.0.27a on Monday to fix a regression that broke smbfs mounts from Linux. Sorry for the problem. Jeremy.

I have matrix data data<-matrix(cbind(c(0,0,0),c(NA,0,1),c(1,1,1),c(0,1,1)),ncol=3) and I want to create a new matrix by checking each element of the data and put value 0 if i have NA and 1 otherwise. For this reason i made the function below pdata<-matrix(NA,ncol=ncol(data),nrow=nrow(data)) pdata<-sapply(1:nrow(data),function(i) sapply (1:ncol(data),function(j) if (is.na(data[i,j])) {pdata[i,j]<-0} else {pdata[i,j]<-1} )) pdata<-matrix(t(pdata),ncol=ncol(data),nrow=nrow(data)) with pdata > pdata [,1] [,2] [,3] [1,]...

...et it run for about three minutes, thinking it was doing some heavy decoding, but I was wrong. Once this is fixed I'll release this player. Thanks for the help! pSOUNDDATA OGGDATA::DecodeOGGLow() { int StreamID; unsigned int Size, Read; unsigned int Location, numStreams; pSOUNDDATA pData; if(pFile == NULL) return NULL; if((pData = (pSOUNDDATA)malloc(sizeof(sSOUNDDATA))) == NULL) return NULL; pData->Length = ov_pcm_total(&sOVFile, -1); if((pData->Buffer = (char*)malloc(pData->Length)) == NULL) { free(pData); return NULL; } pOVInfo =...

Am I missing a Package? I'm not sure why is won't read the functions. Any help is much appreciated. > PData = Data[,3:10] > Spec = portfolioSpec() Error: could not find function "portfolioSpec" > setTargetReturn(Spec) = mean(colMeans(PData)) Error in setTargetReturn(Spec) = mean(colMeans(PData)) : object 'Spec' not found > Constraints = "LongOnly" > efficientPo...

...r's thesis i need to simulate a panel data with the fixed effects correlated with the predicor, so i run the the following code: set.seed(1970) #######################Panel data simulation with alphai correlated with xi##################################### n <- 5 t <- 4 nt <- n*t pData <- data.frame(id = rep(paste("JohnDoe", 1:n, sep = "_"), each = t),time = rep(1981:1984, n)) rho <-0.95 alphai <- rnorm(n,mean=0,sd=1)#alphai simulation x<- as.matrix(rnorm(nt,1))#xi simulation akro <- kronecker(alphai ,matrix(1,t,1))#kronecker of alphai cormat...

...r's thesis i need to simulate a panel data with the fixed effects correlated with the predicor, so i run the the following code: set.seed(1970) #######################Panel data simulation with alphai correlated with xi##################################### n <- 5 t <- 4 nt <- n*t pData <- data.frame(id = rep(paste("JohnDoe", 1:n, sep = "_"), each = t),time = rep(1981:1984, n)) rho <-0.95 alphai <- rnorm(n,mean=0,sd=1)#alphai simulation x<- as.matrix(rnorm(nt,1))#xi simulation akro <- kronecker(alphai ,matrix(1,t,1))#kronecker of alphai cormat...

...master's thesis i need to simulate a panel data with the fixed effects correlated with the predicor, so i run the the following code: set.seed(1970) #######################Panel data simulation with alphai correlated with xi##################################### n <- 5 t <- 4 nt <- n*t pData <- data.frame(id = rep(paste("JohnDoe", 1:n, sep = "_"), each = t),time = rep(1981:1984, n)) rho <-0.95 alphai <- rnorm(n,mean=0,sd=1)#alphai simulation x<- as.matrix(rnorm(nt,1))#xi simulation akro <- kronecker(alphai ,matrix(1,t,1))#kronecker of alphai cormat&l...

Hello, I am trying to use the plm package for panel econometrics. I am just trying to get started and load my data. It seems from most of the sample documentation that I need to use the pdata.frame function to get my data loaded. However, even after installing the "plm" package, my R installation cannot find the function. I am trying to follow the example in plmEN.pdf ( cran.mirroring.de/doc/vignettes/plm/plmEN.pdf ) > library(plm) Loading required package: kinship Loading...

...hesis i need to simulate a panel data with the fixed effects correlated with the predicor, so i run the the following code: set.seed(1970) #######################Panel data simulation with alphai correlated with xi##################################### n <- 5 t <- 4 nt <- n*t pData <- data.frame(id = rep(paste("JohnDoe", 1:n, sep = "_"), each = t),time = rep(1981:1984, n)) rho <-0.95 alphai <- rnorm(n,mean=0,sd=1)#alphai simulation x<- as.matrix(rnorm(nt,1))#xi simulation akro <- kronecker(alphai ,matrix(1,t,1))#kronecker of alphai...

...l work up some sample data and code if no solutions are found. ...snip > # Get info on first and last pair sets and lab names in last pair > nlabs <- labdata[["LABNUMBER"]] #e.g. 1, 2, 5, 6, ... > nolabs <- length(nlabs) #Total number of labs in last pair of labdata > #dpdata <- labdata[which(labdata[["LABNUMBER"]] == 13),] #Dummy paired data row > dpdata <- rbind(labdata[0,], NA) > pdata <- dpdata > # Loop through pairs of labdata to build pdata for pairs (last 10 years or less) > k=0 > for(j in nlabs){ + for(i in ifpair:ilpair){ +...

...fclose(pFile); return false; case OV_EBADHEADER: fclose(pFile); return false; case OV_EFAULT: fclose(pFile); return false; } return true; } pSOUNDDATA OGGDATA::DecodeOGGLow() { int StreamID; unsigned int Size, Read, Location; pSOUNDDATA pData; if((pFile == NULL) || (pOVFile == NULL)) return NULL; if((pData = (pSOUNDDATA)malloc(sizeof(sSOUNDDATA))) == NULL) return NULL; Size = ov_raw_total(pOVFile, -1); if((pData->Buffer = (char*)malloc(Size)) == NULL) { free(pData); return NULL; } StreamID = 0;...

Hi all, I installed all required packages for plm, but pdata.frame is not a recognized function on R to start with. I followed the first page of "introduction to plm" by Croissant&Millo, Here are the messages: * > library(plm) Loading required package: kinship Loading required package: survival Loading required package: splines Loading requ...

...ctors_per_unit = bsize/bytes_per_sector; + + svfs.BlockSize = bsize; + svfs.TotalBlocks = dsize; + svfs.BlocksAvail = dfree; + svfs.UserBlocksAvail = dfree; + svfs.TotalFileNodes = 0; + svfs.FreeFileNodes = 0; + } + } + data_len = 56; SIVAL(pdata,0,svfs.OptimalTransferSize); SIVAL(pdata,4,svfs.BlockSize); SBIG_UINT(pdata,8,svfs.TotalBlocks); SBIG_UINT(pdata,16,svfs.BlocksAvail); SBIG_UINT(pdata,24,svfs.UserBlocksAvail); SBIG_UINT(pdata,32,svfs.TotalFileNodes); SBIG_UINT(pdata,40,svfs.FreeFileNodes);...

....cc 2005-03-31 13:55:11 -05:00 @@ -51,7 +51,7 @@ if (req->state == STATE_IOREQ_READY) { req->state = STATE_IOREQ_INPROCESS; } else { - BX_INFO(("False I/O requrest ... in-service already: %lx, pvalid: %lx,port: %lx, data: %lx, count: %lx, size: %lx\n", req->state, req->pdata_valid, req->addr, req->u.data, req->count, req->size)); + BX_INFO(("False I/O request ... in-service already: %lx, pvalid: %lx,port: %lx, data: %lx, count: %lx, size: %lx\n", req->state, req->pdata_valid, req->addr, req->u.data, req->count, req->size)); r...

...ode and output for your kind reference. data file is also attached rm(list=ls()) library(MASS) library(bdsmatrix) library(zoo) library(nlme) library(sandwich) library(car) library(lmtest) library(plm) data1<- read.csv(file.choose(),header=TRUE,sep=",") D<-na.omit(data1) attach(D) Pdata<-plm.data(D, index=c("CT","T")) pool11<- plm(Y~X1+X2,data = Pdata, model="pooling") # pool21<- plm(Equity.dividend.1~ Profit.after.tax.1+ LaggedDivd.1+ log(Size.1)+factor(CompanyName)-1,data = Pdata, model="pooling") fixed.mod1<- plm(Y~X1+X2,data...

...tss) } ## prepare the data set.seed(1) n <- 2000 T <- 15 x=rnorm(T*n) x1=runif(T*n) fe=rep(rnorm(1*n),each=T) id=rep(1:(1*n),each=T) ti=rep(1:T,1*n) e=rnorm(T*n) y=x+x1+fe+e subs <- as.logical(rbinom(T*n, 1, .6)); sum(subs) data <- data.frame(y,x,x1,id,ti) dim(data) #[1] 30000 5 pdata <- plm.data(data[subs , ], index = c("id", "ti")) dim(pdata) #[1] 18018 5 paste("n=", length(unique(pdata$id))); paste("T=", length(unique(pdata$ti))); paste("N=", length((pdata$ti))) #[1] "n= 2000" #[1] "T= 15" #[1] &quo...

...> rm(list=ls()) > library(MASS) > library(bdsmatrix) > library(zoo) > library(nlme) > library(sandwich) > library(car) > library(lmtest) > library(plm) > > data1<- read.csv(file.choose(),header=TRUE,sep=",") > D<-na.omit(data1) > attach(D) > Pdata<-plm.data(D, index=c("CT","T")) > > pool11<- plm(Y~X1+X2,data = Pdata, model="pooling") > # pool21<- plm(Equity.dividend.1~ Profit.after.tax.1+ LaggedDivd.1+ > log(Size.1)+factor(CompanyName)-1,data = Pdata, model="pooling") > fixed.mo...

...comparisons were confused without being obvious as to why or where. Our solution: specify that filename is as.character so assignment of file to target is correct(after correcting $Filename) now that using read.AnnotatedDataFrame rather than readphenoData. Data<-ReadAffy(filenames=as.character(pData(pd)$Filename),phenoData=pd) Hurrah! It may be beneficial to others, that if the filename argument isn't specified, that filenames are read from the phenoData object if included here. Thanks! -----Original Message----- From: Martin Morgan [mailto:mtmorgan at fhcrc.org] Sent: Thursday, 26 Ju...

...r? I was thinking that dir() and file.info() with isdir==TRUE being something that might work. These folders or directories are numerically named with no dot extension names or other characters. Typically, these are 3132, 3334, ... Here is what I tried. The last line is where I need more work. pData="C:/Myfiles/R/Data/" setwd (pData) x <- dir(pData, all.files=F) y <- file.info(x, x$isdir==T) TIA Kenneth Ray Hobson, P.E. Oklahoma DOT - QA & IAS Manager 200 N.E. 21st Street Oklahoma City, OK 73105-3204

How do you calculate the percentage of data within 2SD, 3SD, 4SD, 5SD, and 6SD of the mean? I used the following link as the data I'm working with: nb10 <- read.table("http://www.adjoint-functors.net/su/web/314/R/NB10") if this helps answer my question. Can you please explain how to calculate the SD's? Please be specific in which part of the function changes when