[R] lmer t value for 3 levels of fixed factor

Hello, I just joined this list today, so am worried about proper protocol, but would like to post a question about lme4. In Baayen, Davidson, and Bates (2008), Mixed-effects modeling with crossed random effects for subjects and items, the authors describe steps for a Latin Square Design (p.

Re: [R] FIML using lavaan returns zeroes for coefficients

Thanks for the helpful explanation. As to your question, I sometimes use lavaan to fit univariate regressions simply because it can handle missing data using FIML rather than listwise deletion. Are there reasons to avoid this? BTW, thanks for the update in the development version. Andrew

[R] FIML using lavaan returns zeroes for coefficients

wishing to reproduce the problem, you can download a sample code file and data frame from the following two links. https://fds.duke.edu/db/aas/Sociology/grad/aam34/files/problem%20code.R https://fds.duke.edu/db/aas/Sociology/grad/aam34/files/sampledat.rdata Thanks for your help! Andrew Miles

[R] MuMIn - assessing variable importance following model averaging, z-stats/p-values or CI?

Andrew Robertson PhD student Centre for Ecology and Conservation University of Exeter, Cornwall Campus Tremough, Cornwall. TR10 9EZ UK Tel: 01326 371852 Email: ar...@exeter.ac.uk Web page: http://biosciences.exeter.ac.uk/staff/postgradresearch/andrewrobertson

Re: [R] Seeking pointers to various regression techniques with R?

) fits just the interaction between x and w without an intercept, and without the main effects for x and w lm(y ~ x/w - 1) I believe this fits the nested factor w inside of x Andrew Miles On Jun 5, 2012, at 4:58 PM, Michael wrote: I read your website but still don't know the difference

Re: [R] How do I obtain the current active path of a function that's being called?

, they are treated as any other object, which is to say that you can only have one function of the same name at a time. Hence whenever you call a source file to load in function A, the old function A gets overwritten. Andrew Miles On Jun 5, 2012, at 4:58 PM, Michael wrote: Hi all, How do I obtain

[R] Automated essay scoring by R

to know which package of R can be used to serve this purpose and how to do then. Grateful if you could give me the guidance. Thank you very much for your help. Regards, Andrew Cheung [[alternative HTML version deleted

[R] Higher log-likelihood in null vs. fitted model

. Andrew Miles [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] Higher log-likelihood in null vs. fitted model

= likelihood). AIC is a function of the likelihood but , as far as I know, likelihood is not a function of the AIC. Thanks for any insight. On Thu, May 31, 2012 at 9:26 AM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 12-05-31 8:53 AM, Andrew Miles wrote: Two related questions. First

[R] Axis mayhem in plot3d

xaxt=n to no avail Any thoughts? Andrew [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] odbcConnectExcelpel() fails to fetch all columns

to work with what I'm given. Kind regards to you both Andrew On 20/04/2012 23:58, andrija djurovic wrote: Hi. I use RODBC for importing Excel files quiet often and never got the similar problem. Have you tried with sqlQuery? z- odbcConnectExcel(./BBaselinePtQaires_apr2011.xls) BQ- sqlQuery(z

[R] odbcConnectExcel() fails to fetch all columns

and cbind to put it back together but that seems inelegant. The worksheet contains 134 rows, 70 columns and is in a spreadsheet that weighs in at 154 KB in total. Can you help unbaffle me? Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Multivariate Multilevel Model: is R the right software for this problem

. Andrew Miles Department of Sociology Duke University On Apr 6, 2012, at 9:48 AM, Eiko Fried wrote: Hello, I've been trying to answer a problem I have had for some months now and came across multivariate multilevel modeling. I know MPLUS and SPSS quite well but these programs could not solve

Re: [R] Multivariate Multilevel Model: is R the right software for this problem

() in the lme4 package, but has the advantage of being able to specify correlations between errors across time, the ability to control for time-invariant effects of time-invariant variables, and allows you to use the missing data maximum likelihood that comes in structural equation modeling packages. Andrew

[R] FIML in R

be interested in updates on this, and in your thoughts on this more generally. Also, please let me know if there is a forum better suited for this kind of discussion. Andrew Miles __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] fixed effects linear model in R

Allison's booklet for more details on this method. Allison, Paul D. 2009. Fixed Effects Regression Models. Los Angeles, C.A.: Sage. Andrew Miles On Feb 7, 2012, at 5:00 PM, caribou...@gmx.fr wrote: Dear R-helpers, First of all, sorry for those who have (eventually) already received

Re: [R] Help need

in seq(0,0.5,0.1)) { sigmasqaured - 1 i = complex(real = 0, imaginary = 1) spectrum - c(spectrum, (sigmasqaured)/(abs(1-2.7607*exp(2*pi*i*f)+3.8106*exp(4*pi*i*f)-2.6535*exp(6*pi*i*f)+0.9258*exp(8*pi*i*f))^2)) } spectrum Andrew Miles On Feb 7, 2012, at 4:08 PM, Jaymin Shah

Re: [R] Problem with GMT+/- time zones

west of GMT have a positive sign and those east have a negative sign in their name (e.g Etc/GMT-14 is 14 hours ahead/east of GMT.) Thanks, Andrew On 2/02/2012, at 20:46 , Jeff Newmiller wrote: Should has nothing to do with it. That is the way the Olsen tz database works. See

[R] Problem with GMT+/- time zones

I'm struggling with time zone version when expressed as hours offset from GMT. Can anyone confirm that the behaviour below is incorrect? It seems that the GMT offsets are backwards: format(as.POSIXct(2011-05-23 17:23:00, tz=Europe/London),tz=America/New_York,usetz=T) [1] 2011-05-23 12:23:00

[R] Change line colors based on data values in Lattice

Say I have x = seq(1, 40) y = seq(1, 40)^2 g = rep(c(A, B, A, B), each = 10) xyplot(y~x, groups = g, type = 'l') Is there a way to have a single line and have it be one color when g == A, and another when g == B? Thanks! Regards, Andrew Liu Fixed Income | Investment Fellow T. Rowe Price 100

Re: [R] Summing rows by years (each time separately)

Thanks for your answer! There is one thing I could not find in aggregate(): I want to have it a sum for each group of 2008,2009 and 2010. In aggregate() I can sum all the rows that have a rowname 2008, all the rows that have a rowname 2009 and all the rows that have a rowname 2010. But I want

[R] splitting strings effriciently

$IPC[i] -lst[[3]] data$IPD[i] -lst[[4]] rm(lst) } Andrew Andrew Roberts Children's Orthopaedic Surgeon RJAH, Oswestry, UK [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] Summing rows by years (each time separately)

10 *2010* 6 5 3 Could You help me with that? Yours, Andrew Gaska -- View this message in context: http://r.789695.n4.nabble.com/Summing-rows-by-years-each-time-separately-tp4276428p4276428.html Sent from the R help mailing list archive at Nabble.com

Re: [R] splitting strings efficiently

Thanks Enrico Jim, The following finished the job in under a minute! res - unlist(strsplit(data[[ComputerName]], \\.)) ii - seq(1, nrow(data)*4, by = 4) data$IPA -res[ii] ## A data$IPB -res[ii+1] ## B data$IPC -res[ii+2] ## C data$IPD -res[ii+3] ## D Andrew On 08/01/2012 13:11, Enrico

[R] Simulate binary time series with known Hurst exponent

is simulate sequences of varying length, probability of the binomial distribution (say from 25-75 % successes), and H. I would very much appreciate any advice on this matter. Thank you, Andrew Andrew MacIntosh Kyoto University __ R-help@r-project.org

Re: [R] Memory getting eaten up with XML

Today I tried the code on a MacBook and experienced the same problem. Which makes me think there is something wrong with the way I am trying to free up the memory...? Andrew Andrew Gormley wrote Hi all. I have an issue that I cannot resolve. I am trying to read in lots of data

[R] Memory getting eaten up with XML

, the memory usage remains the same as it was when I first read it in... Any ideas/solutions? I am using a 32bit version of R 2.14.0 on windows XP, and the latest version of XML (3.6.1). Many thanks Andrew (apologies for the large footer my work appends to all my emails...) Please consider

[R] Lattice graph strips and axes

= list(y = c(0.5),x = c(0.5)),scales = list(x = list(alternating = 1, tck = c(1,0)), y = list(alternating = 1, tck = c(1,0))), xlab = count (n), ylab = titer) useOuterStrips(x) Kind regards Andrew McFadden MVS BVSc | Veterinary Epidemiologist, Investigation and Diagnostic Centre | Biosecurity New

[R] Lattice graph help

) Kind regards Andrew McFadden MVS BVSc | Veterinary Epidemiologist, Investigation and Diagnostic Centre | Biosecurity New Zealand Ministry of Agriculture and Forestry | 66 Ward St, Wallaceville | PO Box 40 742 | Upper Hutt | New Zealand Telephone: 64-4-894 5611 | Facsimile: 64-4-894 4973| Mobile

Re: [R] Compare clustering solutions to a correct one

It sounds like you want to do supervised classification, so maybe a supervised classification algorithm would be more appropriate? Consider logistic regression, rpart, ctree, earth, etc. Andrew -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r

Re: [R] R in batch mode packages loading question

Try the fork() function in the package multicore (if your system supports it) Andrew -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of PALMIER Patrick (Responsable de groupe) - CETE NP/TM/ST Sent: Monday, November 07, 2011 3:49 AM

Re: [R] install.packages problem

. On Fedora or RedHat you would do something like 'sudo yum -y install curl-devel' and on Ubuntu it may be 'sudo apt-get install curl-dev' Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Decision tree model using rpart ( classification

, and without this model that doesn't finish fitting you should already be able to assign a decile to every customer. Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

[R] bug calculating ROC with caret and earth?

and specificity would also not vary, but they do vary. Also, I verified nprune is not too high. I am attaching sample output from R 2.14.0 on Windows 7 64-bit with earth 3.2 and caret 5.07. I don't have this problem with caret and ctree. Andrew R version 2.14.0 (2011-10-31) Copyright (C) 2011

[R] text orientation in persp

I am doing a 3D plot using persp {graphics} but cannot figure out how to change the text direction on the z-axis. The default is for it to read down. Also, I cannot figure out how to change the position of the text itself. Many thanks Andrew Please consider

[R] Sample size calculations for one sided binomial exact test

a prespecified value, but cannot find any functions for computing sample size. do any exist? Thanks, Andrew [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] plotting issues with PCA

cant get the lines to colour code the same way as the points did ordiellipse(graph,site.codings,kind=sd,conf=0.90,draw=polygon) I see there is a command called show.groups but I cant work out how to use it to access the plotcolor file. Any help appreciated. -- Andrew Halford Ph.D Associate

Re: [R] Wide to long form conversion

Take a look here. http://stackoverflow.com/questions/2185252/reshaping-data-frame-from-wide-to-long-format Andrew Miles Department of Sociology Duke University On Oct 6, 2011, at 4:28 PM, Gang Chen wrote: I have some data 'myData' in wide form (attached at the end), and would like to convert

[R] predict.rma (metafor package)

model3 newdat-expand.grid(c(1,2,3,4,5),c(1,2,3)) predict(model1,newmods=c(1,2,3,4,5)) predict(model2,newmods=c(1,2,3)) predict(model3,newmods=newdat) -- Andrew Beckerman Sent with Sparrow (http://bit.ly/sigsprw) [[alternative HTML version deleted

Re: [R] Lmer plot help

* the IVs? I have 5 voices (and 4 questions), but it lists the effects of 4 voices/3questions, similarly if I run summary (model). Best Andrew E. MacFarlane PhD student New Zealand Institute of Language, Brain and Behaviour University of Canterbury | Private Bag 4800 Christchurch

[R] Fitting the negative hypergeometric distribution

, as estimated by another program, for the start values... Maybe I am doing something obviously wrong here, or perhaps there is a better solution. Many thanks, Andrew Wilson __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] excluding models during dredge and model averaging in MuMIn

it if somebody had any ideas about how to tackle this. Also, Cheers, Andrew Andrew JJ MacIntosh, D. Sc. Center for International Collaboration and Advanced Studies in Primatology Kyoto University Primate Research Institute Kanrin 41-2 Inuyama Aichi Japan 484-8506 website: http://www.wix.com/andyjmac/andrew

Re: [R] Calculating p-value for 1-tailed test in a linear model

...@comcast.net] Sent: Saturday, August 20, 2011 6:02 PM To: Andrew Campomizzi Cc: r-help@r-project.org Subject: Re: [R] Calculating p-value for 1-tailed test in a linear model On Aug 19, 2011, at 6:20 PM, Andrew Campomizzi wrote: Hello, I'm having trouble figuring out how to calculate a p

Re: [R] Calculating p-value for 1-tailed test in a linear model

but it doesn't show the understanding. -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Monday, August 22, 2011 9:12 AM To: Andrew Campomizzi Cc: r-help@r-project.org Subject: Re: [R] Calculating p-value for 1-tailed test in a linear model On Aug 22, 2011, at 9:44

Re: [R] Calculating p-value for 1-tailed test in a linear model

. Thanks, Andy - Original Message - From: Rolf Turner rolf.tur...@xtra.co.nz To: Andrew Campomizzi acampomi...@tamu.edu Cc: r-help@r-project.org Sent: Friday, August 19, 2011 9:18:07 PM GMT -06:00 US/Canada Central Subject: Re: [R] Calculating p-value for 1-tailed test in a linear model The r

[R] Calculating p-value for 1-tailed test in a linear model

Hello, I'm having trouble figuring out how to calculate a p-value for a 1-tailed test of beta_1 in a linear model fit using command lm. My model has only 1 continuous, predictor variable. I want to test the null hypothesis beta_1 is = 0. I can calculate the p-value for a 2-tailed test using

[R] plot 3d info in 2d

axes would be larger while sites that are actually further away in 3 d space would be proportionally smaller. any help/advice appreciated Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948 [[alternative HTML version deleted

[R] reflecting a PCA biplot

appreciated. cheers Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Games-Howell post-hoc testing

and possibly unequal variance among samples. -- Andrew D. Steen, Ph.D. Center for Geomicrobiology, Aarhus University Ny Munkegade 114 8000 Århus C Denmark andrew.st...@biology.au.dk __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] xlsx error

get an error message: Error in .jcall(row[[ir]], Lorg/apache/poi/xssf/usermodel/XSSFCell;, : method getCell with signature (I)Lorg/apache/poi/xssf/usermodel/XSSFCell; not found Anyone have any idea how to fix this issue? Thanks, Andrew Winterman [[alternative HTML version deleted

[R] Plotting question

Sincerely Andrew McCulloch __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] list of data.frames

Hey there, I've got a for loop in my code which I'm having difficulty seeing how to avoid. Here's the code, with some sample data. You can safely copy and paste into the interpreter if you want. #I have some data frames. a - as.data.frame(matrix(runif(50), nrow = 10, ncol = 5)) b -

[R] S4 group generic with additional argument

: i486-pc-linux-gnu (32-bit) and     R version 2.13.0 (2011-04-13)     Platform: i386-pc-mingw32/i386 (32-bit) Many thanks, Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Off-topic: (Simple?) Random Sampling when n is a random variable

On Tue, Jun 14, 2011 at 11:02:52AM +1000, Andrew Robinson wrote: Hi everyone, I'm involved in a discussion with a colleague. He suggested a sample design for a finite-sized process that (to all intents and purposes) involves tossing a coin and examining the unit if the coin shows Heads

Re: [R] Off-topic: (Simple?) Random Sampling when n is a random variable

Thanks Greg! Andrew On Tue, Jun 14, 2011 at 04:13:52PM -0600, Greg Snow wrote: This sounds like what is called domains in survey sampling (possibly other names, but that is what I learned it as). The idea is that you take a random sample (or the population) then ask a question to determine

[R] Off-topic: (Simple?) Random Sampling when n is a random variable

some citations that provide guidance either way? Thanks for any assistance! Andrew -- Andrew Robinson Program Manager, ACERA Department of Mathematics and StatisticsTel: +61-3-8344-6410 University of Melbourne, VIC 3010 Australia (prefer email) http

[R] ordiellipse

assuming that 2 points is not enough to perform the relevant calculations here, however I would like to plot one if I could, if only for the sake of pictorial consistency. Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948

Re: [R] gsub() issue...

Hi Charles, It's not clear to me what you mean by doesn't work. test - Interesting 1\nPoint\n cat(test) Interesting 1 Point test1 - gsub(ing 1\nP,ing 3\nP, test) cat(test1) Interesting 3 Point Cheers Andrew On Tue, May 17, 2011 at 10:45:31AM +0200, Thibault Charles wrote: Hello R

[R] Dealing with null values Aggregate function

)) aggregate(x = testDF, by = list(testDF$v3), FUN = mean) For this example I would like to get a mean value for red (v1=4.5; v1=5) and blue (v1=49.5; v2=55) Thank you in advance. Andy Andrew McFadden MVS BVSc Incursion Investigator Investigation Diagnostic Centres - Wallaceville Biosecurity New

[R] mlogit error message Error in X[omitlines, ] - NA : subscript out of bounds

that an answer to this question will be documented. Andrew Miles I am using the mlogit packages and get a data problem, for which I can't find any clue from R archive. code below shows my related code all the way to the error

Re: [R] functions pandit and treebase in the package apTreeshape

Thanks! As noted by an earlier poster, these functions seem to have been removed from the package. It might be worth contacting the maintainer to ask if any equivalents are available ... Cheers Andrew On Sat, May 07, 2011 at 09:56:49AM +0200, Arnau Mir wrote: El 07/05/11 09:49, Arnau Mir

Re: [R] nls problem with R

Apologies, but I don't see a question here ... am I missing something obvious? Andrew On Thu, May 05, 2011 at 01:20:33AM -0700, sterlesser wrote: ID1 ID2 t V(t) 1 1 0 6.053078443 2 1 0.3403 5.56937391 3 1 0.4181 5.45484486 4

Re: [R] functions pandit and treebase in the package apTreeshape

Hi Arnau, please send the output of sessionInfo() and the exact commands and response that you used to install and load apTreeshape. Cheers Andrew On Thu, May 05, 2011 at 04:42:58PM +0200, Arnau Mir wrote: Hello. I'm trying to use the functions pandit and treebase. They are in the package

Re: [R] Vermunt's LEM in R

Hi David, you might have more luck with your request if you tell us what Vermunt's LEM *does*, and provided some links to introductory reading material ... Cheers Andrew On Thu, May 05, 2011 at 03:34:02PM +, David Joubert wrote: Hello- Does anyone know of packages that could emulate

Re: [R] two-way group mean prediction in survreg with three factors

of predicting A and B is just to remove it from the model, or alternatively to make it a random effect. Neither idea seems rock solid at this point. Cheers Andrew On Thu, May 05, 2011 at 09:37:15AM -0400, Pang Du wrote: Oops, I hope not too. Don't know why I had the brackets around B+C. My

Re: [R] Installing rgdal in R: correct -configure flags for GDAL install on Linux Redhat

suggested below. I hope that this helps. Andrew On Thu, May 05, 2011 at 05:23:56PM -0800, Katrina Bennett wrote: Hi, I'm installing rgdal but I keep having failures because I have not been able to find a good source of information for the correct configuration settings when installing GDAL. My

Re: [R] nls problem with R

help you. Cheers Andrew On Tue, May 03, 2011 at 09:08:03PM -0700, sterlesser wrote: the original data are V2 =c(371000,285000 ,156000, 20600, 4420, 3870, 5500 ) T2=c( 0.3403 ,0.4181 ,0.4986 ,0.7451 ,1.0069 ,1.553) nls2=nls(V2~v0*(1-epi+epi*exp(-cl*(T2-t0))),start=list(v0=10^7,epi=0.9,cl=6.2

Re: [R] problem with package adapt for R in Mac

Hi, Is there such a package? I can't find it on CRAN. Can you let us know exactly how you tried to install it, and what the error message was (if any)? Cheers Andrew On Wed, May 04, 2011 at 01:29:37AM -0300, Mat?as Ram?rez Salgado wrote: Hi, How i can install the package adapt in some

[R] what happens when I store linear models in an array?

(lms.ASP[1,1]) NULL It seems like something less than the actual linear model object is being stored in the array, but I don't understand what's happening, or how to easily batch-extract parameters of linear models. Any advice? Andrew D. Steen, Ph.D

Re: [R] select value from a column depending on a value in another column

Try subset(). Andrew On Wed, May 04, 2011 at 02:01:44AM -0700, tornanddesperate wrote: Hi everybody I couldn't find the solution to what must be quite a simple problem. Maybe you can help? treatment session period stage wage_accepted market 1 11 1

Re: [R] what happens when I store linear models in an array?

Hi Andrew, try fitted(lms.ASP[1,1][[1]]) Cheers Andrew On Wed, May 04, 2011 at 01:49:45PM +0200, Andrew D. Steen wrote: I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all

Re: [R] what happens when I store linear models in an array?

Thanks all, this is very helpful. --Andrew Steen -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, May 04, 2011 4:35 PM To: Andrew D. Steen Cc: r-help@r-project.org Subject: Re: [R] what happens when I store linear models in an array

Re: [R] help with the maxBHHH routine

, then simplify it futher. If you get the error with the simplest possible problem, then share it. If you don't , then try to figure out what the changes were that resolved the problem, and scale those back up to your original problem. Does that make sense? Cheers Andrew On Thu, May 05, 2011 at 03:22

Re: [R] Problems saving ff objects

I wonder if this question should be directed to the package maintainer? Best wishes, Andrew On Wed, May 04, 2011 at 02:31:51PM +0100, Jannis wrote: Just did some more testing.May the problem be due to the fact that I am using a windows machine? I just ran the same code on a Linux machine

Re: [R] Panels order in lattice graphs

Hi Cristina, you can probably hack your own solution using the index.cond argument. Cheers Andrew On Wed, May 04, 2011 at 04:50:53PM +0100, Cristina Silva wrote: Hi all, In lattice graphs, panels are drawn from left to right and bottom to top. The flag as.table=TRUE changes to left

Re: [R] fGarch

. Maybe even solve it. Even if you don't, then you enable us to make a much better contribution, because we can actually try out our suggestions before sending them. With what you sent here, all we can do is speculate. Cheers, Andrew On Wed, May 04, 2011 at 04:33:12PM -0400, Paul Ossenbruggen

Re: [R] two-way group mean prediction in survreg with three factors

its levels. What does that summary really mean? Can you meaningfully average across the levels of a predictor that is included in the model as a main and an interaction term? Best wishes Andrew On Wed, May 04, 2011 at 12:24:50PM -0400, Pang Du wrote: I'm fitting a regression model for censored

Re: [R] split character vector by multiple keywords simultaneously

A hack would be to use gsub() to prepend e.g. XXX to the keywords that you want, perform a strsplit() to break the lines into component strings, and then substr() to extract the pieces that you want from those strings. Cheers Andrew On Wed, May 04, 2011 at 04:08:40PM -0700, sunny wrote: Hi. I

Re: [R] Uniform Gaussian Kernel

I can't honestly say that I grasp what you're trying to do, but that said, I wonder if the curve() function will help you? Cheers Andrew On Wed, May 04, 2011 at 07:30:20AM -0700, blutack wrote: I have a vector with lots of different numbers. I need to make a graph showing the Uniform

Re: [R] bivariate linear interpolation

The one that gives results that you trust and uses algorithms that you understand! Cheers Andrew On Wed, May 04, 2011 at 12:52:03PM +, Halldór Björnsson wrote: Hi, I have three matrices (X,Y,P) with the same dimension. The X,Y grid is regular and I want to perform linear interpolation

[R] Overlapping x axes using Lattice

, stack=TRUE,type = count, xlab=VNT,rot=c(180,180),draw=FALSE) Thank you in advance. Andy Andrew McFadden MVS BVSc Incursion Investigator Investigation Diagnostic Centres - Wallaceville Biosecurity New Zealand Ministry of Agriculture and Forestry Phone 04 894 5600 Fax 04 894 4973 Mobile 029 894

Re: [R] Unexp. behavior from boot with multiple statistics

Your interpretation of what the output is supposed to look like is actually correct. Take a look at the estimates of the bias in the BootStrap Statistics. You will see that they are the same as the difference between the location of colMeans of t and t0. I hope that this helps, Andrew On Tue

Re: [R] help with the maxBHHH routine

I suggest that you provide some commented, minimal, self-contained, reproducible code. Cheers Andrew On Wed, May 04, 2011 at 02:23:29AM +0530, Rohit Pandey wrote: Hello R community, I have been using R's inbuilt maximum likelihood functions, for the different methods (NR, BFGS, etc). I

Re: [R] delete excel id automatically generated

Try the function rownames() Andrew On Tue, May 03, 2011 at 03:29:37AM -0700, agent dunham wrote: Dear community, I uploaded an excel with read.xls. My xls file actually have a column which is an id, (plot is the id) : plot height area 347.6 5.4 853.2 4.1 895.4

Re: [R] na.omit - Is it working properly?

in commented, minimal, self-contained, reproducible code ... Cheers Andrew On Tue, May 03, 2011 at 12:18:03PM -0700, Kalicin, Sarah wrote: I have a work around for this, but can someone explain why the first example does not work properly? I believed it worked in the previous version of R

Re: [R] Watts Strogatz game

Hi, I have no familiarity with these functions --- I see that they are not in base R --- so I suggest that at very least you identify the package that you are using. Better would be to contact the package maintainer directly. Sometimes maintainers do not read R-help. Cheers Andrew On Tue

Re: [R] subseting data

I wonder if grep() will help you? Cheers Andrew On Mon, May 02, 2011 at 11:03:52AM +0200, Matev? Pavli? wrote: Hi, Is it possible (i am sure it is) to subset data from a data.frame on the basis of SQL LIKE operator. I.e., i would like to subset a data where only values which

Re: [R] UNIX-like cut command in R

Hi Mike, try substr() Cheers Andrew On Mon, May 02, 2011 at 03:53:58PM -0500, Mike Miller wrote: The R cut command is entirely different from the UNIX cut command. The latter retains selected fields in a line of text. I can do that kind of manipulation using sub() or gsub

Re: [R] Optimization - n dimension matrix

Hello, optim() works for more than one dimension. You might also find this page helpful: http://cran.r-project.org/web/views/Optimization.html Cheers Andrew On Mon, May 02, 2011 at 12:41:19PM -0700, petrolmaniac wrote: Dear all, I am facing the following problem in optimization: w

Re: [R] Help with coloring segments on a plot

more luck. Best wishes Andrew On Mon, May 02, 2011 at 02:26:16PM -0400, Paul Davison wrote: Hi. I need a very short piece of help regarding colouring segments plotted on a graph. When I am plotting segments for the graph, I am using red and darkgreen for the values 1 and 2 respectively

Re: [R] Simulation Questions

+ beta_1 * x + rnorm(n, 0, sigma^2) and use that as the basis of generating random numbers. Not sure how to interpret the second question ... Cheers Andrew On Sun, May 01, 2011 at 12:33:41AM -0400, Shane Phillips wrote: I have the following script for generating a dataset. It works like a champ

Re: [R] Lasso with Categorical Variables

, for X1, if lasso drops the coefficient for levels A and B, but not C and D, does this mean that X1 should be included in the model? It means that X1 should be recoded to be C, D, and the rest. Cheers Andrew Thanks. On Mon, May 2, 2011 at 2:47 PM, David Winsemius dwinsem...@comcast.net

[R] Urgent: conditional formula for nls

I have data vectors x and y both with 179 observations. I'm trying to fit a nonlinear model with five parameters using nls. The formula is only defined within a range of x-values, it should be zero otherwise, thus my attempted use of ifelse: df-data.frame(x,y)

Re: [R] using tapply with multiple variables

, other target variables listed here, ...), by = list(Subject = Subject, Group = Group), FUN = function(x) sum(x == 'C'))) Cheers Andrew On Sat, Apr 30, 2011 at 10:03:24PM -0700, Dennis Murphy wrote: Hi

Re: [R] Different results of coefficients by packages penalized and glmnet

--- that would be best. If not, you should write code to generate an example dataset, or be prepared to share your own data. I hope that this helps, Andrew On Sun, May 01, 2011 at 05:01:54PM +0800, zhu yao wrote: Dear R users: Recently, I learn to use penalized logistic regression. Two packages

Re: [R] Analysis and graphics by groups

to construct a minimal, reproducible, commented example to help us explain what you need to do. I don't have the gapply function and I don't know which package it is in (perhaps you could provide that information next time) so I can't help more than that. Andrew On Fri, Apr 29, 2011 at 04:31:51PM

Re: [R] Variance

easier to work with commented, minimal, self-contained, reproducible code. Please consider submitting that with future questions. I hope that this helps, Andrew. On Thu, Apr 28, 2011 at 05:58:24PM -0400, Dat Mai wrote: I'm trying to find the variance of various outputs in a matrix: for(l in 2

Re: [R] Problem installing package sp in R 2.13.0

Hi Arnaud, the error is telling you that you don't have the make command. This might be because you haven't installed the necessary software to compile R packages. I suggest that you check the FAQ for Macintosh to see how to do that. Best wishes Andrew On Thu, Apr 28, 2011 at 01:30:58PM

Re: [R] Generating a best fit line for non linear data

I think that you probably need to provide the x values to nls. Try, for example, fit - nls(species ~ a *(1 - exp(-b*samples)),start = list(a = 27, b = .15)) I hope that this helps, Andrew On Thu, Apr 28, 2011 at 01:56:43PM -0700, BornSurvivor wrote: I have the following data set, and I have

Re: [R] Problem loading quantreg

/tools/ I hope that this helps, Andrew On Thu, Apr 28, 2011 at 03:15:41PM -0700, derekverley wrote: Hi all, I'm trying to load the quantreg package but keep running into problems no matter which method I have tried. Does anybody know what this error (below) means in plain language and what

<    1   2   3   4   5   6   7   8   >