search for: datac

Greetings List, I'm having some trouble with the use of the names function in a formula. Below is an example of something that works (i.e first line), and the second line is the same formula which doesn't. I want to be able to reference the column of the dataC table so I can run the variogram iteratively with a loop. > v<-variogram(A1~1,loc=~x+y, dataC) > v<-variogram(names(dataC[3])~1,loc=~x+y, dataC) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type Where dataC looks like:...

Hi R-helpers, I have been struggling with calculating row and column statistics, e.g. standard deviation. I know that > datac$Mean<-rowMeans(datac,na.rm=TRUE) will give me row means. I have tried to replicate those row means with the apply function: > datac$Mean2<-apply(datac,2,mean) so that I can replace the function argument with "sd" (instead of mean) to get standard deviations. But, I'm runni...

...A's: if na.rm==T, don't count them length2 <- function (x, na.rm=FALSE) { if (na.rm) sum(!is.na(x)) else length(x) } # Collapse the data formula <- as.formula(paste(measurevar, paste(groupvars, collapse=" + "), sep=" ~ ")) datac <- summaryBy(formula, data=data, FUN=c(length2,mean,sd), na.rm=na.rm) # Rename columns names(datac)[ names(datac) == paste(measurevar, ".mean", sep="") ] <- measurevar names(datac)[ names(datac) == paste(measurevar, ".sd", sep="") ] <- &...

Hi all, I am trying to fit t distribution using the function "tFit" in the library(fBasics). I am using the code tFit(datac[[2]]) and it returns the following list. Title: Student-t Parameter Estimation Call: tFit(x = datac[[2]]) Model: Student-t Distribution Estimated Parameter(s): df 78.4428 I just wonder how can I refer to the estimated parameters. I tried tFit(datac[[2]]) $df,tFit(datac[[2]])@df, but...

...2,9,2,4,3,8,5,10,5,4,4,7,4,7,7,4,16,4,3,12,4,3,10,27,4,14,7,10,4,9,8,3,5,3,8,5,2,8,38,3,3,3,4); plot(x=c(1,2),y=rep(1,2),ylim=range(data),xlim=c(0.5,2.5),type="n",main=c("Adjusted boxplot vs. Classical boxplot"),ylab=c("Loss Data"),xlab=c(""),axes=F) adjbox(x=datac,pars=list(outpch=16,outcex=2,outcol="light blue"),add=T,at=1,col="light blue") boxplot(x=datac,pars=list(outpch=16,outcex=2,outcol="light green"),add=T,at=2,col="light green") axis(2,at=seq(0,max(data)+10,by=10),col="orange",labels=F) axis(1,at=c(1,...

Dear All- I am trying to merge two data files - they have different date formats and different times zones. I need to match up the date/time of the datasets and then invoke a conditional statement, such as: if dataC$mph is >= 12 then keep dataM$co23 for the corresponding time/date stamp. snippets of data files: *dataC.txt* LST in mph Deg DegF DegF2 % volts Deg mph2 w/m2 0506010000 0.00 13.6 218.1 36.8 -999 65.1 -999 -999 18.2 0.2 0506010005 0.00 12.9 21...

Celliax development and download site: www.celliax.org Celliax is a channel driver for the Asterisk Free PBX that manages GSM and CDMA cellular phones through an adapter, composed by a datacable (for commands) and an audiocable (for the voice) interfacing the computer soundcard. chan_celliax is also capable of making and receiving Skype calls through the cellphone, and has an app like app_directory that let you choose which one of your Skype contacts you want to call. Celliax runs on...

Hi, I have two data sets look like below: ========================== state count1 percent1 CA 19 0.34 TX 22 0.35 FL 11 0.24 OR 34 0.42 GA 52 0.62 MN 12 0.17 NC 19 0.34 state count2 percent2 FL 22 0.35 MN 22 0.35 CA 11 0.24 TX 52 0.62 ========================== How to combine these two data set and make it look like below?

...eferences and has the same citation as a pers.comm. The problem for my paper can be resolved by removing the citation from references to text, but this seems unfair to R and the package authors for AGU publications. AGU does recognize some permanent data archives listed here http://www.agu.org/pubs/datacent.html - maybe it would make sense to have a mirror at one of these permanent archives for packages and versions. Regards, Steve McIntyre ----- Original Message ----- From: "Achim Zeileis" <Achim.Zeileis at wu-wien.ac.at> To: "Steve McIntyre" <stephen.mcintyre at u...

I'd like to set up a Samba DC, an additional DC and a Samba file server. Everything is configured with RAID 1. I encrypt all partitions except partition /boot (which is also on RAID). What partitions, partition sizes, file system types do you use on your Samba DC, additional DC and file server? Is there something like "best practice"? kind regards, birgit

...your > Samba DC, additional DC and file server? Is there something like "best > practice"? I usually have the following default partitions (LV) on my production servers / 5 GB /usr 5 GB /var 5 GB /tmp (1-n GB) /home (1-n GB) swap (4 GB) + every shared folder or area datac can grow is on a separate LV, too All filesystems are XFS. Hardware raid level 6. You should keep in mind, that you may have folders/areas that can be filled by services or users with data. And if the volume they are on run out of space, it could happen, that other services stop running. An examp...

...s there something like "best >> practice"? > > >I usually have the following default partitions (LV) on my production >servers > >/ 5 GB >/usr 5 GB >/var 5 GB >/tmp (1-n GB) >/home (1-n GB) >swap (4 GB) >+ every shared folder or area datac can grow is on a separate LV, too >All filesystems are XFS. Hardware raid level 6. > >You should keep in mind, that you may have folders/areas that can be >filled by services or users with data. And if the volume they are on run >out of space, it could happen, that other services sto...

Hi, I am trying to fit t copula to some data, and I am using the following function in the library(QRMlib). Udatac <- apply(datac, 2, edf,adjust=1) tcopulac <- fit.tcopula.rank(Udatac) But the error message come out "Error in fit.tcopula.rank(Udatac) : Non p.s.d. covariance matrix" Could anyone give me some advice? In fact, I am not sure what the "adjust=1" is used for. Many thanks....

Hello! If I have a matrix as 1 2 2 3 and I want to change the value 2 in 0, what can I do? Thank you -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)