Thanks. Dat verhielp het probleem inderdaad.quote:Op vrijdag 21 juli 2017 17:44 schreef Z het volgende:
De boel even opnieuw opstarten wellicht? Ik weet niks van SAS maar met Macro's weet je soms niet welke er hoe in het geheugen staat.
Optie 1 valt af omdat je een aanname schendt, namelijk die van onafhankelijke waarnemingen. Dat geldt ook voor optie 2. Repeated measures anova houdt hier wel rekening mee, dus dat lijkt me de voorkeur hebben.quote:Op woensdag 26 juli 2017 16:13 schreef Sonyanijntje het volgende:
Hallo!
Ik hoop dat jullie mij kunnen helpen! Ik heb data van mijn onderzoek maar ik twijfel over statistische test die ik moet gebruiken.
[...]
Iemand een suggestie welke van mijn opties ik het beste kan kiezen?
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | %combineddata(y); Mergeddata&y; merge a b; run; %if &y<2003 %then %do; data test&y; merge mergeddata&y d; by db32; run; %end; %if 2003 le &y le 2007 %then %do; data test&y; merge mergeddata&y f; by db45; run; %end; |
Kun je ze niet los van elkaar runnen en die van voor 2003 net zo schrijven als 2003-2007?quote:Op woensdag 2 augustus 2017 12:16 schreef Mynheer007 het volgende:
Ik heb een sas-macro maar er gaat iets niet helemaal naar behoren
[ code verwijderd ]
Ik wil dus verschillende acties uitvoeren voor de jaren voor 2003 en de jaren 2003 tot en met 2007.
Bij het runnen van de macro voor de jaren 2003-2007 gaat alles naar behoren, alleen bij het runnen van de jaren voor 2003 krijg ik een error:
Daaruit blijkt dat de macro op de data van de jaren kleiner dan 2003 ook de actie voor de jaren 2003-2007 uitvoert. Ik zie alleen niet waar de fout in mijn script zit
Ik heb hem herschrevenquote:Op woensdag 2 augustus 2017 12:22 schreef MCH het volgende:
[..]
Kun je ze niet los van elkaar runnen en die van voor 2003 net zo schrijven als 2003-2007?
Ik heb de volgende regel herschreven:quote:Op woensdag 2 augustus 2017 12:22 schreef MCH het volgende:
[..]
Kun je ze niet los van elkaar runnen en die van voor 2003 net zo schrijven als 2003-2007?
1 | %if 2003 le &y le 2007 %then %do; |
1 | %if 2003 le &y AND &y le 2007 %then %do; |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 | makelms <- function(){ # Store the coefficient of linear models with different independent variables cf <- c(coef(lm(Fertility ~ Agriculture, swiss))[2], coef(lm(Fertility ~ Agriculture + Catholic,swiss))[2], coef(lm(Fertility ~ Agriculture + Catholic + Education,swiss))[2], coef(lm(Fertility ~ Agriculture + Catholic + Education + Examination,swiss))[2], coef(lm(Fertility ~ Agriculture + Catholic + Education + Examination +Infant.Mortality, swiss))[2]) print(cf) } # Regressor generation process 1. rgp1 <- function(){ print("Processing. Please wait.") # number of samples per simulation n <- 100 # number of simulations nosim <- 1000 # set seed for reproducability set.seed(4321) # Point A: x1 <- rnorm(n) x2 <- rnorm(n) x3 <- rnorm(n) # Point B: betas <- sapply(1 : nosim, function(i)makelms(x1, x2, x3)) round(apply(betas, 1, var), 5) } # Regressor generation process 2. rgp2 <- function(){ print("Processing. Please wait.") # number of samples per simulation n <- 100 # number of simulations nosim <- 1000 # set seed for reproducability set.seed(4321) # Point C: x1 <- rnorm(n) x2 <- x1/sqrt(2) + rnorm(n) /sqrt(2) x3 <- x1 * 0.95 + rnorm(n) * sqrt(1 - 0.95^2) # Point D: betas <- sapply(1 : nosim, function(i)makelms(x1, x2, x3)) round(apply(betas, 1, var), 5) } betas |
1 2 3 | dataset2000<-read_sas("dataset2000.sas7dbat") plot(dataset2000, x,y) etc.... |
quote:Op zaterdag 26 augustus 2017 15:38 schreef Z het volgende:
Moet kunnen. "Putting all the data frames in a list and looping over that list with lapply".
Bedankt voor de hulp, maar dit gaat specifiek over de dataset. Maar ik wil ook dat de functie de titel van de grafiek aanpast etc., maar aangezien die tussen "" staat, doet de functie daar niets mee.quote:Op zaterdag 26 augustus 2017 22:22 schreef ralfie het volgende:
ehh, paste0("dataset", n, ".sas7dbat") ?
Mann-whitneytoets.quote:Op woensdag 30 augustus 2017 16:49 schreef Lyrebird het volgende:
Populatie a: 500.000 samples, mean = 19, stdv = 10
Populatie b: 500 samples, mean = 23, stdv = 11
Populaties zijn niet normaal verdeeld.
Wat voor een test kan ik het beste gebruiken om aan te tonen of populatie b binnen/buiten populatie a valt?
Nee, want je schendt de assumptie van onafhankelijke waarnemingen, wanneer er sprake is van meerdere metingen per persoon. Een waarneming is dan niet meer onafhankelijk want het is immers afhankelijk van de persoon.quote:Op donderdag 21 september 2017 19:23 schreef No-P het volgende:
Graag uw hulp.
Ik heb van opeenvolgende jaren een populatie gevolgt waarvan een proportie positief scored op een test (dichotoom: positief of negatief).
Graag wil ik een trend analyse doen om aan te tonen dat de proportie positieve testen toeneemt in de tijd. De gescreende individuen in de populatie die worden getest verschillen per jaar.
Welke methode is hier geschikt voor? Ik denk dat ik met logistic regression een heel eind kom.
Wat wil je uberhaupt onderzoeken?quote:Op vrijdag 13 oktober 2017 21:21 schreef icecreamfarmer_NL het volgende:
Meer statistiek/onderzoek maar het heeft een wiskunde component en het is niet de moeite om een apart topic te openen.
Ik zit met het probleem dat ik niet weet of ik een t-test moet doen of een regressie. In het verleden heb ik het allemaal gehad maar het is weg gezakt.
In het kort het onderzoek:
Op dag 1 wordt gevraagd naar de mening over A. (Via een likert schaal).
Op dag 2 wordt onder een compleet andere groep mensen gevraagd naar de mening over B.
Het is trouwens onbekend of de variatie van beiden gelijk zijn.
Nu is de vraag moet ik dit onderzoeken met een double tail independent two sample T test of Welchers T test. Of dat ik het beter kan doen met een regressie (least squares).
Zo ja welke moet ik kiezen en vooral waarom. Mijn gevoel en volgens wiki zegt de T test echter kom ik niet echt achter de voordelen van een T-test over een regressie.
SPOILEROm spoilers te kunnen lezen moet je zijn ingelogd. Je moet je daarvoor eerst gratis Registreren. Ook kun je spoilers niet lezen als je een ban hebt.
[ Bericht 57% gewijzigd door peperkoekmannetje op 30-10-2017 14:16:42 (toevoeging) ]
Forum Opties | |
---|---|
Forumhop: | |
Hop naar: |