search for: 0.3333

Hi, Does anybody know, what is going on here? > diag(sqrt(1)) [,1] [1,] 1 > diag(sqrt(0.3333)) <0 x 0 matrix> > sqrt(1) [1] 1 > sqrt(0.3333) [1] 0.5773214 BR, Markku Karhunen researcher University of Helsinki

Hello, My name is Rodrigo, I am using R program and I have a trouble. I am trying to do a dendrogram with genetics information. Let me explain... The Similarity Matrix was already did, and with this matrix I want to construct a dendrogram. So, the distance is done. I need to transform this matrix (that I have) in a dendrogram, I woud be very grateful if someone could help me. PS: I am sending

Hi all, I noticed something strange when I ran aov and anova. vtot=c(7.29917, 7.29917, 7.29917) #identical values fac=as.factor(c(1,1,2)) #group 1 has first two elements, group 2 has the 3rd element When I run: > anova(lm(vtot~fac)) Analysis of Variance Table Response: vtot Df Sum Sq Mean Sq F value Pr(>F) fac 1 1.6818e-30 1.6818e-30 0.3333 0.6667 Residuals 1

Hi, I am using ccf but I could not figure out how to calculate the actual lag in number of periods from the returned results. The documentation for ccf says:"The lag is returned and plotted in units of time". What does "units of time" mean? For example: > x=ldeaths > x1=lag(ldeaths,1) > results=ccf(x,x1) > results Autocorrelations of series 'X', by lag

Hi all, I have a package that is specific to a task I was repetitively using a few years ago. I now needed to run it again with new data. However I am told it was built with an older version or R and will not work. How can I tweak the package so it will run on 11.1? It was a one-off product and has not been maintained. Is there a way to "unpackage" it and repackage it to work? I

signif() is primitive in 2.6.x, and so uses positional matching (almost all primitives do). See the NEWS file for the announcement of the change. On Sat, 2 Feb 2008, ben.hansen at umich.edu wrote: > Full_Name: Ben Hansen > Version: 2.6.1 > OS: Windows > Submission from: (NULL) (66.93.3.101) > > > Hi, > > In R 2.6.0 or 2.6.1 on Windows, I get the following upon

Full_Name: Marek Ancukiewicz Version: 2.10.1 OS: Linux Submission from: (NULL) (74.0.49.2) Both cor() and cor.test() incorrectly handle ordered variables with method="kendall", cor() incorrectly handles ordered variables for method="spearman" (method="person" always works correctly, while method="spearman" works for cor.test, but not for cor()). In

I have data as below.Please let me know how the ACF and Pacf used to determine the order od arima model. Is there any rules need to be followed to determine order.Please advise > turkey.price.ts Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2001 1.58 1.75 1.63 1.45 1.56 2.07 1.81 1.74 1.54 1.45 0.57 1.15 2002 1.50 1.66 1.34 1.67 1.81 1.60 1.70 1.87 1.47 1.59 0.74 0.82

Dear all I am having a linear system of the form A*X=B and I want to find the X by using least squares. For example my A is of dimension [205,3] and my B is of dimension[205,1]  I am looking for the X matrix which is of the size [3,1]. In the matlab I was doing that by the function  X = LSCOV(A,B) returns the ordinary least squares solution to the     linear system of equations A*X = B, i.e., X

Full_Name: t. avery Version: 2.0.0 OS: windows xp / Linux Submission from: (NULL) (131.162.134.159) ks.test does not produce the correct output. If given the script: d1 <- c(53.63984674,0.383141762,1.915708812,0.383141762,10.72796935,6.896551724,20.30651341,5.747126437,0) d1 d2 <- c(76.43312102,15.2866242,3.821656051,1.27388535,0,0.636942675,1.27388535,0.636942675,0.636942675) d2

I ran a negative binomial logit hurdle model and am now trying to plot the effects of a continuous predictor variable (the only variable in my model) on the count and zero component and the overall mean response. I'm confused because for some values, the predicted overall mean is higher than the mean of the non-zero counts (range of predicted overall means=2.2-11.0; range of non-zero count

Full_Name: Ben Hansen Version: 2.6.1 OS: Windows Submission from: (NULL) (66.93.3.101) Hi, In R 2.6.0 or 2.6.1 on Windows, I get the following upon opening the GUI (no previous commands or special settings): > signif(digits=4, x=1/3) [1] 4 It seems to be taking 4 to be the "x" argument, the number to be rounded. However, my understanding (perhaps mistaken) was that it should

Hello R community, I have a question about the logistic regression function. Specifically, when the predictor variable has not just 0's and 1's, but also fractional values (between zero and one). I get a warning when I use the "glm(formula = ... , family = binomial(link = "logit"))" which says: "In eval(expr, envir, enclos) : non-integer #successes in a binomial

I have created a tree and want to save some of the data so that I can create a html table from it. I would like to save the output from data.ltr (see example below) to a file, but haven't found a way to do that, keeping the nice format that typing data.ltr gives me (see output below). Is there a way to do this? Example: library (maptree) library (tree)

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100625/44fbea79/attachment.pl>

Hi list, I have the following code to compute the acf of a time series acfresid <- acf(residfit), where residfit is the series when I type acfresid at the prompt the follwoing is displayed Autocorrelations of series ?residfit?, by lag 0.0000 0.0833 0.1667 0.2500 0.3333 0.4167 0.5000 0.5833 0.6667 0.7500 0.8333 1.000 -0.015 0.010 0.099 0.048 -0.014 -0.039 -0.019 0.040 0.018

I am involved in a study where, as in most of life, men demonstrate themselves to be recalcitrant. So while we have many probands and most of their mothers we only have about 50% of the trios being complete. I have been running tdt and trio.types. It appears as if it is ignoring the duos. Sometimes a duo can be informative. For instance Father ..missing Mother 1/2 Proband 1/1 This duo shows that

Hello everyone, I am using *mlogit* to analyse my choice experiment data. I have *3 alternatives* for each individual and for each individual I have *9 questions*. I have a response from *516 individuals*. So it is a panel of 9*516 observations. I have arranged the data in long format (it contains 100 columns indicating different variables and identifiers). In mlogit I tried the following

This was filed as a bug report on the Debian r-base package. It is more properly a bug report on the ctest package in R. The default method for wilcox.test manipulates x and y without checking the class or data.class of these objects. Possible solutions are - create wilcox.test.factor (if appropriate) - check the class and/or data.class of x and y in wilcox.test.default and produce error

To whom it may concern, I am a student at Macaleste College, and next semester Macalester is going to offer a course for CellBio that is mainly statistically based. For the most part the students will be using R for analysis. The problem is there will be some simple differential equations for the students to solve. The committee that in charge of the classes corriculam would like only to