Using tapply when you have different n for each treatment


Using tapply when you have different n for each treatment



I am an R beginner trying to fit my data to a non-linear regression. Specifically I want to fit the rate of insect development over different temperatures to a power function. Below is my code, modified from the examples at the bottom of this page:
https://docs.tibco.com/pub/enterprise-runtime-for-R/3.1.0/doc/html/Language_Reference/stats/selfStart.html



You should be able to paste the code below right into R to get what I get.


structure(list(Larvae_temp = c(10L, 10L, 10L, 10L, 10L, 10L),
Larvae_rate = c(0.047757234, 0.04174518, 0.04174518, 0.04174518, 0.04174518, 0.04174518)), .Names = c("Larvae_temp", "Larvae_rate"), row.names = c(NA, 6L), class = "data.frame")


SSpower<-(selfStart(~A*(Larvae_temp^B),
initial=function(mCall,data,LHS)
{xy<-sortedXyData(mCall[["Larvae_temp"]],mCall[["Larvae_rate"]],LarvaeDevelopment) ##I think the error appears in this line of code, or the one below
z<-xy[["Larvae_rate"]]
aux<-coef(lm(Larvae_temp~z,LarvaeDevelopment))
pars<-as.vector(coef(nls(Larvae_rate~A*(Larvae_temp^B))),start=list(A=aux[1],B=aux[2]),data=LarvaeDevelopment,algorithm="power")
value<-list(pars[1],pars[2])
names(value)<-mCall[c("A","B")]
value},
parameters=c("A","B")))
getInitial(Larvae_rate ~ SSpower(Larvae_temp, A, B), data=LarvaeDevelopment)



When I show traceback, these are the steps that appear:


Error in tapply(y, x, mean, na.rm = TRUE) : arguments must have same length
9. stop("arguments must have same length")
8. tapply(y, x, mean, na.rm = TRUE)
7. sortedXyData.default(mCall[["Larvae_temp"]], mCall[["Larvae_rate"]],
LarvaeDevelopment)
6. sortedXyData(mCall[["Larvae_temp"]], mCall[["Larvae_rate"]],
LarvaeDevelopment)
5. (attr(object, "initial"))(mCall = mCall, data = data, LHS = LHS)
4. getInitial.selfStart(func, data, mCall = as.list(match.call(func,
call = object[[3L]])), LHS = object[[2L]], ...)
3. getInitial(func, data, mCall = as.list(match.call(func, call = object[[3L]])),
LHS = object[[2L]], ...)
2. getInitial.formula(Larvae_rate ~ SSpower(Larvae_temp, A, B),
data = LarvaeDevelopment)
1. getInitial(Larvae_rate ~ SSpower(Larvae_temp, A, B), data = LarvaeDevelopment)



I think the problem might come from the fact that the n of each treatment is unequal. For example, n= 58,165,113,26 for temperatures 10, 15, 20, 30. The corresponding y-values (Larvae_rate) match each x-values (Larvae_temp). Does anybody know how to fix this?



I'd also really appreciate it if someone could look over my code and see if it makes sense! I am really new to R still and debugging this code has taken over my life.



Please let me know if I can provide more information. Thanks!!









By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

-

Popular posts from this blog

How to input without newline? (Python)

How to input without newline? (Python) I am trying to create a simple file editor in Python via the console (called PyCons). Here is my current program: def PyCons(f): file = open(f, "r") appe = open(f, "a") print("n=====nPyCons Editor, v1.2nPython 3.6.1n") global i i = 0 for line in file: i += 1 print(" {}| ".format(i) + line, end='') for k in range(10000000000000): print("n {}| ".format(i+1), end='') inp = input("") if inp == "@PyCons.save()": print("n=====") break else: i += 1 appe.write("n" + inp) I used the end method to make sure that the existing code in the file was printed properly in its line format, but the input function does not have an end attribute. So, when I go to enter code in the editor: end input end PyCons("magic.html") ... ===== PyCons Editor, v1.2 Python 3.6.1 1| <p>Hello!...
Read more

C++ thread error: no type named ‘type’ MINGW

C++ thread error: no type named ‘type’ MINGW I'm currently doing a project for my operational systems class and I have to make a program to find prime numbers to run in threads. So I did this: #include <iostream> #include <cmath> #include <thread> #define THREADNUMBER 100 using namespace std; int CONTADOR = 0; int CONTADORTHREADEXECUTADA = 0; int primeRange(int v1, int v2) { int a, limit; bool isprime; for (int i = v1; i <= v2; i++) { if (i % 2 != 0 && i % 3 != 0) { isprime = true; limit = i / 2; if(i == 1) isprime =false; limit = (int)sqrt(i); //General case for(a=2; a <= limit; a++){ if(i % i == 0 && i != 2){ isprime = false; break; } } if (isprime) { CONTADOR++; } } } CONTADORTHREADEXECUTADA++; return 1; } int main(int argc, char *argv[ ] ) { int number1 = atoi(argv[1]); int num...
Read more

Analog for TagView in flutter

Image
Analog for TagView in flutter I want to have UI like this: With native Android it can be done with libraries like https://github.com/whilu/AndroidTagView, how can it be done with flutter? 2 Answers 2 you can make the same UI by combining the Wrap and Chip widget as @wasyl montioned . but this is a full example about what you need Wrap Chip Notes : spacing deleteIcon avatar Chip Shape new Column( crossAxisAlignment: CrossAxisAlignment.stretch, children: <Widget>[ new Wrap( children: <Widget>[ Chip( label: new Text("Java"), onDeleted: () {}, labelPadding: EdgeInsets.all(2.0), deleteIcon: Icon(Icons.clear), ), Chip( label: new Text("C++"), onDeleted: () {}, labelPadding: EdgeInse...
Read more