Re: [R] how to keep very small or large number?

- plogis(z, log.p=TRUE, lower.tail=as.logical(y)) ln.p. - sum(ln.p) log10.p. - ln.p./log(10) printSmall - function(ln.x){ l10x - ln.x/log(10) paste(10^(l10x%%1), e, l10x%/%1, sep=) } printSmall(ln.p.) 6.43113656279925e-5355 How's this? spencer graves Cunningham Kerry wrote

Re: [R] FDR analyses: minimum number of features

this, you could estimate probabilities for obtaining results that were bogus vs. marginal, barely useful vs. highly accurate and plot them vs. alternative budgets, etc. I hope this comment makes more sense than my earlier nonsense. spencer graves Dupont, William wrote: I

Re: [R] FDR analyses: minimum number of features

this. spencer graves Dupont, William wrote: Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how

Re: [R] factor as seq() in for loop

Does the follow help? aFactor - factor(rep(letters[1:3], 1:3)) aFactor [1] a b b c c c Levels: a b c for(af in levels(aFactor)){ + print(sum(aFactor==af)) + } [1] 1 [1] 2 [1] 3 spencer graves Tord Snall wrote: Dear all, I would like to use the values

Re: [R] Are least-squares means useful or appropriate?

, spencer graves Felipe wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi. I have been reading about the convenience of using least-squares means (a. k. a. adjusted means) in multiple comparisons (I used to resort to them when using SAS). I even read a post in this list warning

Re: [R] logistic regression with nominal predictors

than one could get in a typical case using Monte Carlo. And Monte Carlo in R is normally fairly easy even for quite complicated situations. spencer graves Ramón Casero Cañas wrote: (Sorry for obvious mistakes, as I am quite a newby with no Statistics background). My question

Re: [R] nls()

that effectively. It's not trivial, but it can be made to work. spencer graves J.M. Breiwick wrote: Hi, I am using nls() with the form: nls(~my.fcn(...)) because I have to iteratively compute the expected y values. The function my.fcn() returns y.obs-y.pred However, I want to fix

Re: [R] Random effect model

in a response to a previous question. Good Luck. spencer graves Hathaikan Chootrakool wrote: Dear R-help group, I would like to model directly following random effect model: Y_ik = M_ik + E_ik where M_ik ~ N(Mew_k,tau_k^2) E_ik ~ N(0,s_ik^2

Re: [R] How to sort data sets?

__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite

Re: [R] help for linear-circular correlation

increase the speed and utility of replies; http://www.R-project.org/posting-guide.html;.) Good Luck, Spencer Graves (a) http://www.pstat.ucsb.edu/faculty/jammalam/html/research%20publication_files/article082003.pdf ecoinfo wrote: Hi R-profs, Maybe my question is a little off

Re: [R] regression with restrictions - optimization problem

. Good luck. spencer graves Mark Hempelmann wrote: Dear WizaRds! I am sorry to ask for some help, but I have come to a complete stop in my efforts. I hope, though, that some of you might find the problem quite interesting to look at. I have been trying to estimate

Re: [R] how to do multiple comparisons in R?

! http://www.R-project.org/posting-guide.html). spencer graves Hongyu Sun wrote: Hi, Sorry I have to bother a question. Does R have the functions to do lsd, tukey, bonferonni, contrast etc. like in SAS? Many thanks, HS __ R-help

Re: [R] Time series ARIMAX and multivariate models

capabililities for time series. A more specific question might produce more useful replies. The posting guide (http://www.R-project.org/posting-guide.html) might help increase the chances that someone else will provide the help you seek. Good Luck, Spencer Graves [EMAIL

Re: [R] Multinomial Logit and p-values

response quicker. I believe that many people have found that the posting guide! (http://www.R-project.org/posting-guide.html) was quite useful for helping them prepare a question that produced more informative responses quicker. Best Wishes, spencer graves Sangick Jeon

Re: [R] test for exponential,lognormal and gammadistribution

of the remaining 72 may help you more. spencer graves [EMAIL PROTECTED] wrote: hello! i don't want to test my sample data for normality, but exponential- lognormal- or gammadistribution. as i've learnt the anderson-darling-test in R is only for normality and i am not supposed

Re: [R] GlmmPQL help

clearly the degrees of freedom in question, what you think the number should be and why. I believe that will increase the chances that you will get a useful reply quickly. I'm sorry I couldn't help more. spencer graves Mark Steer wrote: Hi, I'm running a GLMM on binomial

Re: [R] Simulate phi-coefficient

: VECTOR1-rep(c(1,0),c(15,10)) n - length(VECTOR1) library(MASS) Sig - array(c(1,.5, .5, 1), dim=c(2,2)) z - mvrnorm(n, mu=c(0,0), Sigma=Sig) p - mean(VECTOR1) q. - qnorm(p) Z - 0+(zq.) hope this helps. spencer graves Bliese, Paul D LTC USAMH wrote: Looking for help

Re: [R] Interpolating / smoothing missing time series data

interested in learning how to reproduce in R all the examples in that book, and I'd happily receive suggestions from others on how to do that. spencer graves Not sure if this answers your question but if you are looking for something simple then na.approx in the zoo package will linearly

Re: [R] Time Series Analysis: book?

feel I need to get started, though I found the vignettes with zoo to be quite helpful. spencer graves Jean-Luc Fontaine wrote: Wensui Liu wrote: TS is a huge topic. The book recomended by statisitcian might be different from the one recommended by econometrician. Finance guy might

Re: [R] Doubt about nested aov output

2 1.818989e-12 1 I don't know if this relates to Ronaldo Reis' question, but it certainly seems plausible. spencer graves Christoph Buser wrote: Hi I think the problem is related to specifying Tratamento both as a fixed factor and in the error term. I

Re: [R] partial association model

that on average, posts prepared following this guide get more useful answers quicker. spencer graves [EMAIL PROTECTED] wrote: my last post was filtered,so I post it again with another title. If I do not make a mistake,the partial association model is an extension of log-linear model.I

Re: [R] question sur R

was not of more help, but I don't have time to trace this problem myself. spencer graves Abdelhafid BERRICHI wrote: hello I've tried to simulate a normal law, like that : X1 = c(rnorm(90,50,5583),rnorm(160,1198,13034597),rnorm(40,13,125)) then, I've regressed my ordinal polytomic

Re: [R] Doubt about nested aov output

in R 2.1.1, could you please provide this list with a small, self contained example that displays the symptoms that concern you? And PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html;. It might increase the speed and utility of replies. spencer graves

Re: [R] Distributional characteristics

.) spencer graves Naji wrote: Hi all I've a continuous variable and I want to test (graphically, plotting observed and theoretical distr, qqplot) whether it follows some formal distribution. (some thing close to Ricci document : Fitting distributions with R, Feb05). The distribution I want

Re: [R] R package for ICA

I just got 48 hits from RSiteSearch(independent component analysis), the first of which was http://finzi.psych.upenn.edu/R/library/mlica/html/mlica.html;. hope this helps. spencer graves p.s. I believe people who follow the posting guide typically get more useful

Re: [R] Question

recommend this book. spencer graves Mahmoud Torabi wrote: Dear Sir/Madam I would be pleased if anybody can help me. I'm using linear mixed model (lme) function.I'm doing some simulation in my research and need to be assigned variance components values during of my program

Re: [R] lme and ordering of terms

the choice of coordinate axes but does NOT change the subspaces involved UNLESS you extract the terms in a different order. For more detail, consult any good book on regression. spencer graves Christoph Scherber wrote: Dear R users, When fitting a lme() object (from the nlme library

Re: [R] loglinear model selection

) spencer graves Bernardo Rangel Tura wrote: Hi R-masters! I have a problem and need your help. I have 9 discrete variables with 2 levels each. In exploratory analisys I generate one matrix with chi-square for tables with 2 ariables each with this script setwd(F:/) dados-read.csv

Re: [R] under sample problem

Does the following help: set.seed(1) pop - 1:100 s1 - sample(pop, 5) s1. - pop[!(pop %in% s1)] s2 - sample(s1., 5) s1 [1] 27 37 57 89 20 s2 [1] 91 94 66 62 6 spencer graves [EMAIL PROTECTED] wrote: Hello, I have a problem to treat my data. I seek

Re: [R] Superassignment (-) and indexing

f1 to call f2 and have f2 change x in f1. However, if your f2 gets called some other way or if the name of x is misspelled or changed in either f1 or f2, we could be back to the situation I just described. spencer graves Brahm, David wrote: In a clean environment under R-2.1.0

Re: [R] Superassignment (-) and indexing

Wow! That's great. Thanks. spencer Thomas Lumley wrote: On Fri, 2 Sep 2005, Spencer Graves wrote: Permit a mild protest on the word appropriate in this context. The global assignment operator - provides, for my tastes, excessive opportunities for problems. If I

Re: [R] Pseudo-Voigt fit

they do, it's rarely helpful. I've reluctantly learned that there is often no substutute for reading the [EMAIL PROTECTED] manual.) I'd be shocked if this answered your question, but I hope it is helpful nonetheless. spencer graves [EMAIL PROTECTED] wrote: Hi, I am

Re: [R] Defining an ex-gaussian PDF

reply quickly. spencer graves Michael Kubovy wrote: How does one define PDFs as yet undefined in R, such as the ex- gaussian, the sum of two RVs, one exponential, one Gaussian? The PDF would then be the convolution of an exponential PDF, dexp(), and a normal, dnorm(). Kindly

Re: [R] basic anova and t-test question

the ensemble is not. spencer graves [EMAIL PROTECTED] wrote: Hello, I'm posting this to receive some comments/hints about a rather statistical than R-technical question ... . In an anova of a lme factor SSPos11 shows up non-significant, but in the t-test of the summay 2 of the 4

Re: [R] question about custom contrasts in ANOVA

to worry about which answer seems most correct, and I also know something about the limits of the conclusions. I know this doesn't answer your question, but I hope it helped with a solution methodology. Best Wishes, Spencer Graves Scot W McNary wrote: Hi, I

Re: [R] question about custom contrasts in ANOVA

Dear Prof. Ripley: Agreed. I thought that coding the factor levels manually and / or looking at model.matrix(y~xma) after further simplifying the problem might help Scot clarify his question. Thanks, Spencer Graves Prof Brian Ripley wrote: I was puzzled

Re: [R] Panel regression in R

, but it seems to me to be quite helpful when followed. spencer graves Haibo Huang wrote: I am currently trying to replicate the results I got from RATS for a panel regression. The codes in RATS looks like this: * Final equation for Office Cap Rate Spread * Regression, Panel Data

Re: [R] PDL model

From 'RSiteSearch(polynomial distributed lag)', I learned that there was a discussion of this issue not quite two years ago, which may interest you. In brief, Thomas Lumley provded code with some noted limitations. spencer graves [EMAIL PROTECTED] wrote: Dear r-help

Re: [R] question sur R

anglais q'en francais. (2) PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Bon chance. Spencer Graves Abdelhafid BERRICHI wrote: bonjour je suis elève à l'ENSAI Rennes et je suis actuellement en stage de fin d'études j'ai une

Re: [R] How to add values on the axes of the 3D bi-variable lrm fit?

that and it worked. 6. In the future, I believe you can increase your chances of getting a useful reply quickly if you PLEASE do read and follow the posting guide! http://www.R-project.org/posting-guide.html;. Best Wishes, spencer graves Jan Verbesselt wrote: Dear

Re: [R] GLM/GAM and unobserved heterogeneity

. spencer graves Kyle G. Lundstedt wrote: Hello, I'm interested in correcting for and measuring unobserved heterogeneity (missing variables) using R. In particular, I'm searching for a simple way to measure the amount of unobserved heterogeneity remaining in a series

Re: [R] How to do log normal regression?

get quicker and more useful replies to your questions if you please read and try to follow the posting guide! http://www.R-project.org/posting-guide.html;. Good luck. spencer graves Haibo Huang wrote: I want to fit a Log-Normal CDF function between two variables

Re: [R] Copying rows from a matrix using a vector of indices

! http://www.R-project.org/posting-guide.html;. It can increase the chances that you will get a useful reply. spencer graves Martin Lam wrote: Hi, I am trying to use a vector of indices to select some rows from a matrix. But before I can do that I somehow need to convert

Re: [R] how to reshape an array avoiding for loops

Have you considered aperm? I found this listed under See Also for ?t. spencer graves Wladimir Eremeev wrote: Dear r-help, I have an array a1 with dimensions [1:660,1:65,1:25] I would like the first dimension to be the last one. That is I want and array [1:65,1

Re: [R] specify seed for Random Number Generator

Have you considered set.seed? This was listed with the first hit for 'RSiteSearch(random number seed)'. spencer graves Dhiren DSouza wrote: I need to generate 100 I.I.D samples from an exponential distribution. I use rexp(100,parameter). Is there anyway to specify

Re: [R] How to assess significance of random effect in lme4

by Googling for likelihood ratio with parameter at a boundary. My literature search on this issue is far from exhaustive, and I would be pleased to hear from others who have other articles to recommend. spencer graves # Andrew Robinson wrote: Hi Spencer, thanks

Re: [R] power of a matrix

Hi, Gabor: Thanks. spencer graves Gabor Grothendieck wrote: Its expm. On 8/17/05, Spencer Graves [EMAIL PROTECTED] wrote: Hi, Peter: I couldn't find mexp in the Matrix package, but I did find it in fMultivar and in Lindsey's rmutil. These are different functions

Re: [R] How to assess significance of random effect in lme4

Hi, Harold: Thanks for the clarification. I thought I had read the original post. Obviously, I had misread it. Thanks again. spencer graves Doran, Harold wrote: Yes, it is a different issue. ranef() extracts the empirical Bayes estimates, which are the empirical posterior modes

Re: [R] error in predict glm (new levels cause problems)

0.5775 0.5719 Residuals 17 9.1050 0.5356 anova(fit2) Analysis of Variance Table Response: y Df Sum Sq Mean Sq F value Pr(F) a 1 0.0008 0.0008 0.0015 0.9698 b 2 0.6186 0.3093 0.5775 0.5719 Residuals 17 9.1050 0.5356 spencer graves K. Steinmann wrote

Re: [R] How to assess significance of random effect in lme4

know. Comments? spencer graves p.s. It looks like [EMAIL PROTECTED] is a list containing vectors of length 29 and 6 in your example. I don't know what they are, but I don't see how they can be standard errors in the usual sense. Doran, Harold wrote: These are the posterior

Re: [R] power of a matrix

for mexp in fMultivar but not in rmutil. spencer graves Peter Dalgaard wrote: Rau, Roland [EMAIL PROTECTED] writes: Thank you very much! Thanks also to the authors of this function, Vincente Canto Cassola and Martin Maechler! This is exactly what I hoped for. look

Re: [R] Using nonlinear regression

. If you want examples of nonlinear regression, have you considered nls and optim? spencer graves Mark Miller wrote: Attached is a copy of my code, the data and the plots obtained by varying terms manually, I was told that nonlinear regression in R could find the values for me, but I

Re: [R] General expression of a unitary matrix

Could you provide an example that can NOT be expressed in that form? spencer graves J. Liu wrote: Thank you, Spencer. I read through the websites you suggested. What I need is how to parameterize a 2\times 2 unitary matrix. Generally, since for a complex 2\times 2 matrix

Re: [R] General expression of a unitary matrix

, etc., please post that to the list. spencer graves ## is.unitary - function(x, eps=prod(dim(x))*.Machine$double.eps){ x2 - (x %*% t(Conj(x))) (abs(mean(x2-diag(dim(x)[1]))) eps) } is.unitary(diag(3)+1e-15) is.unitary(diag(3)+1e-16

Re: [R] lme, fixed and random

lme(logratio ~ range + I(range^2), random=~1|fmla, data=data.frame.contianing.range.and.fmla) spencer graves p.s. I highly recommend Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer). Erin Hodgess wrote: Dear R People: I need to set up an lme. I have logratio

Re: [R] General expression of a unitary matrix

unitary matrix, which seems to require three real numbers, not four; with a fourth, you could get a general unitary matrix. spencer graves J. Liu wrote: Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries

Re: [R] help on cross hedge optimal hedge variance ratio

and respond to. On average, this tends to increase the speed, volume and utility of replies. spencer graves Krishna wrote: Hi everyone I am trying to estimate the optimal hedge variance ratio for cross hedging two commodities. the price levels are used (compared to price change

Re: [R] converting a t statistic to r2

for it. spencer graves -- Shaw, Philip (NIH/NIMH) wrote: HI I wonder if anyone can help. I have a longitudinal sample of 100 subjects: 200 data points were acquired starting at different ages and at irregular intervals (subjects have different numbers of repeated data points, so some

Re: [R] Problems with numeric variable containing ? values

to understand what you are asking and provide a useful reply. spencer graves [EMAIL PROTECTED] wrote: Hello, I have a problem with the values for one column in a table. The variable represented in this column is numeric (I get TRUE when I ask is.numeric(x)). However, the values are listed

Re: [R] question

different from this, PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html;. It may help you formulate your question to increase the chances of getting a more useful reply. Best Wishes, spencer graves WFJ van IJcken wrote: Hi, I have a problem

Re: [R] p-values

the posting guide. People report having found answers to difficult questions in the process of preparing a question following that guide, and when they do post a question, they are much more likely to get a useful reply. spencer graves Peter Ho wrote: Spencer, Thank you

Re: [R] get the wald chi square in binary logistic regression

applications where 2*log(likelihood ratio) may not be feasible. If you have evidence raising questions about the above, I'd like to know. spencer graves [EMAIL PROTECTED] wrote: th,ks for your help, i don't have this package on my R, do you know an other package that have

Re: [R] qq.loglogistic

with mine, I'd be pleased to hear about it. spencer graves Petr Mandys wrote: Hi, is there any similar function in R to S function qq.loglogistic, which produces a Q-Q plot? Thanks a lot Pete __ R-help@stat.math.ethz.ch mailing list

Re: [R] question regarding logit regression using glm

can't get from the factor coefficients. spencer graves Haibo Huang wrote: I got the following warning messages when I did a binomial logit regression using glm(): Warning messages: 1: Algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart

Re: [R] Using nonlinear regression

helpful reply -- if you don't find an answer in the course of preparing a question following the guide. spencer graves Mark Miller wrote: Hi, I have been trying to figure out how to use the nonlinear regression to fit the cumulative lognormal distribution to a number of data points I

Re: [R] some thoughts on outlier detection, need help!

question and inspire suggestions. TINSTAFL = There is no such thing as a free lunch (Heinlein, The Moon is a Harsh Mistress) spencer graves Weiwei Shi wrote: Dear listers: I have an idea to do the outlier detection and I need to use R to implement it first. Here I hope I can

Re: [R] exact goodness-of-fit test

-project.org/posting-guide.html;, show us your code and where you got stuck. spencer graves Christine Adrion wrote: Hello, I have a question concerning the R-function chisq.test. For example, I have some count data which can be categorized as follows class1: 15 observations class2: 0

Re: [R] p-values

Hi, Peter: Please see my reply of a few minutes ago subject: exact goodness-of-fit test. I don't know Rayner and Best, but the same method, I think, should apply. spencer graves Peter Ho wrote: HI R-users, I am trying to repeat an example from Rayner and Best A contingency

Re: [R] How to set the floating point precision beyond e-22?

. I rescaled each column by subtracting some number roughly close to the mean or median and dividing by something close to the standard deviation, and the problem went away. spencer graves Douglas Bates wrote: On 8/5/05, Roy Nitze [EMAIL PROTECTED] wrote: We have a problem inverting

Re: [R] interpolation function

I hope you don't take offense at anything said on this list. My philosophy about this is summarized in something I wrote last December that has since been immoratlized in the fortunes package: library(fortunes) fortune(Spencer Graves) Our great-great grandchilren as yet unborn

[R] Counterintuitive Simulation Results

)) # Thanks for your time in reading this. Best Wishes, Spencer Graves Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA [EMAIL PROTECTED] www.pdf.com http

Re: [R] Counterintuitive Simulation Results

Hi, Robert: Yes, I understand now. Thanks very much for your insight. Best Wishes, Spencer Graves McGehee, Robert wrote: Spencer, On the first iteration of your simulation, all of the Qp.t + z.t 0, so you're adding a vector of rep(4.5, 2) to a random

Re: [R] Help: how to stop the process when there is a mishandling

If you are using XEmacs, it's not that easy. I'm not completely clear on the protocol in that case, but I think you have to first click in the R window, then ctrl-G + ctrl-C + ctrl-C. Perhaps an Emacs wizard will enlighten us. spencer graves Duncan Murdoch wrote: On 8

Re: [R] Contour plot crest line

Have you studied the examples on the help pages for contour and persp? Also, if my memory is correct, there is some useful information this on Venables and Ripley, Modern Applied Statistics with S, 4th ed. (Springer). spencer graves Michael Kubovy wrote: Any suggestions

Re: [R] breaking command in command line in R for Mac Aqua

My favorite syntax error for such situations is [}, but as Peter said, that won't work if you are in the middle of a quote. Then one might try '[}'. If that failed, '[} should work. Thanks, Peter. spencer graves Peter Dalgaard wrote: Adaikalavan

Re: [R] problem with for()

Hint: help.search() - An Introduction to R - Simple manipulations numbers and vectors - ... spencer graves Peter Dalgaard wrote: Simone Gabbriellini [EMAIL PROTECTED] writes: Dear list, can someone tell me why this two pieces of code give me the same results? for(i in 0:5){ sum[i

Re: [R] can we manage memory usage to increase speed?

And you can identify inefficient code fairly easily taking snapshots from proc.time and computing elapsed time for sections of your code. spencer graves Tuszynski, Jaroslaw W. wrote: If you have a code that takes 2 weeks to run, than it might be a case of inefficient

Re: [R] How to hiding code for a package

you to encrypt your code. Talk with Insightful if you want to protect your intellectual property. spencer graves bogdan romocea wrote: There's something else you could try - since you can't hide the code, obfuscate it. Hide the real thing in a large pile of useless, complicated

Re: [R] cut.Date functionality for chron date/time objects

How about the following: set.seed(123) mychron - chron(sort(runif(10, 0, 10))) (breaks - chron(pretty(quantile(mychron [1] 01/01/70 01/03/70 01/05/70 01/07/70 01/09/70 01/11/70 spencer graves Sebastian Luque wrote: Gabor Grothendieck [EMAIL PROTECTED] wrote: Assuming, as in your

Re: [R] How to hiding code for a package

Dear Gary Wong: Are you trying to defeat the GNU License? If so, why not make it available as an S-Plus library? There, you don't have to provide the source; binary files will suffice. spencer graves bogdan romocea wrote: There's something else you could try - since

Re: [R] Minimum-chi-square estimation

about including a small, self-contained example that someone can copy from your email into R and try a few things before replying. It can help you find more things yourself AND increase the chances of getting useful reponses to posts. spencer graves -- Dragos Ilie wrote: Is there any

Re: [R] catching errors in a loop

the posting guide! http://www.R-project.org/posting-guide.html -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA [EMAIL PROTECTED] www.pdf.com http://www.pdf.com Tel: 408-938-4420 Fax: 408-280-7915

Re: [R] Forcing coefficents in lm(), recursive residuals, etc.

I know nothing about recresid, but does the following help: x - 1:4 set.seed(1) DF - data.frame(x=x, y=x+rnorm(4)) fit - lm(y~offset(x), DF) recresid(fit) [1] 1.2799320 0.7232336 3.4826581 spencer graves Rick Ram wrote: Hi all, Just to clarify, I know

Re: [R] Errbar()-function, cap and logarithmic scaling

, reproducible example, telling us which package you are using and which operating system. If someone can copy a few lines of code from your email and paste into R and see what you got, you will have a greater likelihood of getting a useful reply. spencer graves -- Ute Visser wrote

Re: [R] replace matrix values with names from a dataframe

-guide.html; and submit another question (if the process of working the posting guide does not itself provide enlightenment). spencer graves [EMAIL PROTECTED] wrote: Hi, I am looking for a way to replace matrix values with names from a dataframe. Let me do this by example: I have

Re: [R] Error Downloading Matrix Package

I've gotten that error message when trying to update a package I had attached. oops. (like do-it-yourself lobotomy.) Then I had to install.packages, because much of it got deleted and it wouldn't work any more. spencer graves Prof Brian Ripley wrote: Possible

Re: [R] fitting extreme value distribution

-project.org/posting-guide.html. It can increase the likelihood you will get a useful reply relatively quickly. spencer graves kunnavrv wrote: hi, rgev function gives me random deviates and I have a data set which I am fitting to an EVD,IS there a way I can plot both observed and ideal

Re: [R] Wishart Density

a useful reply. spencer graves Wang, Meihua wrote: Dear R users, I am doing MCMC using Metropolis-Hastings. My model is bivariate-log-normal and the prior for variance-covariance is wishart distribution. I am wondering if there are some simple codes about how to get the density

Re: [R] SETAR Estimation

(threshold autoregressive) produced 6 hits, the first of which was for wle.ar in package wle. hope this helps. spencer graves p.s. If this does not solve your problem, PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html. It can increase your chances

Re: [R] a class matrix with class ordered data

If you still are interested in help with your problem, PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html. It can increase the likelihood of you getting a useful reply. See other comments in line. spencer graves Werner Bier wrote

Re: [R] Convert quarterly data to monthly data

it is. If you need it, please PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html and try this list again.] spencer graves Haibo Huang wrote: Hi, I am new to use R, but can anyone tell me how to tranform quarterly data to monthly data? I know SAS has

Re: [R] Question on glm for Poisson distribution.

contr.treatment(outcome) and got something useful. If this is not adequate and you still could use help, PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html and try this list again. spencer graves Ghislain Vieilledent wrote: Good afternoon, I REALLY try

Re: [R] any package to fit such model?

Have you looked at ?lm, including the examples? The link in your email described a standard one-way, fixed effects ANOVA, and the lm help page includes a worked example for that. spencer graves ronggui wrote: I have search the internet but none are found.The lme function

Re: [R] find confounder in covariates

It was easier than I thought. To your example, just add one line: alias(f) Model : y ~ z1 + z2 + z3 Complete : (Intercept) z1B z1C z3 z2Bb 1 z2Cc 1 spencer graves Young Cho wrote: Hi, I was wondering if there is a way, or function

Re: [R] Clustered standard errors in a panel

Have you considered lmer in library(lme4)? If you are interested in this, you may want to check the article by Doug Bates in the latest R news, www.r-project.org - Documentation: Newsletter. spencer graves Thomas Davidoff wrote: I want to do the following: glm(y ~ x1

Re: [R] .gct file

Try ?read.table or args(read.table). Might skip=2 do what you want? spencer graves p.s. I routinely readLines(File, n=11) to see how many headers there are AND identify the sep character. Then I quantile(count.fields(File, ...)) to see if all records have the same number

Re: [R] definition of index.array and boot.return in the code for boot

Excellent question. Try 'getAnywhere(index.array)'. It's hidden in namespace:boot. Ditto for boot.return. spencer graves Obrien, Josh wrote: Dear R friends, I am reading the code for the function boot in package:boot in an attempt to learn how and where

Re: [R] Time Series Count Models

Mirner Klein (http://www.stat.sdu.dk/publications/monographs/m001/KleinPhdThesis.pdf and http://genetics.agrsci.dk/~bmk/IEKS.R). spencer graves Brett Gordon wrote: Thanks for the suggestion. Is such a model appropriate for count data? The library you reference seems to just be form

Re: [R] Confidence Intervals for Arbitrary Functions

is not even defined, being Inf-Inf.) In sum, you will likely get better answers in less time if you can find it politically acceptable in your professional efforts to work in decibels or logaritms, where ratios become differences. spencer graves Gabor Grothendieck wrote: On 7/16

Re: [R] Proportion test in three-chices experiment

problem. You may answer your own question in the course of working through that guide, and if you don't, the exercise could make it easier for someone else to suggest something you actually find useful. spencer graves Jonathan Baron wrote: On 07/17/05 20:12, Rafael Laboissiere wrote

<    2   3   4   5   6   7   8   9   10   11   >