similar to: Minor bug in lillie.test from nortest package (PR#9784)

Hi, I'm a newbie and i am unable to use lillie.test in nortest I have a message: "Couldn't find function "lillie.test" I am under windows2000 with R1.8.1 nortest is listed with .packages(TRUE) How to do to use lillie.test function? Laurent Houdusse Analyste Programmeur

Hi all, I have a dataset of 2000 numbers ( it's noise measured with a scoop ) Now i want to know of my data is normal distributed (Gaussian distribution). I did already: - 68-95-99.7 test - Q-Q-plot and now i used "nortest library" and the Lilli.test() However i don't understad the output? lillie.test(z) Lilliefors (Kolmogorov-Smirnov) normality test data: z D =

I believe there to be a bug in the cvm.test module of the nortest package authored by Juergen Gross. I do not know how to contact the author directly. I've been running some normality tests using the nortest package. For some of my datasets the Cramer-von Mises normality test generates an extremely high probability (e.g., 1.637e+31) and indicates normality when the other tests do

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1,

I am reading regressors from an excel file (I have no control over the file) and some of the element names have spaces: i.e. "Small Bank Aquired" but I have found that lm(SourceData ~ . - "Small Bank Aquired", mcReg) doesn't work (mcReg = modelCurrentRegressors) As they are toggles I have ran them through factor() to be treated propertly as 0 or 1 but due to the fact I

Hi, this is my first help request so please bear with me. I've been running some normality tests using the nortest package. For some of my datasets the Cramer-von Mises normality test generates an extremely high probability (e.g., 1.637e+31) and indicates normality when the other tests do not. Is there something I'm misunderstanding or potentially a bug in the code? Below are the

Good Night I made different test to check normality and multinormality in my dataset, but I don´t know which test is better. To verify univariate normality I checked: shapiro.test, cvm.test, ad.test, lillie.test, sf.test or jaque.bera.test and To verify multivariate normal distribution I use mardia, mvShapiro.Test, mvsf, mshapiro.test, mvnorm.e. I have a dataset with almost 1000 data and 9

Hi experts of R, Are there any functions in R to test a univariate series for long memory effects, structural breaks and time reversability? I've found for ARCH effects(ArchTest), for normal (Shapiro.test, KS.test(comparing with randn) and lillie.test) but not for the above mentioned. Where can I find a comprehensive list of functions available by type? Thank you Renato Costa -- View this

I don't know from the nortest package, but it should ***always*** be the case that you test hypotheses H_0: The data have a normal distribution. vs. H_a: The data do not have a normal distribution. So if you get a p-value < 0.05 you can say that ***there is evidence*** (at the 0.05 significance level) that the data are not from a normal distribution. If the nortest package does

Hi, I use ks.test or lillie.test to verify a normal distribution. It's performed for a group My users use SigmaStat software and a One Way ANOVA on several groups In the result page there is a probability value to determine if Normality test is failed or passed So, how can i retrieve this probability value on several groups? Is there another function in R to verify normality on several

Hello all, Would somebody be kind enough to show me how to do a KS test in R for a lognormal distribution with ESTIMATED parameters. The R function ks.test()says "the parameters specified must be prespecified and not estimated from the data" Is there a way to correct this when one uses estimated data? Regards, Kwabena. -------------------------------------------- Kwabena Adusei-Poku

Just to inform: I posted that before in R-sig-ecology but as it might be interesting also for other useRs, I post it also to the general r-user list: Hello Alexandre, thank you very much. I also found another way to extract summarizing information from lm results over e.g. 1000 repeated model runs: results2 <- t(as.data.frame(results)) summary(results2) Although some questions popped up in

Dear List, I have a problem I'm finding it difficult to make headway with. Say I have 6 ordered observations, and I want to find all combinations of splitting these 6 ordered observations in g groups, where g = 1, ..., 6. Groups can only be formed by adjacent observations, so observations 1 and 4 can't be in a group on their own, only if 1,2,3&4 are all in the group. For example,

There are probably many mistakes with this code. I am used to coding in C with a debugger, so I am very new to coding in R without one and different syntax. My code in the upper left panel of R Studio isz<-c(3,1,4,5,2,1,0,3,5,8)zlength(z)y<-c(0,0,0,0,0,0,0,0,0,0)ylength(y)zdiff = diff(z)zdiffn<-length(zdiff)nx<-zdifff<-function(x) { for (k in 1:n){

Is there a way to have users of the samba server, but not add them by smbpasswd -a <UserID>? I want the samba server to be a domain member and the users to only authentic from the PDC. These are the steps that I have attempted: Users are in both the Windows domain and the UNIX NIS account 1. smbpasswd -j <Domain> -r PDC -U <admin> Joined the Domain 2. edited the

Hi, I am new to R and need to read in a file with 19 columns and 7000 rows and make it into a list of 7000 lists with 19 items each. For a simpler case of 10 by 10 table I used x <-scan("file", list(0,0,0,0,0,0,0,0,0,0)), perhaps clumsy, but it did the job. However with the large 19x7000 (which needs to be transposed) I am not sure how to go about it. Coudl somebody suggest a way?

Hi R community I have a little problem, and I tried to solve it by myself but I couldn't. I building an if loop, and I want to check a value inside an interval. This would be the case: pvalue=0,2999 if(pvalue>0.05 or pvalue<0.1) as you can see I would like to check in that if loop if my pvalue is inside of that interval(from 0.05 to 0.1), and I tried these options: if(pvalue>0.05

Hi there, I want to do a nonparametric covariance analysis and I have tried to use the package "sm" function "sm.ancova" but it didn't work for me because I have more then one covariates (I have 18 covariates and 3 factors). I want to analyse for one factor (who has 13 levels) where the differences for my response are, using the explanatory (covariates). (The other two

Hi all, If you take a look at the service.rb file in the win32-service repository (the new one in the toplevel repository path), I''ve got this bit of code, which succeeds, but I can''t seem to unpack the data structure properly. Did I pack it wrong to begin with? I should know this but I''m spacing out. proc_status =

Dear R-devel, My current research focuses upon sensitivity analyses which require [1] clusterings of patients in a baseline covariate X-space and [2] examining the distribution of within-cluster treatment differences in outcome. I have generated some primitive R code for this, but I really need to be able to interact with graphical displays of within-cluster information (local average