search for: norm

...t I am not sure. Version 1 ##t test len_sig_test = length(sig_test) if(len_sig_test%%2 > 0) message("Number of groups for t-test is not paired.") if(max(sig_test) > group) message("Group number for t-test is out of range.") data.norm.test = data.norm.zo #normalized data nrow <- nrow(data.norm.test) ncol <- ncol(data.norm.test) p.ttest = rep(NA, nrow) #parametric test -- t test p.wilcox = rep(NA, nrow) #non-parametric test --...

...is no insect count at all. So, I tried fitting 'zero inflated Poisson model' to the data using cozigam function as follows: > chf.res1 <- cozigam(chafer ~ s(longitude, latitude) + s(PC1) + s(PC2) + s(PC3) + s(PC4) + (year) + paddock2, family=poisson, data=df.chf09) iteration = 2 norm = 1.023552 iteration = 3 norm = 0.4258558 iteration = 4 norm = 0.003208923 iteration = 5 norm = 0.006238185 iteration = 6 norm = 0.006029638 iteration = 7 norm = 0.005037579 iteration = 8 norm = 0.004390153 iteration = 9 norm = 0.003970208 iteration = 10 norm = 0.003670730 it...

This issue ties loosely into other recent S4 topics on this board. The methods package defines a number of implicit generics for linear algebra related functions (rcond, norm, backsolve) that, when used, interfere with base package operations. Here is the cut-and-paste version of the code the illustrates the problem: # rcond x1 <- cbind(1, 1:10) rcond(x1) setGeneric("rcond") rcond(x1) # norm example(norm) setGeneric("norm") example(norm) # ba...

dear list members, there seems to be a problem with the prelim.norm function (package norm) as number of items in the dataset increases. the output of prelim.norm() is a list with different summary statistics, one of them is the missingness indicator matrix "r". it lists all patterns of missing data and a count of how often each pattern occured in the da...

Dear alltogether, I experience very strange behavior of imputation of NA's with the NORM library. I use R 2.2.0, win32. The code is below and the same dataset was also tried with MICE and aregImpute() from HMISC _without_ any problem. The problem is as follows: (1) using the whole dataset results in very strange imputations - values far beyond the maximum of the respective column,...

Hi, I was trying to read in a file test.txt, which has the following data: norm norm norm class class class a 1 2 3 4 5 6 b 3 4 5 6 7 8 c 5 6 7 8 9 10 in my R code, I do the following: --------------------------------------------------- > mat <- read.table('test.txt',header=T,row.names=1,sep='\t') > mat norm norm.1 norm.2 class class.1 class.2 a...

Dear All, Trying to install package 'norm' under linux-mandrake8.2 I've got the following messages: R CMD INSTALL norm_1.0-9.tar.gz * Installing *source* package 'norm' ... ** libs g77 -mieee-fp -O3 -fomit-frame-pointer -pipe -mcpu=pentiumpro -march=i586 -fno-fast-math -fno-strength-reduce -O3 -fomit-frame-pointer -pipe...

...brary(outliers) filename<-"timediff.csv" time<-read.csv(filename, sep=";",header=TRUE,dec=".") # Esto es lo que has de cambiar chisq.out.test(time$TimeDiff) # Ayuda adicional #--------------------------------------------------------------------- # Si no supones normalidad #--------------------------------------------------------------------- time.without.outs <- data.frame(row.names= 1:dim(time)[1], dif = time ) x <- boxplot(time.without.outs$TimeDiff,range= 1.5) # Quita los outliers que se van más de 1.5 rangos itercuartilicios. Puedes poner otro crit...

Dear All, I have the following data coming out from s <- with(final, summarize(norm, llist(gtt,fdiab), function(norm) { n <- sum(!is.na(norm)) s <- sum(norm, na.rm=T) binconf(s, n) }, type='matrix') ) ie gtt fdiab norm.norm norm.norm2 norm.norm3 18...

Hello all, I am trying write a fairly simple function that provide a quick way to calculate several distributions for a data set. I am trying to provide a function that has a argument that specifies which distribution is outputted (here "norm" or "cumu"). I also have a melt argument but that seems to be working fine. I have been able to get my function working well for just one distribution but when I add another and try to add a "dist.type" argument (with potential values "cumu" and "norm"),...

Greetings! I am bootstrapping and I am using EM in the norm package to fill in missing data for a financial time series with each step of the loop. For the most part EM works fine for me, but the following error message is guaranteed before I hit the 200th scenario: Iterations of EM: 1...2...3........348...349...Error: NA/NaN/Inf in foreign function call (...

Hello, I apologize if an answer to my questions is available, or if I submitted this question incorrectly. I have read the mailing lists, as well as the R Project and CRAN homepages. However, I may have missed something. I noticed the package 'norm' has been removed. Its page http://cran.r-project.org/web/packages/norm/index.html now reads: "Package ?norm? was removed from the CRAN repository. Formerly available versions can be obtained from the archive." My questions are: 1. Why was norm removed? I used the package and foun...

...n't found anything relevant to my problem. I guess it is something simple...I just cannot see it. I have the following data frame, named "data", which is a subset of a microarray experiment. The columns are the samples and the rows are the probes. I binded the first line, called "norm", which represents the estimated output. I want to create a linear model which shows the relationship between the gene expressions (rows) and the output (norm). *data* GSM276723.CEL GSM276724.CEL GSM276725.CEL GSM276726.CEL norm 0.897000 0.590000 0.683000...

Hi Folks, I'm using the 'norm' package (based on Shafer's NORM) on some data. In outline, (X,Y) are bivariate normal, var(X)=0.29, var(Y)=24.4, cov(X,Y)=-0.277, there are some 900 cases, and some 170 values of Y have been set "missing" (NA). The puzzle is that, repeating the multiple imputation starting from...

...little green book on missing data (Sage 2002). The dataset for which I am having problems, "usnews", can be found at: http://www.ats.ucla.edu/stat/books/md/default.htm. I am working on a Windows machine with R 2.5 installed, all packages up-to-date. The problem has to do with the prelim.norm() function of the package "norm". Specifically, I need to use this pre-processing function to later use the EM algorithm and DA procedures in the norm package. I am getting an error with the following code. ---------------------- > pre <- prelim.norm(usnews) Error in as.double.d...

Hey eveyone! I hope someone can help wiht this question. I have a matirux of all zeros and ones and I would like to indentify all unique patterns in the rows andthe number of times the pattern occurs. I changed all zeros to NA tried to use prelim.norm to identify all patterns of missing data in the rows. I got the message Warning message: NAs introduced by coercion Any ideas of how to get this to work? Or are there any way to indentify all the unique patterns in a huge matrix? ( 10000 x 71) Thanks for any suggestions!! E...

...;-read.csv(filename, sep=";",header=TRUE,dec=".") # Esto es lo que >> has de cambiar >> chisq.out.test(time$TimeDiff) >> >> >> # Ayuda adicional >> #--------------------------------------------------------------------- >> # Si no supones normalidad >> #--------------------------------------------------------------------- >> time.without.outs <- data.frame(row.names= 1:dim(time)[1], dif = time ) >> >> >> x <- boxplot(time.without.outs$TimeDiff,range= 1.5) # Quita los outliers >> que se van más de...

Given the following data, I am plotting log.det ~ norm.beta, where the points depend on a parameter, lambda (but there is no functional form). I want to find the (x,y) positions along this curve corresponding to two special values of lambda lambda.HKB <- 0.004275357 lambda.LW <- 0.03229531 and draw reference lines at ~ -45 degrees (or normal...

On Sun, 10 Oct 2004, RenE J.V. Bertin wrote: > Could I ask some hints/help in interpreting the following ANOVA results, > please? This concerns an experiment where I study the incidence and > severity of motion sickness. I have Sickness.norm, a subjective > discomfort/sickness estimate, normalised to 0..1, the session time T > (normalised to 0..1 and binned in 0.2 wide bins) and a qualitive > indicator if a given Subject was sick or not. For instance, to see if > there is an effect of time, I do > > > > summar...

Hello all I'm running da.norm function in R for climate data rngseed(1234567) theta1=da.norm(mydata, thetahat, steps=1000,showits=T) param1=getparam.norm(mydata,theta1) As I understand the 1000 steps represent the markov chain values. Is there a way to plot them? Something like plot(1:1000, param1$mu[]). I just can't fin...