[R] Plotting a curve for a Holling Type III Functional Response

Hey, So I have a scatter plot and I am trying to plot a curve to fit the data based on a Holling Type III functional response. My function is this: nll2<-function(a,b) { conefun<-(a*DBH^2)/(b^2+DBH^2) nlls2<-dnbinom(x=cones ,size=DBH, mu=conefun,log=TRUE) -sum(nlls) } and my plot is this:

[R] File system issue since upgrading macOS to High Sierra

Hi, Since upgrading macOS High Sierra, R will stop responding, requiring me to force R to quit, when I am trying to open a script using File -> Open Document. I don't have any problems opening files by double-clicking them in Finder. Sometimes R gives me the message in the attached screenshot:

[R] Storing tableGrobs in a list

Hi, I'm trying to generate tableGrobs in a loop, store them in a list so I can use it in a call to gtable_combine(). L1<-list() for (i in seq( ... )) { L1[i] <-tableGrob( ... ) } gtable_combine(L1, along=1) On the assignment inside the loop, I get "number of items to replace is not a multiple

Re: [R] Storing tableGrobs in a list

gt;L1<-vector( "list", N ) >for (i in seq.int( N )) { > L1[[i]] <-tableGrob( ... ) >} > >PS Post using your email program "plain text" mode... HTML gets >stripped anyway and that often leads to partial corruption of your >message. Read the Posting Guide. &

Re: [R] Storing tableGrobs in a list

mpt someone to dig a little. > >On June 14, 2018 8:07:19 AM HST, Stats Student > wrote: >>Thanks for the replies. Wasn't aware that Gmail on Android sent HTML >by >>default, apologies. >> >>Storing the tableGrob-s in a list worked but for some reason >>

[R] Adding lines to the page

Hi, I'm looking for a way to add lines to a report. To be clear, I don't want to add lines to any specific plot, but instead to add line(s) to the page itself - e.g. add a line to the footer area, above the actual footer text. Any thoughts on how to do this? Many thanks. __

Re: [R] Adding lines to the page

Rnw files) is very powerful, but there are many >other tools as well (e.g. bookdown) depending on your preferences. > >[1] https://cran.r-project.org/web/views/ReproducibleResearch.html > >On June 27, 2018 1:53:58 PM PDT, Stats Student > wrote: >>Hi, I'm looking for a way

Re: [R] Adding lines to the page

; > > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along and > sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > On Wed, Jun 27, 2018 at 3:56 PM, Stats Student <

[R] scale_y_continuous with sec.axis

Hi, I'm using scale_y_continuous with sec.axis and it's doing what I need but I don't understand how it picks which of the two series becomes the secondary. Does anyone have any insight into this? Thanks! __ R-help@r-project.org mailing list -- To U

[R] Question about correlation

Greeting. Dear Mr/Mrs/Miss, OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05 Above show the first 6 data sets, may I ask that the reason of R show the error like "Error in cor(data) : 'x' must be numeric" ? And how

[R] Conduct Network Analysis

Greeting. Dear Mr/Mrs/Miss, I want to create a network by using R but I only have a table that contain OTU ID and the abundance value of two samples ONLY. Isn't possible? If can, which package can be used? Greatly appreciated to any suggestions and helps. Thank you. Best regards, Kang Chin Yi

[R] Changing PDF orientation midstream

Hi, I'm wondering whether it is possible to change the orientation of the PDF in the middle of the document. In other words, pages 1,2,3 - portrait, pages 4,5 - landscape, etc. This is how I call it - pdf (file, paper="US") or USr for landscape Thanks! _

[R] Ordering of facet_wrap() panels

Hi, I am generating multiple charts with facet_wrap() and what what I see, R/ggplot sorts the panels by the facet variable. So adding an index to the facet variable (1 - bucket, 2 - bucket, etc) does solve the sorting issue but it's ugly. I also read this post which, if I understand correctly,

Re: [R] Ordering of facet_wrap() panels

nique(f)) >[1] c b a c b a >Levels: c b a ## not your desired order > >Again, please consult the docs and perhaps a tutorial or two as >necessary. > >-- Bert > > > >On Wed, Aug 15, 2018 at 8:22 AM, Stats Student > >wrote: > >> Many thanks, Bert. >&g

[R] Survival questions

Hi - I am using the survfit() function to produce Kaplan-Meier survival curves for several different groups. survfit (Surv() ~ cohort, data=d) Everything works fine, but I'd like to do something different. I have a 12 month survival curve (base) and I also have two survival curves based on some

[R] Addressing results of modFit() of the FME package

Hello, I don´t manage to address the results I get from modFit(). Would be really nice if you could help me to do so. This is my loop: for (i in 1:10) { CO2 <- CO2.list[[i]] #Beobachtungen Obs <- data.frame("x=t"=numeric(),"y=Respi"=numeric()) Obs <- cbind(CO2[,1],CO2[,3]) colnames(Obs) <- c(

Re: [R] help me with holt-winter model

Please help me with this i need to submit my thesis . Thanks In advance -- View this message in context: http://r.789695.n4.nabble.com/help-me-with-holt-winter-model-tp2295552p2298464.html Sent from the R help mailing list archive at Nabble.com. _

[R] Please advise acf and pacf in order to determine order of Arima

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

Re: [R] Odp: Help me with prediction in linear model

Thanks Murphy and pikal, I need another help,for fitting first fourier transformation ,i used following thing .Please advise on this beer_monthl has 400+ records EXample: > head(beer_monthly) beer 1 93.2 2 96.0 3 95.2 4 77.1 5 70.9 6 64.8 time<-seq(1956,1995.2,length=length(beer_monthly)

[R] density plot with frequency units

I would like to create a kernal density plot, but rather than show density units on the vertical axis I would like frequencies. I know histograms do this but I don't want the bars, just the density curve. Thanks! __ R-help@r-project.org mailing list h

Re: [R] density plot with frequency units

(Solar/Batteries O.O#. #.O#. with > /Software/Embedded Controllers) .OO#. .OO#. rocks...1k > --- > Sent from my phone. Please excuse my brevity. > > r student wrote: >> >> I would like to create a kern

[R] t-test on a data-frame.

Dear R-helpers, In a data frame I have 100 securities,monthly closing value,from 1995 to present,which I have to 1. Sampling with replacement,make 50 samples of 10 securities each,each sample hence will be a data frame with 10 columns. 2. With uniform probabilty,mark a month from 2000 onwards as

[R] help with algorithm

I'm wondering if anyone can give some basic advice about how to approach a specific task in R. I'm new to R but have used SAS for many years, and while I can muscle through a lot of the code details, I'm unsure of a few things. Specific questions: If I have to perform a set of actions on a grou

Re: [R] help with algorithm

roups > > dat <- mtcars > dat$gm <- with(mtcars, ave(mpg, cyl, FUN = mean)) I tried, but I think the NAs are giving me trouble. tmp <- oh tmp$GM <- with(oh, ave(FINCP, PUMA, FUN=mean)) summary(tmp$GM) Min. 1st Qu. MedianMean 3rd Qu.Max.NA's

[R] density plot for weighted data

I'm trying to create a density plot using census data, where the weights don't sum to 1. >plot(density(oh$FINCP,weights=oh$PWGTP)) Warning message: In density.default(oh$FINCP, weights = oh$PWGTP) : sum(weights) != 1 -- will not get true density How would I go about doing this? Thanks!

Re: [R] density plot for weighted data

Like below? plot(density(oh$FINCP,weights=oh$PWGTP/sum(oh$PWGTP))) On Tue, Aug 2, 2011 at 10:06 AM, David Winsemius wrote: > > On Aug 2, 2011, at 12:51 PM, r student wrote: > >> I'm trying to create a density plot using census data, where the >> weights don&

[R] Error with cForest

All -- I have been trying to work with the 'Party' package using R v2.15.1 and have cobbled together a (somewhat) functioning code from examples on the web. I need to run a series of unbiased, conditional, cForest tests on several subsets of data which I have made into a loop. The results ide

[R] lme4 2 factor factorial model with random factors

Using lme4 how does one define a 2 factor factorial model with both factors being random? Specifically I am just trying to recreate the results from Montgomery's Design of Experiments book (7th edition), example 13.2. In this example there are 2 random factors and I want to include the interaction

[R] LogLikelihood of a Distribution Given Fixed Parameters

I'm trying to figure out if there is a way in R to get the loglikelihood of a distribution fit to a set of data where the parameter values are fixed. For example, I want to simulate data from a given alternate lognormal distribution and then I will fit it to a lognormal distribution with null param

Re: [R] LogLikelihood of a Distribution Given Fixed Parameters

heers, > > Rolf Turner > > > On 25/04/14 07:27, Jacob Warren (RIT Student) wrote: > >> I'm trying to figure out if there is a way in R to get the loglikelihood >> of >> a distribution fit to a set of data where the parameter values are fixed. >> For ex

[R] Unbalanced Design Power Analysis

I have an unbalanced design I would like to run a power analysis on. What I have been able to find has pointed me to using the pwr.f2.test function as described below. My problem is that I don't know how to appropriately define the numerator and demoninator df. If someone can help here is some mor

[R] Need help for R install

Dear R committee: I am Renzhi, Ph.D student in computer science in the University of Missouri. I have one question for you. I try to install R in the linux server, but I don't have the root permission, is there any way to install the R locally? Thank you very much for helpi

[R] Help With Permutations

I have a problem with permutations functions in R I just started using R so maybe it is an easy question I need to obtain all the 9.somthingExp 157 permutations that can be given from the number from 1 to 100 I wrote the following commands: > library(gregmisc) >options(expressions=1e5)

Re: [R] error in data.farme--duplicate row.names error

ael Dewey [mailto:li...@dewey.myzen.co.uk] Sent: Thursday, May 19, 2016 11:05 AM To: Rees, Lisa Marie (MU-Student); r-help@r-project.org Subject: Re: [R] error in data.farme--duplicate row.names error Dear Lisa What does Bmat contain? Perhaps try table(table(Bmat)) and see if any entries are greate

Re: [R] error in data.farme--duplicate row.names error

Don, Thank you for your helpful response. At this point, I do believe it is a package error and I have contacted the developer. Thanks, Lisa -Original Message- From: MacQueen, Don [mailto:macque...@llnl.gov] Sent: Thursday, May 19, 2016 4:27 PM To: Rees, Lisa Marie (MU-Student); r

[R] GGplot annotate by facet

Hi, I have a dataframe which I need to plot in ggplot2 it looks like this :- head(nodelta_firstexon) Value Type Histone 1 0.06 high H3K27ac 2 0.12 low H3K27ac 4 0.04 high H3K27me3 5 0.16 low H3K27me3 7 0.02 high H3K36me3 8 0.13 low H3K36me3 I have another data frame with p-v alues

[R] Visualizing and clustering one half of a symmetric matrix

Hi all, I have a distance matrix (symmetric) which looks somewhat like this (only a small portion shown) ENSG0101413 ENSG0176884 ENSG0185532 ENSG0106829 ENSG0101413 1.000 1.000 1.000 1.000 ENSG0176884 0.328

Re: [R] Visualizing and clustering one half of a symmetric matrix

I do want to cluster it and only plot the lower half of the matrix. From: Peter Langfelder Sent: Thursday, September 15, 2016 11:33:13 PM To: Khan, Saad M. (MU-Student) Cc: r-help@R-project.org Subject: Re: [R] Visualizing and clustering one half of a symmetric

[R] error in data.farme--duplicate row.names error

I'm using the "GameTheory" package--- DefineGame(14,values) and values is equal to 16,383 observations. I keep getting the following error- [Error in data.frame(rep(0, 2^n - 1), row.names = Bmat) : duplicate row.names: 1, 11, 111, 12, 112, 1112, 2, 13, 113, 1113, 3, 1213, 11213, 111213

[R] (no subject)

What is plot.new? How can I fix this data or add plot.new so it works? > library(maps) > library(splancs) > area = 6*4 > lambda = 1.5 > N = rpois(1,lambda*area) > u = runif(N,-2,4) > v = runif(N,0,4) > plot(u,v,asp=1) > h = chull(u,v) > h = c(h,h[1]) > plot(u[h],v[h],"1",asp=1) Error in plot.xy(

[R] (no subject)

What is plot.new? and how can i get it to work so i can load other data? > library(splancs) > area = 6*4 > lambda = 1.5 > N = rpois(1,lambda*area) > u = runif(N,-2,4) > v = runif(N,0,4) > plot(u,v,asp=1) > h = chull(u,v) > h = c(h,h[1]) > plot(u[h],v[h],"1",asp=1) Error in plot.xy(xy, type, ...) :

[R] Moran's I test- Ordinal Logistic Regression Model

Hi, I am trying to do a Moran's I test on an ordinal logistic regression model. I have a simple spatial weights matrix listed below I would like to use. Y= 10 0 0 0 0 0 0 0 01 1 0 0 0 0 1 1

[R] Optimizing with constraints using alabama

variable available to be examined by my objective function. Currently, I am calling hin again within my objective function - certainly not an optimum solution. Thanks, Roland DePratti Graduate Student Central Connecticut State University Data Science Program Here is my inequality constraint